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1. Introduction and Summary 


1. 1 INTRODUCTION 

The geosynchronous microwave atmospheric sounding radiometer (MASR) concept 
has great merit in meteorological applications. This study outlines an instrument 
design which meets needed system performance goals, and shows that the instru- 
m e nt technology is feasible. 

The microwave atmospheric sounding radiometer is a four-band, millimeter wave 
radiometer designed to accomplish three-dimensional atmospheric temperature 
and relative humidity mapping from geosynchronous orbit. The MASR mission 
collects data to aid in the observation and prediction of severe storms. The geo- 
synchronous orbit allows the continuous atmospheric measurement needed to 
resolve mesoscale dynamics. The instrument may operate in conjunction with 
muitispectral imaging and infrared sounding sensors, and will be particularly 
valuable in collecting data in overcast regions. The anticipated ability to predict 
severe thunderstorms and larger cyclonic storms from MASR data adds to the 
economic and social benefit. 

Angular resolution of MASR is limited by the beam width of a shuttle launched 
antenna (4.4 meter diameter) at the 118 and 183 GHz oxygen and water vapor sound- 
ing frequencies. Temperature resolution is a function of system noise temperature 
and integration time. Resolution element size and required integration time 
combine to limit the mapping coverage rate. From geosynchronous orbit, the 
4. 4 meter antenna provides nadir earth surface resolution of about 20 km at 
183 GHz. The antenna is mechanically scanned in raster fashion to cover an area 
with contiguous resolution cells. The nominal dwell time per cell is 1 second at 
183 GHz, allowing an area 800 by 800 km to be mapped in about 30 minutes. 

Temperature sounding is achieved by obtaining a profile of the thermal radiance 
of the atmosphere in the frequency region of the 118 GHz oxygen absorption line. 
There are 11 double sideband (DSB) channels centered at the peak of this line and 
located various distances from the line center. Each channel will receive radia- 
tion predominately from a given altitude layer of the atmosphere, and the atmo- 
spheric temperature profile can be obtained from the 11 measured brightness 
temperatures using an inversion algorithm. 

Similarly, there are six DSB channels centered at the 183 GHz water vapor line 
to determine the humidity profile. In addition, there are two window bands at 
104 and 140 GHz to correct the measurements at 1 1 8 and 183 GHz for ground 
emissions (see Figures 1 and 2). 

A 4.4 meter symmetric Cassegrain antenna brings all four frequency bands to a 
common focus. The radiometer receiver package is designed to demultiplex the 
four bands through the use of quasi-optical filters and to direct the individual 
beams to four separate feed horns, followed by four separate receivers. The 
118 and 183 GHz receivers have very wide band channelized IF systems for the 
purpose of profiling the microwave spectrum of the atmosphere at these frequen- 
cies. Window bands measure the temperature of the earth's surface and there- 
fore, each requires only one IF channel. A common chopper wheel and calibration 
load* serve all four bands. Synchronous detection and other signal processing are 
achieved through a single onboard digital processor. Figure 3 is a simplified 
block diagram of the radiometer receiver subsystem. 
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Millimeter wave technology status offers several approaches to MASR, each having 
"softness” but currently under rapid development. A preferred technical approac 1 
is difficult to recommend at present, but several technical approaches can meet 
the MASR system requirements. The present study effort has selected a tentative 
configuration which appears to have least technical risk. Technology advances in 
alternate approaches could change the selection but the present choice serves as 
a model to estimate weight, power, size, and cost of the flight instrument. 


This report addresses the tradeoff and design study of the MASR radiometer. 
Volume III gives details of the 4.4 meter antenna subsystem. 


The authors are grateful to the following for their assistance in preparing this 
final report: D. H. Staelin, Massachusetts Institute of Technology D. Held 

and A R. Kerr, NASA-Goddard Institute of Space Sciences. T. S. Chu »\eLl 
To Ieoh one Laboratory; J. M. Schuchardt, J. R. Langley, and J. A. Stratigos, 
Ceoraia Institute of Technology; J. Payne, B. L. t’lich. National Radio Astronomy 
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1 . Introduction and Summary 
1.2 SUMMARY 


Design analyses and highlights of technical tradeoffs surrounding a radiometer 
coniiguration are described in this summary. Estimated AT's will meet antici 
pa ted MASR needs, 


MASR Study objectives are to make technical tradeoffs, determine feasibility, and 
isolate key technology areas. These objectives were achieved with in-depth tech- 
nology surveys, concluding with a receiver hardware definition of sufficient detail 
to allow weight and power estimates and the preparation of a preliminary 
specification. J 

The study revealed the importance of quasi-optical techniques to reduce feed 
diplexer circuit, and local oscillator (LO) filter losses. It also indentified the 
mixer as the "softest" technology, including the mixer configuration, mixer diodes 
if "P e dance match, mixer bandwidth. and mixer local oscillator. 

While the best state ol the art mixer performance was reported to be sufficient 
to meet MASR requirements, typical performance is in general marginal. Con- 
clusions are that MASR is feasible, AT (radiometer temperature resolution) 
requirements generally can be achieved, quasi-optical design techniques are 
recommended, and mixer development me hodologv will pace the program 
program. “ * 


This technical report consists of sections on Design Analysis 
and Antenna Interface; appendices deal with the MASR concept 
formalism, reference averaging, and properties of materials/ 
paragraphs summarize the highlights of each section. 


Receiver Design, 
and optical design 
The following 


D e sign Anal ysis. Radiometer temperature resolution is commonly expr 
the minimum temperature difference which can he detected 


essed as 


AT KT s //I7 


where 


T s v system noise temperature 
B - IF bandwidth 
T - integration time 

L r nTj Sy f e , m f na ^ SiS ? f Ule '' adiometer - was determined that the effective 

tarn) I! !! ^ °J l' 6 raC ' i0meter (normall y K a for a Dicke switched sys- 

tem) could be reduced to a value approaching 1.4, depending on the gain stability 

of the system, by using the technique of reference averaging. The improvement 
in temperature resolution results by integrating the reference channeled- a 

™j', er , of pJ: ; el w- rl ,° dS ’ Wher0 a pixul purioci equals the integration time of the 
signal channel. With practical limitations, improvement factor offered by refer- 
ence averaging is about equivalent to that achievable with a double Dicke (Graham's 
receiver) system, but does not require the complexity of two receivers! 


1-6 






WV 1 ... 


It is recommended that calibration of the radiometer receiver be performed at the 
antenna feed with a movable reflector that periodically switch c* the receiver input 
to a reference hot load and to cold space. The calibration accuracy expected is 
of the order of ±1°K in ail channels. The radiometer system calibration, including 
the main antenna, must be performed in the far field of the main antenna. Pos- 
sible methods of achieving total radiometer system calibration are discussed in 
Topic 5. 4. 

It is recommended that quasi-optical filters be used both for LO filters and for 
frequency multiplexing in MASR. Analytical techniques commonly used for optical 
system design were found to be advantageous; ray matrices and Gaussian beam 
formalism provide an extremely versatile design approach for quasi -optics . 
Unconditional stability of the beam optical system assures negligible loss due to 
diffraction and radiation from these open structures. Appendix B provides the 
fundamental reference work of Kogelnik and Li which details the mathematical 
formalism; Topic 3.7 summarizes the work and illustrates its applications to 
quasi-optical systems. 

Appendix C presents an analysis that determines the effects of gain variation 
errors in a reference averaging radiometer and the maximum improvement 
achievable. This analysis shows that for ypical gain stability the K value can be 
reduced from K = l (Dicke case with square wave modulation and demodulation) 
to K = 1.4. The analysis is summarized and results presented in Topic 3.4. 

A theoretical treatment of reference averaging was conducted which includes 
general considerations of chopping frequency, asymmetrical chopping, and 
chopper transition effects. Because this work is highly theoretical and lengthy, 
it is included in Volume IV of this report. 

Another theoretical treatment was developed of chopper frequency selection in a 
reference averaging Dicke radiometer. One of the more significant results is 
that second detector low frequency noise (1/f), which is suppressed in a conven- 
tional Dicke radiometer by proper selection of chopper frequency, cannot be sup- 
pressed in the same way in a reference averaged system. Rather, the reference 
averaged system ha s a floor on the detector noise feedthrough. Becaus e this 
analysis is highly theoretical and lengthy, it is also included in Volume IV. 

Receiver Design. A tentative receiver subsystem block diagram is shown in 
Figure 1. The basic four-band radiometer receiver consists of a quasi-optical 
quadraplexer which separates the bands and directs the beams to separate feed 
horns and mixers. The 118 and 183 GHz mixers are matched into IF networks 
that provide IF channel separation for the 11 temperature channels and six water 
vapor channels. The window bands at 104 and 140 GHz are amplified after mixing 
over nominal IF bandwidths of about 1 GHz and detected. A signal processing 
and control subsystem manages all functions, including the antenna drive, calibra- 
tion control, chopper synchronization, and synchronous detection of the 19 radi- 
ometer channels. In addition, the signal processing and control subsystem 
provides the capability to perform reference averaging to enhance the performance 
of the radiometer. 

Radiometer temperature resolution or AT performance is a direct function of sys- 
tem noise temperature. To facilitate the following discussion, the contributing 
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factors are given in Table 1, together with example values for parameters for 
each radiometer band and a computation of its respective system noise 
temperature. 

Mixer Tradeoff. The heart of the radiometer receiver is the low noise mixer. 
Reported mixer conversion loss and noise temperature data are . scattered over 
a wide range of values for the frequencies of interest. It is obvious that con- 
sistent mixer performance will be achieved only with great effort. A recom- 
mended approach is to begin mixer development early and to build a selected 
inventory of mixers suitable for MASR. 

The quality of mixer performance is conventionally expressed in terms of 
single -sideband noise temperature TmSSB atlc * mixer conversion loss Lm- 
Mixer performance is affected by the quality of the mixer diode material, 
precision contacting, microwave circuit design, and impedance match for 
both the signal input and IF output circuits. Further, it is also a function of 
the quality and level of local oscillator signal. Of all the factors influencing 
the radiometer system noise temperature, TjyiSSB is the term over which 
there is least control. State of the art values of TmSSB reported in the liter- 
ature vary over a wide range. Figure Z shows a number of demonstrated 
values. The solid line shows the probably minimum achievable values 
expected, and also represents approximate MASR needs for the more difficult 

bands . 

Balanced mixers have the desirable property of canceling LO noise and are 
recommended for MASR to provide a degree of immunity to the noise generated 
by solid state LO sources such as IMPATTS. Waveguide balanced mixers are the 
recommended baseline design. However, conventional balanced mixers at these 
frequencies are difficult to build. Subharmonically pumped mixers a form of 
balanced mixer, also have LO noise cancellation properties. Further, these 
devices have the advantage of using a local oscillator at half the signal frequency. 

Ouasi-Ootical Balanced Mixer. This study has produced an alternate approach to 
ac hieving the noise cancellatio n effects of a balanced mixer with two single-ended 
mixers in a quasi-optical assembly (Figure 3). The signal from the antenna is 
diplexed with a wire grid polarizer 1 and the i83 GHz band is reflected by 90 . 

The signal polarization is upward (not shown in figure) while the local oscillator 
polarization is in the plane (see Figure 3). Both signal and Local oscillator 
encounter a rotation plane of 45° and are split by a second wire grid polarizer Z 
and directed to separate single-ended mixer mounts. The signal and LO lields 
are in phase in one mixer and out of phase in the other. Their electrical outputs 
combine in such a way as to cancel the dc and noise terms while the IF terms add 
in phase. This approach was determined to be unique and was reported as a new 

technology item. 

Local Oscillator Sources. It was determined that IMPaTT and Gunn solid state 
sources are available for use as local oscillators. IMPATTS produce sufficient 
power in all four bands but are noisy and must be filtered with about 40 dB of 
filtering over each band of frequencies, which for the 183 GHz humidity band is 
±10 GHz from the center frequency. The basic approach to accommodating t e 
IMPATT is to utilize balanced mixers, that provide about 20 dB of immunity to 
the noise, and Fabry Perot quasi-optical filters that provide at least another 25 dB 
of filtering. The combined effects of the balanced mixer and quasi-optical filter 
provide sufficient margin to assure that a good low noise mixer is achieved. 
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TABLE 1. CONTRIBUTORS TO To YSTEIVI 


T SYSTEM “ T A + (L F ‘ T 0 + L F t r 


T R - 2 (T MSSB + L M T IF ) 


(Lc-1) T 


T r = double sideband receiver temperature 

t a = antenna temperature 

L r = feed loss 

Tq = feed temperature 

L|^ “ mixer conversion loss 

T|p = IF noise temperature. 

F|p = IF noise figure 
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I'he Gunn oscillator (or the transferred electron oscillator) is somewhat quieter 
than an IMPATT, but has an upper frequency limit near 100 GHz. With noise 
characteristics lower than that of an IMPATT and adequate power below 100 GHz, 
the Gunn oscillator is an attractive alternative to the IMPATT when used with a 
doubler or with a subharmonic mixer. 

Local Oscillator Filters. Because of relatively high unloaded Q, Fabry-Perot 
quasi -optical filters have fundamental advantages over waveguide filters for 
selectivity and suppression of LO noise in the IF sidebands. Of the two typos 
demonstrated, the parallel plate Fabry-Perot has the advantage of arbitrarily 
large rejection bands, while the ring configuration offers a unique LO injection 
scheme. Because of the wide IF bandwidth requirements however, the recom- 
mended filter for MASK local oscillators is the parallel plate configuration. 
Specific parallel plate filter designs were developed for each of the MASR bands. 
The most difficult filter is that for the 1 18 GHz band which needs a loaded O of 
6000 to accommodate a minimum IF frequency of only 20 MHz. For a 118 GHz 
minimum IF frequency of SO MHz instead of 20 MHz, a loaded Q of 2400 would be 
adequate, resulting in LO filter losses of 4. 1 dB instead of 8. 1 dB. The l 18 ( Hz 
LO power is marginal with 8. 1 dB filter loss, and raising the minimum 118 GHz 
IF would eliminate this problem. 

IF Design . Difficult requirements for MASR are the 0.6 to 10 GHz IF network 
for the Il^O band and the 20 to 4000 MHz IF network for the band. Two 
approaches were analyzed; one utilizes available components and techniques while 
the other assumes the development of contiguous filter banks. The analyses show 
that while both approaches are feasible and technically suitable, the conventional 
approach requires less development. In either case, the ultimate performance of 
the IF network will be determined by how well the mixer-to-IF impedance match 
is maintained over the required bandwidth. 

Antenna Interface. The antenna interface consists of an RF Interface, thermal/ 
mechanical interface, and calibration interface. In a typical antenna, the RF 
interface is defined by a waveguide that connects the receiver to the antenna feed. 
In the recommended quasi-optical design, the RF interface is defined by the com- 
plex beam parameters for the four bands. The illumination patterns for the 
individual bands are generated separately and combined in the quasi-optical 
quadraplexer . 

The calibration Interface deserves special comment. The calibration procedure 
for the radiometer receiver is self-contained. A positionable mirror at the 
antenna feed allows the receiver to look at the antenna, cold space, or a calibrated 
reference "hot" load. This two-point receiver calibration sequence could he per- 
formed in a few seconds and repeated often, allowing the receiver to measure 
antenna temperature T*\ to an accuracy of about 1'K. 

Knowing T^ accurately is only on e step in obtaining mesoscale temperature data 
since the value of T,\ is a function of the antenna beam efficiency, sidelobes, field 
brightness, and ohmic losses. The ultimate calibration of the actual target tem- 
perature against l ^ is a ditticult task which must he carefully performed in space. 
Topic 5.4 is devoted to the discussion of antenna calibration. 

^Iblp. -Qp^al Antenna Feed. The tradeoff study shows that feed, diplexing, and 
LO filter Tosses could be significantly reduced by using quasi-optical components 
rather than waveguide components. Typically, a ten-fold reduction in losses is 
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each horn, are given in Figure 4. 
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1, Introduction and Summary 

1.3 RECOMMENDED RADIOMETER CONFIGURATION 

The present study effort has resulted in the selection of an overall configuration 
which offers excellent theoretical performance and affords the least t echnical risk. 

^^Z^l teChn0,08y S ‘ atUS ° ffers severaI a PP roac hes to MASR, each having 
aoDroach H f currently under rapid development. A preferred technical 

me P nts The P lesent - b “ t several configurations can meet MASR require- 

ments. The study effort has resulted in the selection of a tentative configuration 
which appears to have least technical risk. Technoiogy advances in alternate 
ppioaches could change the selection but the present design serves as a model to 
estimate weight, power, size, and cost of the flight instrument. 

The baseline approach utilizes balanced mixers in all bands. IMPATT local 
oscillators with quasi-opttcal filters are recommended for 118, 140, and 183 GHz 

Quasi' o'^c^ a H °H ?S | CllUt0r/d0Ub,er U recomme *ded the 104 GHz band. 

P * f!i de Vo U J tl ?!.« Xer8 are recommended in the form of band dropping 

wire arid i* 1 ?- 14 °’ aad 1 1 bancIs * These dropping mirrors are either 

gud polarizers, resonant grid dichroic mirrors, or Fabry-Perot filters An 
isometric view of the optomechanical assembly is shown in Figure 1. The incom- 
reflert™ the * nteana reaches a minimum beam waist at the calibration sub- 

rmatched n off r av; ^ f ° CU8 ™ d a ** T P'Y to the chopper plane and recoUimated by 
ff " aX1S mirPor * T l ie beam is formed and directed through a series of 
... ( rmirors and to the balanced mixers for each band The local 

reapecth-e *IF "n e t work sT m b 1 n 6 ^ with the si Snat output, are directed the 



































2. RECOMMENDED RADIOMETER CONFIGURATION 

2. 1 PERFORMANCE ESTIMATES 

2.2 RADIOMETER MECHANICAL SUBSYSTEM 

2. 3 MASR WEIGHT AND POWER REQUIREMENTS 
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Channel 

Number 


Average 


118.75 GHz 

Waveguide balanced 

1 100°K 

5dB 

IMPATT 

Fabry-Perot parallel plate 
19 to 4090 MHz 
T riplexer/9 way power divider 
1 dB 
250° K 
1 .55 sec 


Center frequency 
Mixer 

Mixer temperature, T^ggg 
Mixer conversion ioss, L m 
Local oscillator 
Local oscillator filter 
Total IF bandwidth 
IF filter 
Feed loss 

Antenna temperature 
Integration time 


2. Recommended Radiometer Configuration 
2. 1 PERFORMANCE ESTIMATES 

Tables 1 through 3 are performance estimates for the four radiometric bands 
based on state of the art technologie s with modest improvements anticipated. 

TABLE 1. PERFORMANCE ESTIMATES FOR 113 GHz BAND 





AT 

§ 

u 

Goal 

Dicke (K = 2) Reference A /erage (K = 1.4) j- 

0.2 

0.06 

0.04 

0.2 

0.08 

0.06 

0.2 

0.11 

0.08 

0.2 

0.14 

0.10 

0.2 

0.16 

0.11 

0.2 

0.18 

0.13 

0.2 

0.22 

0.16 

0.2 

0.27 

0.19 

0.2 

C.34 

0.24 

0.2 

0.40 

0.29 

0.2 

0.40 

0.29 

0.2 

0.21 

0.15 


. 



u 

l 

i. 


I 

5 

i 

i. 


»■ 



TABLE 2. PERFORMANCE ESTIMATES FOR 183 GHz BAND 


Channel 

Number 

f|F* 

GHz 

B|F. 

GHz 

f IF, 

dB 

T|F, 

°K 

AT 

Goal 

Dicke (K = 2) 

Reference Average (K - 1.4) 

H*1 

8.625 

2.5 

5.5 

738 

1 ' 

0.18 

0.13 

H-2 

6.25 

2.5 

5.1 

648 

US 

0.16 

0.11 

H-3 

4.375 

1.25 

5.3 

692 

0.2 

0.24 

0.17 

H-4 

3.125 

1.25 

4.9 

606 

0.2 

0.22 

0.16 

H>5 

1.875 

1.25 

4.3 

490 

0.2 

0.19 

0.14 

H-6 

0.9375 

0.625 

3.7 

389 

0.2 

0.24 

0.17 

Average 

0.2 

0.21 

0.15 


Center frequency 


183.30 GHz 


Mixer 

Mixer temperature, Tj^ggg 

Mixer conversion loss, I 

Local oscillator 


Waveguide balanced 
2100°K 
8.5 dB 
IMPATT 


Local oscillator filter 
IF bandwidth 
I F filter network 
Feed loss, Lp 

Antenna temperature, T A 
Integration time 


Fabry-Perot parallel plate 
0.625- 10 GHz 
Cascaded diplexers 
0.5 dB 
250°K 
1 .0 sec 


xe-zeso* 













TABLE 3. PERFORMANCE ESTIMA TES FOR 104 AND 140 GHz BANDS 
Z Ft, ' - AT 


channel 

f|F. 

B|p» 

F | p» 

Tip 


Mumber 

MHz 

MHz 

dB 

OK 

Goal 

Dicke (K = 2) 

W-1 

500 

1000 

3.0 

290 

0.2 

0.05 

W-2 

500 

1000 

3.0 

290 

0.2 

0.08 


W-1 center frequency 
W-2 center frequency 
Mixers 

Mixer temperature, T MSSB 
W-1 
W-2 

Mixer conversion loss, L m 
W-1 
W-2 

Local oscillator 
W-1 
W-2 

Local oscillator filter 
W-1 
W-2 

Total IF bandwidth, each band 
Feed loss, each band 
Antenna temperature, T A 
Integration time 
W-1 
W-2 


104 GHz 
140 GHz 

Wavequidp halanced 

600°K 

600°K 

5dB 
6.5 dB 

Gunn doubler 
IMPATT 

Waveguide 

Fabry-Perot parallel plate 
1 GHz 
0.95 dB 
250°K 

1.76 sec 
1.31 sec 
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2. Recommended Radiometer Configuration 
2. 2 RADIOMETER MECHANICAL SUBSYSTEM 
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2. Recommended Radiometer Configuration 

2. 3 MASR WEIGHT AND POWER REQUIREMENTS 

The weight, power, and thermal /mechanical interface data of the MASR payload 
are tabulated in Table 1. 


TABLE 1. PAYLOAD WEIGHT AND POWER REQUIREMENTS 


Weight Tabulation, lb 

Optomechanical subsystem 

15.5 


183 GHz receiver subsystem 

4.2 


118 GHz receiver subsystem 

3.7 


104/140 GHz receiver subsystem 

1.8 


Signal processing subsystem 

1.0 


Miscellaneous electronics 

1.4 


Cables and connectors 

4.6 


Chassis 

18.0 


Total MASR weight estimate, lb 

50.2 


Power Tabulation, W 

Chopper drive 

5 


Calibration mirror drive 

15 

0 

Local oscillator sources 

45.3 


IF amplifiers 

3.2 


Signal processing 

13.2 


Power conditioning 

16.3 

14.3 

Total power required, W 

98.0 max 

81.0 min 

Package Dimensions, in. 

16 x 14 x 32 

Thermal Load, W 

81 to 98 
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3. 2 NOISE TEMPERATURE 

3.2.1 System Noise Temperature 

3.2.2 Mixer Noise Temperature 

3.2.3 IF Contribution to Noise Temperature 
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3. Design Analysis 
3.1 MASR CONCEPT 

SSiESaS-rSaKSSS: 

ature information as a function of altitude. — — - 

sounding, the brightness \ersus ircq y . . f rorn the radiative transfer 

^3S=»:sssf 
Si a; sssf rt^SS-. 

The radiometer measures directly the antenna temperature given by 


T A< E > 


= J -I a(f.0,<S>)T (£,e,<t>)dn 
4 w 4r B 


calibrations based on the antenna gain function are generally inaccurate. 

1 i • j. , • | to 7 GHz and measurements at tins frequency are 

£* fSSSIfiiVSSS! -2 »"- |d 

about temperature as a function of altitude. 

Appendix A presents a more detailed physical analysis of the MASR concept. 
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3. Design Analysis 
3. 2 Noise Temperature 

3 2 1 SYSTEM NOISE TEMPERATURE 

sistent with convention, is pi esented. - ■ 


The receiver noise temperature T R 1 3 ^ b ^he are^fo generate the noise 

resistance of an otherwise noiseless 1 f ce . connected to an idealized noise- 

less^ ou re' e .'"it 6 die r e^ore^ r ep^r e s e nts °th e noise generated within the receiver 
referred to its input port. 

Th, system «»,— Js th. 

'[“bm.'S?-K'‘T R 1 *.”eS P e..d o, ie.tntsetm.. *. mi*.r s.d 

IF, and for a single sideband system 


RSSB 


= t mssb + l m t if 


T MSSB ^^T F Se f feV a encTe U s r ^ ™ ise , 

1. muttipltea b, >b. miser »»»" 


temp 

Lm- 


to a single sideband receiver system (see Figure 1) the imag e^eband 

is reiected, all sources of internally generated teceivei 
the mixer signal input and the system temperature 


T SSSB = t a + t rssb = t a + t mssb + l m t if 


Physically, image rejection can be implemented using a sharp RF filter at the 
mixer input. 

"seftsssp 

“ the^ame « in the single sideband case, but now its effect is divided between 
two signaL channels with the result that 


T R = I t rssb 


T = T. + 1/2 (T MSS b + h M T 1F ) 



Feed losses between the antenna and the mixer raise the system noise temperature 
by attenuating the signal and by emitting noise. For equal signal and image chan- 
nel losses L F , the system temperature becomes! 2 ) 


F 


T s = T A = (L F ■ 11 T o + 2 (t mssb + l m t if ) 


where T Q is the physical temperature of the lossy element. For T Q = 290 ®K, 
about 7 °K noise is emitted for each 0. 1 dB loss. The factor Lt F also multiplies 
the receiver noise temperature, which results in an added contribution for the 
highest frequency channels of 


AT S = *“^~2 (T MSSB + h U T IF } 4 


1. 023 - 1 


(4000 °K) a 47 °K 


Thus, each 0. 1 dB loss between the antenna and mixer adds approximately 54 °K 
to the system noise temperature. 


References 


1. J.D. Kraus, "Radio Astronomy, " McGraw Hill, New York, 1966, pp 265-266. 

2. A.J. Giger. et. al. , BST.T, 42, No. 4, P1096, July 1963. 
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IMAGE CHANNEL 


double-channel 

RADIOMETER 

T s = t a + i 

(t mssb + l m t if ) 
t im = t a 
l im d l m 




Lfc =L fs 
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t A - ANTENNA TEMPERATURE FOR IMAGE CHANNEL 
1M 

FIGURE 1. MIXER MODEL. EQUIVALENT CIRCUITS 
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3. 2 Noise Temperature 

3, 2. 2 toER 

A millimeter wave diode mixer can be characterized by its conversion loss Lu 

i “P ut “ ise temperature T MSSB . Experiment shows that simple 
mixer models relating the properties of a diode mixer to those of passive devices 
have limited usefulness. ^ 

Neglecting feed losses, the system temperature of a double sideband receiver was 
shown to be 


T c = T. +i(T 


*S ' 2 (J -MSSB + L M T IF ) 

t MSSB adds directly to system temperature while multiplies the IF noise tern- 

perature Tjy, In the high, frequency MASK channels these two contributions to T„ 
are of comparable size. 5 

TMSSB can be determined using the relationship for receiver noise temperature 
T -R = 2 fT MSSB + L M T IF ) 

and measuring Tr, L^, and Tip. Rapid and accurate techniques for measuring 

mixer performance have been developed. (1) measmmg 

^ratur^ritif^^effnef ti0 « ° f Performance is the mixer output noise tern- 

1 1 4 t( defmed as , th f ratl ° of mixer output noise temperature to its 
physical temperature when both the signal and image channels are terminated at 
X Q . V When = Lj M 


T out = tT o =Lf/if; + 

M M 


MSSB 

:u. 


MSSB 


= (L*, t - 2) T. 


similar concept is that of diode noise temperature T n , defined as the physical 
temperature of an attenuator having insertion loss and at the signal and 
image frequencies and producing the same output temperatures the mixer* 2 ) 


T . = t T 
out o 


+ - 1 - 


For L, „ « . 



For fixed values of t or Td* these mixer characterizations imply similar 
relationships between T MSSB and L M , and it is tempting to use these relationships 
in a tradeoff study. However, data taken with a mixer operating at 87 GHz and at 
115 GHz show that these simple models do not adequately describe mixer oper- 
ation. From Figure 1 the data taken at 87 GHz agrees with the mixer noise 
temperature ratio model, but data at 115 GHz bears no resemblance to the per- 
formance predicted by the simple model. 

Attempts to relate and TmsSB using a simple model of the millimeter wave 
mixer are unproductive due to the complex nature ol the noise and mixing pro- 
cesses. L M and T M SSB musfc both be considered i i evaluating mixer performance, 

References 

1. S. Weinreb and A. R. Kerr, ’’Cryogenic Cooling of Mixers for Millimeter and 
Centimeter Wavelengths, " IEEE J, Solid State Circuit s (special issue on 
microwave integrated circuits). Volume SC-8, February 1973, pp. 58-63. 

2. A. R. Kerr, "Low Noise Room Temperature and Cryogenic Mixers for 
80-120 GHz, " IEEE Transactions on MTT , Volume MTT-23, No. 10, 

October 1975, pp. 781-787. 

3. D.N. Held, "Analysis of Room Temperature Millimeter Wave Mixers Using 
GaAs Schottky Barrier Diodes," PhD thesis. Institute for Space Studies, 
Goddard Space Flight Center X- 130-77 -6, January 1977. 

4. Also see D. N. Held and A. Kerr, IEEE Transactions on MTT , Parts I and II, 
February 1978. 
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3a 2*3 IF CONTRIBUTION TO NOISE TEMPERATURE 

Contributions to the system noise temperature from the IF are determined bv the 
mixer conversion loss as well as the effective noise tempe rature of the IF network. 

The system noise temperature for a radiometer (DSB) was shown to be 


To = T A +I(T 


A ' 2 v A MSSB + L M T IF } 


It can be seen from the above relation that the contribution to the system noise 
from the IF is dependent on the conversion loss. The receiver noise temperature 


T 


R 


2 {T MSSB * L M T IF } 


and 


T = 
l S 


t a +t r 


Now the IF contribution to the receiver noise temperature is simply L^Ttf and 

vllueTa/] ™n£:V iValen \ t0 ^SSB- L M TlF is pl ° tted in Fi ^ e 1 *** 

iSShLiLJtW ?v as a function of conversion loss. It should be 

2ed . th f' t th ® contribution to the system noise temperature for a radiometer 
application ls just half of the plotted values. 

The IF noise temperature is generally dominated by the noise figure of the first 
stage of the IF amplifier. When the gain of the first stage is high and Ls followed 

Howev!r er - St M ge ° f gaiD J ith reasonably low noise figure, the above rule holds. 
However, in the proposed radiometer application, the IF may consist of a network 
of low gain amplifiers and filter banks where the overall effective noise tempera- 
ture is not entirely determined by the noise figure of the first stage. Figure 2 

hrofdh a e d ° ne P° ssibl c e , IF «<*work -here the first stage amplifier is a low gain 
then "d amplifier, followed by a demultiplexing filter bank where each baLl is 
then amplified separately. Let the noise figure of the first stage be F a and that 

mte h r lo^s oTl 1 r? X iS the individual demultiplexedtend. A 

is gTven by X lDtr ° dUCed for the Xth band - The noise figui-e for the network 


F AB = F A + 


Li x f bx~ i 

L G. J 
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IIB'JTION 
MPE RATI 


and the noise temperature is 


t ab = < f ab- 1)T o 


= t a + at e = ( f a- 1,t o + 




where ATjr is the excess noise temperature of the network over and above the 
noise temperature of the first stage. Values of excess noise temperature are 
plotted as a function of the gain-to-loss ratio of the first stage gain and filter loss 
for a second stage amplifier noise figure of 2. 


IF NETWORK 


mvAwm 

RMflK 


L x F BX “ 1 T 
AT - = O. 



10 11 12 


MIXER CONVERSION LOSS. L M . DC 

FIGURE 1. IF CONTRIBUTION TO RECEIVER 
NOISE TEMPERATURE 


FIGURE 2. EXCESS NOISE TEMPERATURE 
CONTRIBUTION DUE TO IF NETWORK 
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3. Design Analysis 




3o 3 RADIOMETER AT PERFORMANCE PARAMETRIC TRADEOFFS 


Radiometer AT performance is directly proportional to the system noise 
temperature. Typical values of Ts are estimated for the four radiometric bands 
and their dependence on feed and mixer conversion losses, and IF noise tem- 
perature are analyzed. 


The temperature i*e solution of the Dicke radiometer can be written 


AT 


where B is the receiver bandwidth and t is the radiometer integration time. Since 
the radiometer performance directly depends on Tg, it is instructive to examine 
the term closely. System noise temperatures for the four radiometric bands can 
be calculated from the expression developed in Topic 3.2. Certain assumptions 
must be made such as antenna noise temperature, mixer conversion loss, IF noise 
figure, and feed losses. Figure 1 lists tentative calculated values. 


The dependence of system noise temperatures and radiometric AT upon the included 
parameters can be analyzed by examining the partial derivatives of Ts with respect 
to the various parameters. From tlie general expression for syseiii noise 
temperature 


T s = T a + (L f -1) T 0 + T 


O ' "f R 


T R = 2 (X MSSB + L M t if ) 


we have 


T S = T A + ^ L F ~ ^ T 0 + 2 


Pmssb 




The relative magnitudes of the partial derivatives tells which parameters are most 
important in achieving and maintaining a low noise temperature. Talcing a parti- 
cular example of the 183 GHz band and the parameter values given in Figure 1, 
tlie partial derivatives and their numerical values are 
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3T< 

w 

all 


= T + • = 3818 

o 2 


bt, 


3F 


IF 


f- L M T o = 1152 


ax, 


9L 


M 


1T T IF = 3936 


9T, 


3F 


MSSB 


T =163° 
2 o 


In two of the expressions, noise figures were used instead of noise temperatures 
to make all of the expressions dimensionally equivalent. The results indicate that 
feed losses are of greatest importance, followed by IF noise figure, mixer con- 
version loss, and mixer noise figure. 

It should be emphasized that these expression are for small incremental effects. 
Large changes in parameters, such as 3 dB increments, give substantially dif- 
ferent results. Figure 1 illustrates how the AT values for the 1S3 GHz band change 
for a larg.* (3 dB) change in mixer conversion loss, IF noise figure, and feed 
losses . In the left column of Figure 1, the AT goals are listed for the six H^O 
channels; in thu right columns the calculated values of AT are listed for each of 
three values t ? mixer conversion loss. The value of L^t = S. 5 dB is selected as 
typical, in the next group of computations, the effects of increasing the IF noise 
figure are listed. In the final grouping of computations, the effects of increased 
feed losses are shown. 

Conclusions are that in striving to achieve ultimate performance, feed losses have 
tlie greatest incremental effect. However, gross changes in any of the parameters 
have almost a direct effect on performance. 
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AT FOR 

l m = 7 

L m = 8.5 

l m = ™ , 


“ RADIOMETER SENSITIVITY 
■ IF NOISE FIGURE, dB 

= IF NOISE TEMPERATURE, °K 
- FEED LOSS, dB 
= MIXER CONVERSION LOSS. dB 
= IF BANDWIDTH, GHz 
“ 250°K 
= 21G0°K 

“ 1.4 (REFERENCE AVERAGING) 
KTs 

%/STT 


; t A + <l f " 11 T o + -- (t mssb + l m t if ) 


FIGURE 1. RADIOMETER AT TRADES AT 183 GHz BAND 
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3. Design Analysis 

3.4 RE F E R E N C E A V E R AG1 N G 


Reference averaging can reduce the Dicke radiometer K factor Irom Z to some 
value of the order 1. 3 to 1.6 depending on chopping fraction and gain stability. A 
summary of the detailed analysis given in Appendix B appears below. __ 


The temperature resolution of a radiometer is determined by errors from two 
sources: l) random fluctuations inherent in the square law detection of noise 

power in a frequency band and Z) temporal variations ii: gain. I he use of refer- 
ence averaging” to improve a Dicke radiometer, first proposed by Bremer and 
Prince, 0) has the effect of reducing only tin* first type ot errors. the technique 
reduces the variance in the reference voltage by increasing its integration time 
for postdetection smoothing, which leads to a decrease in the overall variance ot 
the differential output voltage. In a switching radiometer with r "0 percent duty 
cycle instead of integrating over one -half the dwell time tT }? ) for a resolution 
element or “pixel,” we can integrate over Ntp/3 by averaging the relerence sam- 
ples in N consecutive pixels. Reference averaging can be implemented using a 
digital signal processor to store the required signal and relerence data, arid to 
perform the necessary averaging numerically. Such an extension ot integration 
time is possible for the reference samples because the reference temperature is 
essentially constant, whereas for the antenna samples the integration time must be 
limited to t /> to resolve the variation of scene temperature from pixel to pixel. 
For this case, the antenna temperature error 6T/^ n due solely to detector 
fluctions is given in the Bremer and Prince report. 1 -) 


6T 


S 


T 


An 


v/B 


K 


S 


/B 


Z I — 
N 


0 ) 


where Tg is the system noise temperature, B is the IF or predetection bandwidth, 
t is the pixel dwell time, N is the number of pixels over which the reference 
samples are averaged, and l\ is the radiometer sensitivite factor. In Equation 1, 
square wave modulation and detection are assumed. 

To determine how large N should be and to estimate the possible improvement, a 
detailed study was made to understand how reference averaging affects gain vari- 
ability errors (see Appendix B). We derived the following expression for the 
antenna temperature error 6 1 Ag due to gain variation in a relerence averaging 
radiometer for symmetrical chopping and averaging 



where G is the average of the time varying gain C'.(t) within a pixel dwell time; \\\, 
Tp. and arc the antenna, Dicke reference, and receiver noise temperatures, 
respectively. is an integral involving Git) defined as 





C N (t)G(t)dt 


( 3 ) 
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This case corresponds to minimum gain variation error 
lSY N t 1 D 4 1 K 1 


'* S SE"” - 

6T Ak »*'«"# drived from Kquation I and Figaro £ W " h * 1 An * nd 

further reduction in detector noise is possible by incorporating asymmetric 
• U T If W 1 reference averaging. If the fraction of time for antenna samnlina 
is eng umed irom 0. ->r to q r The reference integration time per pixel is 

nXuIo^" *t^l t0 bU< ,hiS iS ^^icted by averaginf otr 

pixels to gut a total reference integration time of N(l-q) n ,. Instead of 

Kquation 1. the error b 1 ,\ n is now a function of N and q 1 
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To do a similar analysis for this case, we need to know also M'Ac as a tvinction 
of N and q. This will be derived in a later study. Knowing 6TAgfN, 9 ) and 
Equation 4, we can then determine values of N and q that minimize the combined 
error 6 T a to find the maximum improvement achievable by reference averaging 
with asymmetric chopping. For any fixed q, we expect K (N) tor asymmetric 
chopping to behave in a manner similar to Figure 4. 

On the basis of Figures 3 rnd 4, we tentatively estimate very roughly that the 
K factor can be reduced by reference averaging to some value of the order ot 
1.3 to 1.6, depending on the chopping fraction used and the gain stability achiev- 
able. The upper value is based on symmetric chopping and the data in l* igure l. 
The lower value is based on asymmetric chopping and the assumption that gain 
variation errors can be kept small up to N Id. 
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The accuracy with which remote temperature can be radiometr ically measured is 
dependent on radiometer calibration accuracy, radiometer linearity, and size and 
rate of variations in radiometer gain, noise level, and bandwidth during the time 
interval between calibration and data collection. The calibration accuracy is 
examined below. 


The dc component of the radiometric output may be represented by 


V = a + bT + cT ' 

X XX 


+. 


where a and b represent the dc bias and gain constants of the radiometric receiver 
and T is the input temperature to the radiometer. Nonlinear terms cT x ~, etc., 
should be negligible. The radiometer will be carefully tested for linearity before 
launch and any required corrections will be determined. ( ^ 


Let T a represent the antenna temperature when the antenna is viewing the earth, 
Tl the antenna temperature when the antenna is viewing cold space, and T^ the 
temperature of the hot reference load. Given the output voltage for these three 
situations, the antenna temperature is 


( 


V 2 - V a\ 


r 


- v,/ l l - Vj/ L Z 


From this and the linear relationship between voltage output and input temperature 
the sensitivity of T a to errors in Ti and Tj may be written 


AT 


«. ■ 


T, - T 


AT 


Examination of this equation shows that the sensitivity to errors in T] is propor- 
tional to T^ - T a , and the term vanishes when T^ = T a . Errors due to uncertainty 
in measuring the radiometer output voltage because of analog-to-digital quanti- 
zation and nonlinearity will be smaller than the radiometer AT and can be 
neglected. 

Radiation at cold space and hot reference temperatures is reflected into the feed 
using a movable calibration mirror (Figures 1 and 2). The temperature of cold 
space is known and the efficiency of the cold look calibration will be measured 
before launch. 




FIGURE 2. ANTENNA SIDE VIEW SHOWING 
RADIOMETER CALIBRATION LOADS 






At frequencies above 100 Gil?., the radiometric temperature scale undergoes a 
constant depression in relation to the physical temperature scale, and at tempera- 
tures below about 10 °K is nonlinear (see Figure 3). This behavior is due to 
inaccuracies in the Rayleigh- Jeans approximation to Planck's radiation law 



whe re 

P = radiation power emitted by matched load 
k r Boltzman's constant 
h = Planck's constant 
T = temperature 
f =■ frequency 
B = bandwidth of interest 

Since we calibrate our system using a thermal emitter, we are calibrating to the 
physical temperature scale for temperatures above 10°K. However, below 10°K 
the concept of a linear temperature scale is no longer valid, and at 183 GHz, the 
zero radiation corresponds to a "temperature" of 4. 1 °K on our scale. Therefore, 
for calibration purposes we assign a temperature Tj to cold space of 4.2' 5 I\ when it 
is radiating at a Rayleigh-Jeans radiometric temperature of 0. 1 °K. 

The radiometer calibration procedure requires the use of a high emissivity micro- 
wave absorber at an accurately known temperature to provide a hot reference. To 
achieve the required performance, Brewster angle impedance matching techniques 
developed by Richard Iwasaki at the Jet Propulsion Laboratory are used. (2) 

The hot target is of a ridge design, using an iron filled, high dielectric constant 
epoxy absorber is cast directly onto an aluminum baseplate (Figure 4). Baseplate 
fins protruding into the absorbing material are used to decrease temperature 
gradients from the tip to the root of the ridges. Platinum resistance sensors 
accurate to ±0. 1 °K are used to monitor target temperature. 

The Brewster angle for the high dielectric constant materials used is relatively 
independent of frequency, allowing the absorber to be used for calibration over 
wide frequency ranges. 

The measured emissivity of the calibration target built for use in Nimbus -F by 
JPL was 99.6 percent. 

A Dicke radiometer using a chopper wheel reference load is used to reduce the 
effect of small variations in radiometer performance between calibrations. O) 
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To calculate an estimate of MASR radiometer temperature accuracy, the root sum 
of the squares of errors in T A due to estimated uncertanty in the calibration tem- 
peratures and due to the radiometer AT will be taken. 


The temperature of cold space is known, but uncertanty in the efficiency of the 
cold look aperture will introduce some error 


Tj = 4. 0 ±2 °K 


Radiometer AT = 0. 2°K 


AT 


T - T 

h iA 

T 2 - T, AT 1 




^ £a g£ 


+ (AT) 


blank ** 


r ? 7 “1 1 /2 

AT a = L (0 * 34) + (°- 2) J = 0. 39 °K 


The emirsivity of the hot absorber will be >99 percent, the temperature will b< 
monitored to ±0. T’K, and temperature gradients will be <0. 2°K. 


T z = -300 ±0. 5 °K 


AT A = 0. 47 °K 


The total uncertainty in T. includes contributions from errors in the two cal- 
ibration temperatures anda contribution from the radiometer AT 

(AT a ) 2 = (0. 39°K) 2 + (0.47°K) 2 + (0.2°K) 2 

AT a = 0.64 

Radiometer calibration accuracy better than 1 °K can be expected in all channels. 
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3. Design Analvsis 


3. 6 QUASI -OPTICAL DESIGN FORMALISM 
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Ray matrices are unimodular, that is, their determinants equal unity. Thus, a 
ray matrix M s of complex optical structures is built up by cascading matrices of 
individual elements M-* 

’a B ' 

M. = DET M. = 1 

1 CD. 1 

J l 


DET M. = i 


s l Z 3 i 


Ray matrices for four common optical elements arc given in Figure 2. 
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and their phase fronts are slightly cur . [ields emitted by lasers are 

tractable and physical. / reasona • tric Waveguides and HE^ modes in 

Gaussian; TEM waves in produced by scalar horns can be 

nation profile is illustrated in Figure 3. 
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the° distance z from the waist according to 


w 2 ( z ) = w 1 + 


tsfl 


The resulting beam contour w(z) 
fraction angle 0 is given by 


is illustrated in the Figure 3. The far field dif- 


e = \/tt w 


ius of curvature R of the wavefront changes according to 


The radius 

R(z) = z 


j\ + (ttWq^/Xz) 


As the beam propagates, R varies 


R = Ztt w q /X 


from infinity at the waist to a minimum of 


at 


z = tt w q / X 


symptotically. Also, the beam experiences a phase shift 0 


and approaches z asymptotically. — - --- 

relative to an ideal uniform plane wave given by 


- 1 2 « 
0 = tan (Xz/ttw q ) 
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3) Complex Beam Parameter . The real beam parameters w and R are 
used to define c complex beam parameter q by 


q 


R 



TTW 

O 


Using the q parameter the laws of beam propagation simplify to 


q 2 = q 2 + Z 

where q ^ and q^ refer to the input and output planes, respectively, and z is the 
distance between these planes. 

A thin lens of focal length f transforms the q parameters according to 


1 1 1 



where q, and q 2 are measured immediately before and after the lens. A beam 
passing through a general optical system characterized by its ABCD matrix are 
related by 


q 2 = ( Aq j 4- B)/(C q] + D) 

a solution which is accurate to within a phase factor 

6 < - j In M - —hJ: T \ . 

\ irW Q ) 

Because of the combination properties of ray matrices described earlier, a compli 
cated optical system of s elements is expressed simply as 


M = M. M,.M, M 

s 12 3 s 

A beam waveguide consisting of a series of lens elements represented by the 
individual matrix is simply 

M system = [ M l] ‘ 

where there a total of i lenses in the system. 
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4) Beam "Smith 11 Char ts. The complex beam parameter can be 
represented In a beam "Smith" chart in much the same way as a complex impedance 
parameter. Where the beam parameter is given by 1/q = 1/R - i2/b, the complex 
plane of i/q = i/R + 2/b is actually plotted. The imaginary axis is a plot of i/R 
while the real axis is a plot of 2/b, where b is defined as the confocal parameter, 
b = 2 ttw 0 ^/ X. . 

Any beam propagating in free space can be described by its radius of curvature at 
some point along the propagation axis and by its confocal parameter. Any propa- 
gating beam mode can be matched to any other through an appropriate set of inter- 
faces, such as lenses or reflectors. The process of mode matching assures 
unconditional stability of the system, which means that the only diffraction losses 
suffered are those introduced by the limited aperture, or truncation loss of the 
Gaussian beam. Figure 4 illustrates the use of the complex beam parameter 
chart to design a quasi-optical feed for a radiometer. In Figure 4, there are 
three modes described by the appropriate confocal parameter for each mode, 
bj, b 7 , and b^. 

The bj mode is determined by the characteristics of the antenna and is represented 
on the beam chart by a circle of radius 1/bj. It is instructive to "walk through" 
the system along the z-axis and to observe the techniques of mode matching. 

• Enter Chart at z = -a>. The point of beginning is from the antenna 
represented by zj = -co on the chart. Moving along the z-axis in 
the positive direction is the same as moving along the circle of 
constant b past zj = 0 to the lens L \ . 

• Lens Lj # The lens is the mode matching element to match the bj 
mode to the bi mode. We encounter the lens at approximately 

2 1 = +1 on the bj chart. The lens has the property of changing the 
radius of curvature of the wavefront. In this case, we want to 
change it from a jjositive value to a negative value suitable to match 
the bo mode. The length of line labeled Lj on the chart tells us the 
power of the lens needed. Its focal length is f = -1/Lj. 

• Jumping to Chart b>. In order to match modes, we must exit from 
chart bj at the same point on the imaginary axis as we enter chart b,. 
The beam chart representing region b^ Has a smaller radius than 
that of bj because the value of l/b^ is smaller. We continue to move 
in the positive z direction until reaching lens Lt. 

• Lens L 7 . The second lens must have sufficient power to match the 
strongly converging beam mode b3, thus the power of the lens is 
greater than that of the first lens. Its focal length is given by 
-1/L 2 . 

• Jumping to Chart Chart b3 represents the converging beam which 

must match the feedhorn to the mixer. The confocal parameter b^ is 
a small number because the beam waist is small and the radius of the 
circle in chart 3 is therefore large. We enter the chart from the 
matching lens at approximately the 7.3 = - 1 point and end at the feed- 
horn where the beam waist is smallest and the wave is a plane wave. 



3 . Design Analysis 

3.7 FABRY -PEROT FILTER DESIGN CONSIDERATIONS 

The use of quasi-opticaL Fabry -Perot tUter, ■ A^WUcL^techniques 

oscillator filters and for frequency multiplexing in MASK. Anaiyttca 
commonly used for optical systems are directly applicable. 


The technique of generating ofi^ 1 ^"^" ' s^ain.e 

plexor filter. 

The following design relationships are extracted by the general analysis given by 
Born and Wolf. ^ 2 > The relationships are given without proof for the purpo 
aiding'i'n the design application to MASK. The filter transmission characteristic 

is given by 


P ./ P. 
out in 


1 

1 + F sin'' ("7* ^ ^ 


and is shown in Figure 1 . The reflection 
= 1 “ p out^ P in» Nvhove the F parameter is 


characteristic R is given by R " P ro f p in 
defined by 


F = 


4R 

( 1 -R)‘ 


The sharpness of the filter vs enhanced for values o^R:y ? roachin, ^ 
th e value of F becomes very large. The term F sin l"/- ^ "/ 1 

f = f Q , but becomes very large for values of f * t Q . thereby dnving out' ln 

zero. 
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than 100,000. 
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FABRY PEROT FILTER RESPONSE 



GEOMETRIC 

LOSS 


WALK OFF LOSSES 



FIGURE 1. FABRY-PEROT PARALLEL PLATE FILTER CHARACTERISTICS 




Finesse factor is defined by the ratio of free spectral range to the 3 dB bandwidth 
of the filter transmission band. The characteristic of the filter is such that the 
rth dB bandwidth is related to the 3 dB bandwidth by 

B r = B 3 (lO Wl ° - l) 

Table 1 lists these characteristics and other important design relationships of the 
Fabry-Perot filter for applications to MASR. 

Practical application of Fabry-Perot filters to the millimeter wave region has been 
achieved by Goldsmith ^ and Gustincic^. 
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4. RECEIVER DESIGN 


Where possible, recommended receiver component designs are described. In 
"soft" technology areas, such as mixers and local oscillator sources, various 
possible approaches are examined. 


Mixer Design Alternatives . Three mixer designs look attractive for MASR. 
Selection of a preferred configuration from among these alternatives depends on 
development progress in waveguide balanced mixers, quasi-optical mixers, and 
subharmonically pumped mixers. 

Mixer Diodes. High quality vapor phase epitaxial gallium arsenide material is 
difficult to obtain, and a batch selection and impound methodology is recommended 
for MASR flight programs. A conventional diode array of 1 to 2 pm dots, whisker 
contacted in a Sharpless wafer, is the recommended diode design. 

Local Oscillator Sources . IMPATTS are primary contenders for use as LO 
sources in MASR.. Quasi-optical filters used with balance mixers will provide 
sufficient LO noise rejection. The Gunn oscillator is unsuitable for use as a 
fundamental local oscillator due to a present high frequency limit of 100 GHz. 
Superior Gunn noise performance makes it an attractive alternative to the IMPATT 
when used with a doubler or with a subharmonic mixer. 

Quasi-Optical Local Oscillator Filters . Fabrey-Perot quasi-optical filters have 
a much higher unloaded Q than waveguide filters, resulting in improved selectivity 
and lower loss. 

IF Design . IF networks consisting of filters, diplexers, and amplifiers are 
described for the 118 and 183 GHz channels. 
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4. Receiver Design 

4. 1 MIXER DESIGN ALTERNATIVES 


Three mixer designs look attractive for MASR. Selection of a preferred 
configuration from among these alternatives depends on development progress 
in waveguide balanced mixers, quasi-optical mixers, and subharmonicallv 
pumped mixers. 


During the course of the present study, a number of mixer configurations were 
analyz'd and evaluated; they include single-ended, waveguide hybrid balanced, 
crossbar balanced, quasi-optical balanced, and subharmonically pumped balanced 
configurations (Table l). The characteristics of each approach are summarized 
below, 

lhe designs treated in detail in the succeeding topics of this subsection were 
selected as the most appropriate for the MASR application Of those configurations 
that provide both reliability and high performance, the waveguide balanced mixer 
probably involves, at this time, the least technical risk. However, further devel- 
opment work in the subharmonic mixer may improve confidence and swing opinion 
in its favor. 

Singte- Ended . Single-ended mixers have been used more than any other type at 
these frequencies in ground based radiometers for radiotelescopes. The lowest 
reported conversion loss above 100 GHz has been achieved with a single-ended 
mixer. The disadvantages of this unit are difficulty in LO injection and extreme 
sensitivity to LO noise. It is estimated that 45 dB of LO filtering is required with 
IMPATT local oscillators. Filters having sufficient noise rejection at vhe IF fre- 
quencies also have high ioss at the local oscillator frequency. Thus, the use of a 
single-ended mixer in the MASR system depends upon the availability of sufficient 
local oscillator power of the required spectral purity. 

Wavegu ide Balanced Mixer . The classic waveguide balanced mixer using a magic 
T or a short-slot hybrid has been successfully used at 118 GHz and can be scaled 
up to 183 GHz. A balanced mixer has the advantages of local oscillator noise sup- 
pression, high local oscillator to RF isolation, and good IF source impedance. 

The main disadvantage is the added complexity due to the two mixer diodes per 
mixer. The use of a waveguide balanced mixer in MASR will be constrained to a 
top-wall hybrid or lolded side- wall hybrid in order to provide contiguous LF out- 
puts of the two diodes, minimizing the capacitance to help achieve a uniform 
impedance over the very wide IF bandwidths. 

Crossbar Balanced Mixer . The crossbar balanced mixer has a good IF impedance 
characteristic that makes it attractive for use with very large IF bandwidths. The 
mam disadvantages of this design are poor local oscillator to RF isolation, poor 
LO noise suppression, and poor reliability. The crossbar design has been used 
exclusively with beam-lead or pill-packaged diodes, limiting its use to the lower 
millimeter wave frequencies. Mechanical and thermal ruggedness of the crossbar 
mixer with whisker contact dioc’es remains to be proven. 

Subharmonically Pumped Mixer. Subharmonic mixers have the distinct advantage 
of using a local oscillator at half the signal frequency. The excellent performance 
reported has been achieved at the lower millimeter wave frequencies using beam- 
lead diodes. Further development is required to prove that the subharmonic 
mixers can perform at the higher signal frequencies. The subharmonic mixers 
have the advantage of local oscillator noise cancellation. 
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Quasi-optical Balanced Mixer . The quasioptical balanced mixer is designed to 

use two single- ended waveguide mixers in a quasioptical hybrid assembly. This 

technique should be considered as a method of achieving local oscillator noise 

suppression without ':he usual concomitant losses. 
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TABLE 1. MIXER DESIGN ALTERNATIVES 


Balanced Waveguide Mixers 

o Topwall hybrid meets conditions of proven design and special 
low output capacitance needed for MASR IF bandwidth 
requirements. 

Quasi-optical balanced mixer 

• Quasioptical balanced mixer permits balanced effect to 
be achieved with two single-ended mixers. 

Subharmonic mixers 

• Subharmonic mixers offer good mixer performance with 
half local oscillator frequency. Local oscillator noise 
cancellation in IF is inherent in system. 
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4. Receiver Design 

4. 1 Mixer Design Alternatives 

4. 1. 1 WAVEGUIDE BALANCED MIXER 

Topwall or folded sidewall hybrids meet the requirement of Local oscillator noise 
cancellation and the special requirement of low output capacitance for wide MASR 
IF bandwidths. 


Balanced mixers have long been used to provide immunity to noise in the local 
oscillator and to provide convenient local oscillator injection. (*) Special wideband 
IF requirements of MASR are a special constraint for the output of the mixer. In 
the balanced mixer, output terminals should be joined and matched in the optimum 
characteristic impedance to the IF circuit. Physical distances between the mixer 
whiskers and the point where they are joined should be minimized. The waveguide 
hybrids that can meet this requirement are the topwall hybrid and the folded side- 
wall hybrid. In both ox these structures, the IF outputs of the mixers are contig- 
uous and can be conveniently matched to a transmission line to the IF circuit. 

Figure 1 illustrates the topwall hybrid assembly as it may be used for the 183 GHz 
mixer. The signal enters through a scalar horn in one arm of the hybrid while the 
local oscillator enters another horn. After passing through the hybrid, the outputs 
are matched to an identical pair of mixer diodes mounted in Sharpless wafers2,3 > 
The wafers are shown separated from the interface flange, their respective out- 
puts in close proximity. The IF outputs are individually filtered in the Sharpless 
mounts through an inverted LP coax filter. When mounted to the interface flange, 
the IF outputs of the individual wafers are joined through a coax feedthrough that 
has a characteristic impedance of about 100 ohms. The 100 ohm output is matched 
to the IF network, which may be a wideband amplifier (buffer) or a multiplexer 
filter network. 

The advantages of the design are 1) a standard design that has been used at the 
lower microwave frequencies for many years; 2) identical diode waves — a fea- 
ture that facilitates the testing, evaluating, and selecting of matched pairs; 

3) optimum MIC geometry from mixers to IF; 4) demountable mixer diodes tor 
easy replacement. 
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4. Receiver Design 

4. 1 Mixer Design Alternatives 

4. 1.2 QU ASI- OPTICAL BALANCED MIXER 


The quasi-optical balanced mixer permits local oscillator noise cancellation to be 
achieved with two single-ended mixers and a quasi-optical hybrid assembly. 


The quasi-optical balanced mixer using two single-ended mixers is made possible 
by combining the received signal and local oscillator through a polarization 
diplexer. Both signal and local oscillator (combined) signals are then divided 
with a second polarization splitter and focussed onto separate mixer mounts in 
such a way that IF variations caused by local oscillator noise are out of ph«. se at 
the mixer outputs and are thereby cancelled with these outputs combined. 

For radiometric applications, the incoming signal polarization is random. Thus, 
selecting either one of the two incoming polarizations results in no difference in 
the signal quality. As shown in Figure 1, the polarization perpendicular to the 
paper is selected and reflected toward the top of the figure while the rejected 
polarizaiion is transmitted by the wire grid polarizer 1. At the plane of rota- 
tion, the signal polarization is rotated 45°. This results in two components, one 
into the plane of the paper and one to the right in the plane of the paper. The local 
oscillator signal, which is incident on polarizer No. 1 from the bottom of the fig- 
ure, is orthogonal in polarization to the signal and is transmitted through polarizer 
with no (or little) attenuation. After the rotation plane, the local oscillator polar- 
ization is shown separated into two components, one out of the plane of the paper 
and one to the right in the plane of the paper. Components of the signal and local 
oscillator in the plane of the paper are transmitted through the second polarizer 
(No. 2), reflected off the mirror at the top to the upper mixer mount. Components 
of the signal and local oscillator perpendicular to the paper are reflected off polar- 
izer 2 to the right, reflected off the mirror at the right and focussed to the other 
mixer mount. 

The signal and local oscillator are in phase in the upper mixer and are out of 
phase at the lower mixer. This results in the production of an IF signal at the 
upper whose noise components are out of phase with the IF signal at the lower 
mixer. Noise cancellation effects can be illustrated with a simplified analysis. 
Consider the simplified schematic diagram of the mixer in Figure 1. The signal 
and local oscillator fields can be represented at point 1 by 


1 ) 


E (t) = E cos (u t + <j> ) 

s s s Y s 

P = 

-90° - 

-45 

E LO (t) " E LO COS ‘"LO 1 + <t> LO ) 

P - 

0° - 

45 


o 


o 


The polarization of these fields with respect to the plane of the paper looking along 
the direction of propagation is indicated at the right. The changes caused by the 
polarization rotator are also indicated. 


"' This technique uses waveguide mounts for the mixer diodes and therefore differs 
in concept from the quasioptical single-ended mixer reported in J. J. Gustincic, 
"A Quasi-Opticai Receiver Design, " Symposium on Microwave Theory and 
Techniques, June 1977. 
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FIGURE 1. QUASI-OPTICAL BALANCED MIXER 
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After passing through the second polarizer, both the signal and local oscillator 
..elds are split evenly and directed to mixei s a and b. The fields prior to enter- 
ing the mixers at points 2 a and 2 b are 
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Upon entering the mixers, these fields produce voltage outputs which are proper 
tional to the square of the vector sum of the signal and local oscillator fields 
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which combined yield 


4) V a > V b + 2KE.E L 0 co.,- if t t ^ 


*o *s • * 1.0 


lh, condition I placed on .this analysis is that t( 5 s - ( 3 j 0 l ( a - b) « */ 2 , and since p* 
and PLO determined by the local oscillator frequency and (highest) signal fre- 
quenev , the maximum allowed difference in path length, a and b, can be estimate* 


(a - b) <*; 


' 4 ll s- f l. C ' 


« 7. S mm 


1 his condition is easily achieved with state of the art mechanical 
techniques. 
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Subharmcnlc mixers using loo ill oscillator pump that us at half the normal l O 

frequency with an antiparallel divide pair ran provide mixer conversion loss that 
closely approximates that of a normal mixer. In addition, this configuration pro- 
vides local oscillator noise suppression* 1 1 • •-) 

I he subharmonic mixer using halt the normal I O frequency, 1/2 x t'j 0 , and an 
antiparallel diode pair can provide very efficient conversion of a signal from the 
frequency fs to 1 1‘ | ~ t's I • lids efficient mixing is accomplished because the l .O 

power directly produces a time-varying conductance with a fundamental frequency 
variation at tj,o. This face is demonstrated in Figure l, which shows the conduc- 
tance waveforms for a single diode and a diode pair pumped at 1/2 x fj q. The diode 
pair has a time variation that is twice the pumping frequency. As in a standard 
mixer, the trequenev conversion process in the subharmonic mixer occurs through 
the interaction of the signal power with the time-dependent conductance of the 
diodes. Since the primary troquency in the conductance waveform is fj q, the prin- 
cipal mixing products are generated at t*LO *l*S* the same as in a standard mixer. 

The signals generated at fjj.* 1 Ia'“ *S l J V the interaction ot the 1.0 noise with 
each diode have phases which cause them to cancel in the IK circuit. Imbalances 
in the diodes which create unequal conversion loss for the two diodes will degrade 
the noise cancellation. Likewise, these imbalances will cause the time-dependent 
conductance to have some half fi^ contribution* This will slightly degrade mixer 
performance, hut should not be too serious in applications where sufficient 1.0 
power is available and where t\K V v t's* 

Hie question of how much 1.0 power is needed has not vet been determined ana- 
lytically. However, experiments by Hell Labs at millimeter waves have shown 
that diodes with low cutott frequency or large parasitic elements require much 
higher pump power levels t~l 0 dm. O.ood mixer performance has been obtained 
with less than 10 dHm, while poorer diodes in the same circuit need 10 d H more 
l .O power. 

Ihe best results have been obtained with divides that have cutoff frequencies, tj c >, 
-20x the signal frequency of the mixer. Kor the ISO CHr. application, this implies 
that diodes with tj about 4000 CHr, are needed, Ibis figure of merit is probably 
beyond the present state of the art and can probably be met by utilizing the ultimate 
capabilities ot Ga As materials technology and millimeter wave diode fabrication 
techniques. 

Given the best diodes, it is still necessary to devise a mixer circuit with parasitica 
low enough not to degrade the device capabilities. The circuits that have performed 
best at millimeter waves have generally used some modification of the Sharpless 
water tv' contain the diode. I he circuits that we propose for the subharmonic 
mixer also use that approach for embedding the diode in the RF circuit. Figure .» 
s how s two proposed mounting schemes for the divides. Common to both approaches 
are the distributed tow pass titters which complete the RF circuit for the diodes 
and allow IF power and bias currents to pass from the RF circuit region. With 
proper diode polarities, the close physical proximity ot the diodes allows the IF 
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4. 2 Mixer Diodes 

4.2. 1 EPI- GAL LIU M ARSENIDE MATERIAL 

Mixer performance is critically dependent on high quality ;vapor phase) epitaxial 
gallium arsenide material. Sources of this high quality material are extremely 
limited. Batch selection-and- impound methodology is recommended for MASR 
flight programs. 


Performance of a mixer diode ir dependent on geometry and upon the material. 
Candidate materials include silicoi , gallium arsenide, and indium phosphide, in 
either a liquid phase or vapor phase epitaxial form. It is generally believed that 
the cutoff frequencies achievable with silicon are far below those achieved either 
with gallium arsenide or indium phosphide. Given the results reported to date, 
gallium arsenide appears to be the preferred material for frequencies from 100 
to 200 GHz. Although indium phosphide shows promise of having even higher cut- 
off frequencies, its state of development is far behind that of gallium arsenide, 
and its selection as a candidate material involves greater risk. 

The two types of gallium arsenide that have been considered are liquid phase and 
vapor phase epitaxial (VPE) material. Mixer conversion loss and noise figure 
are a critical function of the junction forward resistance, and whisker contacting 
of vapor phase material produces lower resistance and lower capacitance junc- 
tions. Some of the producers of vapor phase material are: 

• Hughes Electron Dynamics Division (Torrance) 

• Raytheon Research Center (Sudbury) 

• GHz Devices (Cambridge) 

• British Drug House (London) 

• Plessey, Ltd. (London) 

Production of quality VPE material suitable for low noise mixers is at present 
achieved accidentally, since full understanding of the growth control process has 
not yet been achieved. However, most of the above laboratories are currently 
improving facilities to accommodate a strong interest in Gunn diodes, FET ampli- 
fiers, and microwave mixers at the lower frequencies. It is reasonable to expect 
that these facilities will be able to produce batches for evaluation for mixers in the 
200 GHz region. Certainly, one identified batch of high quality material is s< *i- 
cient to provide all the mixers needed for a MASR flight program. 

Evaluation of tin material is not a simple process; it involves fabrication of diodes, 
contacting, mounting, and testing of samples from each batch produced. The 
methodology of EPI gallium arsenide evaluation and selection is illustrated in 
Figure 1. The recommended approach is to procure batches from a number of 
sources, fabricate sample diodes, and evaluate each of the batches from each 
source. When the quality needed for MASR has been achieved, the high quality 
batch is impounded for the exclusive fabrication of mixers for the program. 
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FIGURE 1. EPI GALLIUM ARSENIDE INVENTORY METHODOLOGY 
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4. Receiver Design 
4.2 Mixer Diodes 


4.2.2 DIODE DESIGN 


^ea^Tith^ography,” is Recommended 1 a ^Sha ' |° 2 d ° ts ’ fabricated with electron 
diodes are the most suitable packages for the^O ?o e 200 C GH^“gion hiSker C ° ntacted 


co T?r ? methods ° f fa -- 

achieve with photolithon raohv ^^7 J smallest dots practical to 

have been made by this method Electron b^am a rVv, ter ’ ^ S ° me excellent diodes 
sub -pm dots of great precision as well T lltho e ra P>W is useful in producing 

stripe geometry These two ronf . ° tlC 8 eometr ies such as the cross 
stripe device has the Xta,. ^of^" ^ V ShOWn in Fi 8 ul ' e 1- The cross 

and thus produces a higher cutoff freoTm vTlIa S *‘ eslstance and shunt capacitance 
However, the simpler dot geometry of -2 L d? the f conventi ° nal dot geometry. (D 
MASR applications. 8 Y of 2 ^ diameter is probably adequate for 


Scroscopes' or’sc^nlng eUcfron m.cros cTp^'^SEMst" Thews’ 

fmiKhic"ed in ^ hi s ^hion areVery' low* 8 SEMs'g . COntactl “?« aad yields'of diodes^ 1 
accurate contacting hut have the dis advantn > ^J ve SU P e . ri01 resolution and permit 
,nd the use of micfomani^a^^ environment 


in the diode RF equivalent circuit/ For lo ^ Y Velated to the P ar asitic elements 

f co> has to be at least ten times the ODeratin/f VG1 * Sl ° n l ° SS * tlie cutoff frequency, 

mixers), and the series resistance should l 8 , eqUenCy (greater for subharmonic 
quency is given by ^ce should be as low as possible. The cutoff fre- 


f co = (2 ttR C)" 1 


whe re 


R 


= forward series resistance (8 ohms for best mixers) and 


C = junction capacitance (0. 007 pF for best mixers). 


(i) 


/• 

The calculated cutoff frequency is thus 2 8 x 10*2 - T , u . , , . . 

suitable mixer at frequencies up to 300 GHz. ’ hlch should produce a 


1F output terminal. Figure 1 

can be provided to thTXsker tSrou.h the* ° , a ° d ° Ut P ut IF Pi". DC bias 

local oscillator filtering in the IF ° s provided bl^ anod ‘ Zed whisker P°st, while 
built into the IF pin. ^ t ^ le inside-out low pass filter 


o 


0 


Reference 


^glon.^S Lj’^artton". ' *> \ »» Wavelength 

vox. No. 12, December 1974. 

A r-> * . 


2 . 


^o^r-yki^S-s.^: 140 *• 
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4. Receiver Design 

4.3 Local Oscillator Sources 

4.3. 1 TRANSFERRED ELECTRON OSCILLATORS 

The Gunn, or transferred electron oscillator is of about 

fundamental local oscillator because of a present hlg attractive alternative 

inn GHz Sut.erior Gunn noise performance makes it an attractive 

to the IMPATT when used with a doubler or with a subharmonic mi- ■ 


Materials. Interest in indium phosphu^ ^ a ^uM matein^l ^s grmsdng^because 

tion tolerances, and improved noise performance. 

fg“r 

may extend this limit to ~150 GHz. 

Reliability . The reliability of Gunn devices U^tego^ over 

by Varian has established that the MTBF at tl e P h (1) Three 55 GHz 

m* T - 

“d P .S miiri.‘“r.b\« 2!er„,U». .. tb. filtered IMPATT-balanced mla.r at 
183 GHz. 

; the meratare. (=•« 1 4 ' 1 
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4 * 3 Local Oscillator Sources 

4.3. 2 IMP A IT SOPRCFS 

Power , frequency ws * d withbal " 

'."fflclent LOnoUe -ejection. 




~ ... *• T\lP ATT oscillators are 

Power output and frequency of . aU channels using the 1°' ' 

Tiiown" in Figure 1. Adequate P^^rto-.e-. A minimum available LO power 
•tat. of the art and considering LO f ^ lo(ie 


Powe r 


shown in Mguu- 1 . 1 , nfi . tt , r i osse s. A minimum ' U,UU1 " ‘ , . 

state of the art and considering LO td ^ - s tho narrow LO Liter that »• 

of o niW occurs at Mb CHa CH , LO power will become available 

^ iuToTtU 

Reliability. IMPATT lifetime is 

- •*« - - Arrhenius relationship 


T ex 
o 


x p(ff) 


where r = 1.6 x 10 -“ hours the e^ti on ship^r edicts a median 



Device defects are responsible for causing ^hip^EUmtaVtion'of these defective 
would be expected from a „d an initial burnt,, will result 

iXSS'Z&Vor the MASR program. 

No,..', a „ pio., w o 

high IF frequencies. The meast u -hich implies a noise spectrum 

(Figure 2) is reasonably flat le also varies with circuit 

flat ‘beyond 12 GHz tor a ISo GHz INU Mi. 
parameters and bias current. 

LO noise rejection of thC 

required from the quasioptical LO W th e 20 MHz IF band edge present. 

Chosen over cavity filters because of 

their much lower losses. 

SI ability . The frequency drift ^o“l 2^^ r ^the 

K^rsrr«sws. «. -» - — »—* 

Day-to-day IMPATT -»« <• “ «• ” *** **>"« ** *" 

to aging i ^ not available. 
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DISTANCE FROM CARRIER. GHZ 

FIGURE 2. AM NOISE AWAY FROM CARRIER 
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Oscillat or Circuit . Figure 3 is a cross section of a Y-band (1/0 to 260 GHz) 

IMP ATT oscillator developed at Hughes.™ It consists of three major sections: a 
tapered waveguide section that transforms full height waveguide to reduced height, 
a reduced height waveguide wafer section that contains the IMP ATT diode, and a 
mechanical tuning short section. The diode mount resembles the Sharpless water 
mount used in mixer diode applications. I’se of the wafer mount facilitates the 
variation of circuit characteristics at the chip level, allowing the development of 
a diode package with optimum parasitics and good thermal properties. A detaile 
description of millimeter-wave IMPATT design and development procedures is 
uiven in the references listed below. 


Phase Lock System. The IMPATT local oscillators will be phase locked to a 
-Hible reference tb assure long-term frequency stability. < 5 > A block diagram of 
the proposed phase lock system is shown in Figure 4. A sample oi the 1\11 A I I 
output is mixed with a harmonic of a lower frequency stabilized Gunn source, 
generating an IF signal close to the crystal reference. This IF signal is then 
locked to the reference by varying the IMPATT bias current, resulting in a phase 
locked IMPATT frequency of t's = Nf LO ± f R* 


Locking to the proper sideband, above or below the nth harmonic of the Gunn is 
assured by the use of two phase detectors working in quadrature. The output o 
the quadrature phase detector is sideband dependent, and is used with the search 
oscillator to move the IMPATT toward the desired sideband. 
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4. Receiver Design 

4.4 Quasi-optical Local Oscillator Filters 

4.4.2 PARALLEL PLATE VERSUS RING RESONATOR CONFIGURATIONS 

Parallel plate Fabry- Perot filters have the advantage of arbitrarily large free 
spectral range while the ring configuration offers a convenient way of combining 
signal and local oscillator. The large IF bandwidth requirements for MASR make 
the use of parallel plate configurations mandatory. 


An important requirement for a local oscillator filter is that the filter passbands 
have sufficient spectral separation to avoid noise bands in the IF of the receiver, 
for the individual MASR bands, the required spectral separations arc: 

183 GHz band, maximum IF frequency 10 GHz 

14 0 GHz band, maximum IF frequency 2 GHz 

118 GHz band, maximum IF frequency 4. ! GHz 

104 GHz band, maximum IF frequency 2 GHz 

While one of the LO filter passbands must lie at the center frequency of each line, 
the first passbands on either side must be* separated from the center frequency by 
at least the values listed above. In a Fabry- Perot filter, the separation of pass- 
bands is called the free spectral range, FSR, and is given by FSR = c/2d, where 
c is the velocity of light and d is the separation between the plates. ^ This relation 
yields the following values of maximum plate separation for the individual bands: 


183 GHz band 


d < 15 mm 


140 GHz band 


d < 7 5 mm 


118 GI Iz band 


d < 37 mm 


104 GHz band 


d < 7 5 mm 


l or a parallel plate tabry- Perot configuration, plate separations of values as low 
as 2 to 3 mm are easily achieved, while a ring resonator must maintain plate 
separation determined by the aperture of the device. 

Figure 1 shows a parallel plate l abry- Perot filter with its mirror surfaces 
slightly concave to correct for diffraction. The effective plate separation, d, is 
very nearly the physical separation at the mirror centers. The application for 
MASR is for the local oscillator source to be connected to the horn labeled P in 
while the P out horn is mounted to the mixer. Since out of band energy is reflected 
back to the source, a slight tilt (typically 3°) is required to avoid affecting the 
performance of the LO source. An isolator may alternately be used where tilt 
losses are too severe. This type of configuration is applicable to anv of the MASR 
bands. 

1 he ring resonator with local oscillator injection and provision for simultaneous 
LO filtering are also shown in Figure 2 (see subsection 3.8). It can be seen that 
the resonator spacing, d, is related to the berm waist, w G , and that the free 
spectral range is determined by the single pass distance, d, through the resonator 
Ihe speciiic example for a 140 GHz configuration is shown where, for the MASR 
tentative design, w Q = 27 mm. The plate separation is thus greater than 2 16 mm, 
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and the free spectral range FSR is less than 1.38 GHz. Since the IF frequencies 
for this band may extend to 2 GHz, the ring resonator is not suitable for the 
MASR application. 

The unique properties of the ring resonator for local oscillator injection should 
not be overlooked. As the signal is incident on the filter, energy at or near the 
140 GHz frequency enters the filter and is lost while energy on either side of 
140 CHz is reflected intc the 140 GHz mixer. From the local oscillator side, 
energy at precisely 140 GHz is transmitted through the filter to ch? mixer while 
local oscillator noise on either side of the 140 GHz frequency is .reflected at 90°. 
This convenient property of ring resonators also prevents reflected energy from 
returning to the IMPATT, which may otherwise cause instabilities. 
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FIGURE 1. PARALLEL PLATE RESONATOR (GOLDSMITH) 













I' 4 ^to^Local Oscillator Elite r. 

4 ' 4 3 n.TEK DESIGN FOR .83 GHz LOCAL OSCILLATOR 

. , , total insertion loss of 1 dB and a 

TI... 1S3 GHz local oscillator Hit c> J ^ ! ' (Tslll , { , ,) and Figure 1. 

iiHhand rejection ot morcto ^MIU 

TABLE 1. 183 GHz LOCAL OSCILLATOR FILTER 

CHARACTERISTICS 


Type 

3 dB bandwidth 
10 dB bandwidth 
Fiee spectral range 
Beam waist 
Aperture diameter 
Walkoff diffraction loss 
Reflectivity 
Plate spacing 
Ohmic loss 
Midband rejection 
I Walkoff geometric loss 


Fabry-Perot parallel plate 

381 MHz 
1200 MHz 
10 GHz 
22.4 mm 
70 mm 

0.14 dB (without correction) 

0.91 
1 3 mm 
0.42 dB 
26.5 dB 

0.4 dB for 3° tilt (no isolator used) 


issa 


LOSS O 4 ilB 


183.3 

FREQUENCY, GHZ 


P.GURE, 183 GHZ LOCAL OSCILLATOR FILTER RESPONSE 




4. Receiver Design 

4.4 Qua s i- opt i cal Vocal Oscillator Kilters 


.1.4.4 OKSIO.NS FOR US GHz l OCA l OSOU.l.A TOR Vll U**K 


Tho most difficult filter for MASK application is tor the l IS C.ll/ l-u al om II u 
where IV frcqvuMic v of CO MU.* must be filtered. I he largest conti lbution to ^ > > 
ilui< front ohmic loss *8. 1 dlD, since diffraction loss can be correcico toi am 
UuZes Z be " voided by u.in B an isolator. . ho filter oharactor.st.es and 
response are suntnta rir.od in lahle l ami 1‘ V^ ul c 


TABLE 1 1 1S GH; LOCAL OSCILLATOR Fit TER CHARAC1 ERISTICS 


Typo 

3 dB bandwidth 
10 dB bandwidth 
Fieo spoon al i anyo 
Brain waist 
Apoi tuio diamotoi 
Walkett ditti action loss 
Ret activity 
Plato spacing 
Ohmic less 
Midband rejection 
Walkett geometric less 


1 abiy Porot paiallol plate (and collected plate) 
26 MU.- 


without 

collection 

optics 


32 A mm 
70 mm 
0 dB 


with 

collection 

optics 


•JO MH: 

•1 GH/ 

48.3 mm 
160 mm 

1.3 dB 
0.99 
38 mm 

8. 1 dB 
54 dB 

8.2 dB toi 3° tilt (isolatoi lecommended lathei 
than tilt) 

4.9 dB toi 1.8° bit 
0 dB toi no tilt 



FIGURE 1. 1 18 GH*: LOCAL OSCILLATOR FILTER RESPONSE 






4. Receiver Design 

4.4 Quasi-optical Local Oscillator Filters 

4.4.5 DESIGN FOR 104/140 GHz LOCAL OSCILLATOR FILTERS 

Filters for the 104 and 140 GHz window bands are of similar design. A midband 
insertion loss of 1.4 and 1.9 dB with walkoff losses of 4. 6 dB is achievable with 
a midband rejection of 26 dB. Filter characteristics are summarized in Table 1. 


TABLE 1. 104/140 GHz LOCAL OSCILLATOR FILTER 


Type 

Fabry-Perot parallel plate 

3 dB bandwidth 

64 MHz 

10 dB bandwidth 

200 MHz 

Free spectral range 

2 GHz 

Beam waist 

22.4 mm 

Aperture diameter 

70 mm 

Walkoff diffraction loss 

1.1 dB (without correction) 

Reflectivity 

0.89 

Plate spacing 

75 mm 

Ohmic loss 

1.4 dB at 104 GHz, 1.9 dB at 140 GHz 

Midband rejection 

25.6 dB 

Walkoff geometric losses 

3.5 dB for 3° tilt 
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4. Receiver Design 
4. 5 IF Design 

4.5. 1 NETWORK FOR 0.625 TO 10 GHz IF 

A network of contiguous diplexers, filters, and amplifiers can provide satisfactory 
li multiplexing and noise figures for each of the six H.,0 channels. 


reLhlT,Z l 'n e ;' rt ,° £ C0 ." tigU0 “ S filtcr banks is progressing rapidly and has 
leached the point where it can be stated with assurance that a six-wav evanescent 

^ch e a C Wte e r U v OI,S |Th ltlpl r er , t0 meCt the 0, 625 to 10 GHz requirement is feasible* 
imoedanee t It” ° W °, SS '' vould present a substantially constant input 

for U f | r ' 1 ’ miXer OVer 4 he entirc band - However, the design procedures 
loi such fillets aie new, and a custom design requires further effort beyond the 
hcope ol this report. 

frTd!rll°Hv l de * us ing contiguous diplexers in cascade are more common and 
examnb^nrf applicable to the problem. Such a design has been worked out as an 

critlcal P a rt of the filter is the impedance match between the mixer and the 
first diplexer. From the output of the hybrid balanced mixer, the filter input 
coax connects to the junction of a high pass low pass diplexer. The nominal 

aXd mix e er T, ^ thc a PPr°*i™te output impedance of a bal- 

cro! d , T ln )P edance of each tiller is also 100 ohms, increasing at the 

of the°two 5 At the crossove r frequency, the shunt impedance 

oi the two parallel branches is also 100 ohms. 

JonH 5 D,!n„ 1 s°d GI l iZ P ° rti , 0n is again s P Ut into two bands in another evanescent mode 
contiguous diplexer whose outputs cross over at 7. 5 GHz; each band feeds an IF 

amplifier. Similarly, the band from 0.6 to 5.0 GHz is also split into two bands 
which cross over at 2. 5 GHz. Finally, the lower bands are each split agafn so 
tha six bands are produced, each with a 5 0 ohm output and each with an IF ampli- 
fiei. Figure 2 shows the complete 183 GHz band IF circuit with six output da- 
tive IF note y fig g ures frequencies covered by each channel and the respec- 

Using the relationships developed in Section 3, and the effective IF noise figures 

fit ar f c f lcu l ated which include the effects of filter losses. The radiometer iT-'s 

coLTrvatile 8 f ° r a convent > ona l Dicke radiometer. These calculations assume a 
consei v ative 8.5 dB mixer conversion loss and an 0.5 dB feed loss A tV nical 

antenna temperature of 250 I< is also assumed. It should be noted that the cal-u- 
latod values ot AT exceed the goals by about 1.7 times. These values may be ~ 

X by achievin 8 a lower mixer conversion loss or by the use of reference 
averaging in the processing of data. 


"'For example, the Filtronics (S 
12 GHz 1.5 dB maximum loss. 


anta Clara, California ) quintuplexer , 0. 6 to 
55 isolation, 2. 5 to 1 maximum VS WR. 
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5 GHz 

LOW PASS FILTER 
0 6 TO 5 0 GHz 


FIGURE 1. WAVEGUIDE BALANCED MIXER OUTPUT INTERFACE 
TO IF NETWORK 



BANDPASS 

EVANESCENT 


71 




TO 10 GHz 
- 5.5 dB 


HI 


BANDPASS 

.50 

EVANESCENT 



5 TO 7 5 GHz H 2 


1 CONTIGUOUS 

DIPLEXER | 


CONTIGUOUS 
L DIPLEXER _J 




BANDPASS 

50 ii 1 

EVANESCENT 

FTl 


, 1 r 

1 BANDPASS 


j .VANESCENT 

1 1 l 


TENTATIVE BANDPASS FILTER 

DESIGNS l_ 

N - 5. 0.01 dB RIPPLE TCHEBYSHEV 
3 dB BANDWIDTH - CHANNEL BANDWIDTH 
SIGNAL - 95% MAIN CHANNEL 

5% ADJACENT CHANNELS 


LOW PASS 
» 1* 2.5 

c 

■h 

r 


i 

CONTIGUOUS 

diplex e ^ | 

i 

i oo n i 



1 


CONTIGUOUS 
DIPLEXER 


BANDPASS 

1 | 50 “4 

EVANESCENT 

l 


■O- 


BANDPASS 

1 1 50 n r 

EVANESCENT 

1 l 


1 


iT* 

CHANNEL 

F.F' dB 

T °K 

t if- k 

GOAL 

FOR K - 2 

H 1 

5.5 

738 

0.2 

0.18 

H-2 

5.1 

648 

0.2 

0 16 

H 3 

5.3 

692 

0.2 

0 24 

H 4 

4.9 

606 

0.2 

0.22 

H 5 

4.3 

490 

0.2 

0 19 

H 6 

3.7 

389 

0.2 

0 24 


CONTIGUOUS 
L . DIPLEXER 1 

MIXER CONVERSION LOSS 
FEED LOSS - 0 5 dB t35°Kl 
ANTENNA TEMPERATURE 


- 8 5 dB 

- 250°K 


AT - RADIOMETER SENSITIVITY 


IF 


- IF NOISE FIGURE 

- IF TEMPERATURE 


FIGURE ’83 GHz IF NETWORK DESIGN 


1.75 TO 5 GHz H-3 
■ T - 5.3 dB 


2.5 TO 3 75 GHz H 4 
F T - 4 9 dB 


1 25 TO 2 5 GHz H 5 
F T - 4 3 dB 


0 625 TO 1.24 GHz H 6 
F T - 3 7 dB 
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4. Receiver Design 
4. 3 IF Design 

4.?. 2 NETWORK FOR 20 to 4000 Mil.*. IF 

Satisfactory IF multiplexing and noise figure performance for the 11 temp erature 
channels can be provided bv direct input to a filter network, but a broadband prf - 
amplifier between the mixer and filter could be beneficial to isolate the mixet 
from the filter network. 


The optimum performance of a mixer is achieved with an optimum termination 
for the mixer at the IF frequency. Broadband Ga As FI. 1 devices are presentK 
available with meet the need of a first stage IF amplifier from 20 to 4000 MHz. 

The following discussion treats two cases: lithe use of a triplexer filter directly 

connected to'the mixer and 2) use of a Ga As FET amplifier as a butter between 
the mixer and filter. In Figure 1 the input to the triplexer is shown with a tenta- 
tive 100 ohm input impedance. Directly below, an optional amplifier is shown. 

The noise figure of the amplifier at the high end of the band is 4. a dB; at the low 
end it is about 2.:' dB. The signal enters the triplexer filter assembly and is 
divided into three bands: 3040 to 40°0 MHz. 1580 to 2280 MHz, and 20 to 
1340 MM/. The lower band is amplified with a 40 dB amplifier and entois a 
wav isolated resistive power divider having a 23 dB loss. The nine outputs are 
separately filtered with bandpass coax filters, amplitied, and then detect ed 1 he 
effective noise figures of each channel F if is given at the output oi each amplifier. 
In parenthesis, a larger noise figure is indicated for the two high frequency chan- 
nels. The noise figure in parenthesis is that seen when the preamphnei is 
inserted between the mixer and the triplexer; without the amplifier, the noise tig- 
ure is determined bv the 3.0 dB IF amplifiers for each of these channels plus the 
effects of losses introduced bv the filter. However, when the amplifier is inserted, 
the noise figure is dominated by the front-end amplifier. Thus, the use ot a 
broadband front-end amplifier is justified only as a buffer to aid in tno broadband 
impedance match between the mixer and IF, and is achieved onlv at tne expense 
of the noise figure of the higher frequency channels. 

Calculations of radiometer AT performance for ea:h of the T channels is given in 
Figure 1. The IF noise figures assumed are those obtained without the broadband 
amplifier. The mixer conversion loss assumed is 3 dB, a value which IS con- 
sidered within thi state of the art at 116 GHz. An antenna feed loss oi 1 dB and an 
antenna temperature of 250°K are assumed. It should be noted that the calculated 
AT values are lower than the values listed as goals. lhe conclusion ol this 
assessment is that the 118 GHz receiver technology is well within the state oi the 

art. 
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[ .. AT* 

CHANNEL 

F| F . dB 

T °K 
1 IF* * 

GOAL 

FOR K - 2 

T 1 

3.3 

330 

0.2 

0.0*5 

T-2 

3.3 

330 

0.2 

0.08 

T-3 

3.0 

359 

0.2 

0.11 

T-4 

3.5 

359 

3.2 

0.14 

T-5 

3.5 

359 

0.2 

0.16 

T-6 

3.0 

290 

0.2 

0.18 

T-7 

3.0 

290 

0.2 

0.22 

T-a 

3.0 

290 

0.2 

0.27 

T -9 

2.5 

225 

0.2 

0.34 

T-10 

2.5 

225 

0.2 

0.40 

T-1 1 

2.5 

225 

0.2 

0.40 


AT - RADIOMETER SENSITIVITY 
F |F - IF NOISE FIGURE 
T |p - IF TEMPERATURE 


MIXER CONVERSION LOSS - 5 dB 

FEED LOSS - 1 dB 

ANTENNA TEMPERATURE - 250°K 


FIGURE 1. 118 GHz IF NETWORK DESIGN 


4. Receiver Design 

4.6 Quasi-optical Diplexers 

4.6. 1 SEPARATION OF 183 GHz BAND 

Wire grid polarizers are the lowest loss elements capable 

wave bands. It is logical to use such a device to separate the 18 * GHz band tio 
the other bands. 


Th P most difficult radiometric band is the 183 GHz band. High mixer conversion 
loss difficult IF, and local oscillator source noise and power limitations combine 
to produce a receiver noise temperature significantly higher than any of the other 
radiometric bands. Input feed losses add further to the input noise temperature. 

Since the 183 GHz band is a key MA.SR sounding band, one approach is to i minimi 
input feed losses in this band even at the expense of increased losses m the other 

bands. 

If it was desirable to minimize losses in the 1 18 GHz band at the expense of the 
183 GHz band, the 118 GHz band would be separated using a polarizer an ^ h feed 
183 GHz band would be removed with a high pass filter resulting . 
loss at 1 18 GHz and 1.0 dB loss at 183 GHz. The resulting A 1 s would be those 
presented in Topic 2. 1, improved at 118 GHz and degraded at 183 GHz by a factor 

of 1. 12 

A wire grid polarizer functions as a diplexer in a radiometer becauS / s b ° t '\ t P i ° 1 e ar ' 

Wire grid polarizers have long been in common use as microwave devices. Imple- 
mentation is routine at the lower microwave frequencies but « ^ ethod 

183 GHz. Two fabrication techniques are commonly used. The p t ffj) F - 2 

consists of photo- etched copper strips on a mylar of Kapton substrata. § 

shows the desion for a midband frequency of 150 GHz. with reactive and ohmic 

losses estimated for 183 GHz. Losses in the mylar or Kapton ^T^OOzTnch) ’ 

tut are kept low because the thickness of the sheet is only 50 pm (0. 002 men). 

Another fabrication method is by hand winding phosphor bronze wire on screws of 
very small pitch, 100 to 250 turns per inch. The tightlv stretched gri o 
formedby the screw pitch is bonded to an elliptical ring, and the excess wire is 
cut ally after the bond is set. The resulting polarizer will provide a circular 
aperture at an incident angle of 45°. Figure 3 illustrates winding the : me wire on 
the screws, bonding to the ring, and finishing the fina, polarizer. Perioimance 
determined by the wire size, spacing, and wavelength, as shown in Figui e . 

Refer ence s 

1 T. S. Chu, M. J. Cans, W. E. Legg, "Quasi-optical Polarization Duplexing 
' Of Microwaves," BSTJ, Vol. 54, No. 10, December 1975. 

2. J. Payne, private communication, National Radio Astronomy Observatory, 

Kitt Peak, Tucson, Arizona, September 1977. 
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183 GHz 
RECEIVER 


ALL BANDS 


ALL 

BANDS 



11 BANDS 


140, 118 
AND 104 GHz 
RECEIVERS 


FIGURE 1. BASIC WIRE GRID POLARIZER 


y- 

-<#>. 


MYLAR 
OR KAPTON 


DESIGN PARAMETERS 

0 - 45° ANGLE OF INCIDENCE 
b - 0.1 mm GRATING PERIOD 
d = 0.04 mm COPPER STRIP WIDTH 
t = 0.005 mm COPPER STRIP THICKNESS 
t'- 0.05 mm MYLAR THICKNESS 
REACTIVE LOSSES 

T n < 0.001 (0.004 dB) 

R1 < 0.001 (0.004 dB) 

OHMIC LOSSES 

T < 0.003 (0.01 dB) 

Rl< 0.001 (0.004 dB) 



SMALL PITCH 
SCREWS 


'PHOSPHOR 

BRONZE 

WIRE 


BOND TO 
ELLIPTICAL 
RING v 




FINAL POLARIZER AFTER 
CUTTING AWAY WIRE 


FIGURE 2. PH0T0ETCH FABRICATION METHOD 


FIGURE 3. WIRE MESH FABRICATION METHOD 
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4. Receiver Design 

4.6 Quasi-optical Diplexers 

4. 6. 2 SEPARATION OF 118 GHz BAND FROM WINDOW BANDS USING 

fabry-perot filter 

minimize geometric walkoff losses. - “ 


The characteristics of Fabry-Perot filte ^ achieved . (1) For some applica- 

mission bands and broad reflection ban for other] such as for a high pass 

tions, the Fabry-Perot is an ideal diple , suited to 4) A anique application 

or low pass application, the device is p Y * ^ u g GHz band from the window 

is seen for the Fabry-Perot in separating the and 140 are about 4 GHz, 

bands at 104 and 140 GHz The require b^andwito^t W an^ Subsection 3 . 8 , a 

permitting a loaded Q of about 30. g , n f on lv 0 1 dB and reflection 

Filter can be designed which has a transmission loss of only 0 dB u mted 

loss of 0. 14 dB The is reflected at a„ angle of 30*. The 

with the incident axis and the 1 . d an d fourth half-wave resonances 

unique property of the filter is tha close to the specified win- 

(M = 3, M = 4) are resonant at frequencies ; that very clo^ ^ ^ s H ssumed , 

dow frequencies. For example if a re P ^ £ourth g is at 14 0 GHz. Since the 
frequencies 5 are not critical, the 105 GHz ^uency would serve 

». - ...» » 

at an etalon tilt angle of 0 - 15 . 

The reflectivity of the etalon need be only °’ 7 ' fYs GHzband, thU amounts 

118 GHz of 15 dB. In terms of reflection loss lor the no un 

to 0 . 14 dB. 

* j • a i An 3 Tiprtur 6 of 70 mm is assumed 

" ufraction losses are neglu 

gible because the plate spacing is only 4. 2 mm. 
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4-44 




4. Receiver Design 

4. 6 Quasi-optical Diplexers 


4. 6. 3 


l 18 GHz BAND from window band using 

RESONANT GRID DICHROIC FILTERS 


band's “om^htn's OH, 1 : es °. nan ‘g rid diplexers can be used to separate the window 
at 45 “ incMence devices have the advantage of operating well 


”ated with P iloi / * dlchr01c fBter ls a thin sheet of conducting material perfo- 

sheets^resonant It Ll™ ™ ^ desired fret i uenc y band. Two or more 

to produce shTroer n, 3ame / r ® qUenC y ° r at separate frequencies may also be used 
‘° pr h ° sharper passband characteristics. Excellent design approximations 

can be obtained from a mathematical model based on a planar array of idenHcal 

SiEL 

^°de mod^s 1 aL^^lpe 2 “ V ° khl f the continuiVof thetogenHaf flelX”' theVave- 

— «=» -sssr^aK * 

J^O^rnrb^h^l 011 a i d re 1 flection characteristics of a resonant grid backed by a 
0 . 002 inch thick material of dielectric con^nf f ~ 7 . i * S » . X a 

ure 2 riparlir • « m^eciric constant € r =3.1 (Kapton) are shown in Fie- 

S £Sr‘~ - “ 

Grating lobe formation is avoided by selecting DY < ?X // 1 j. cir. a\ u , . , 

=';rr g ; 

addition of identical sheets. P may G made stee Per, however, by the 
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5. ANTENNA INTERFACE 


Interfaces between the radiometer receiver and antenna consist of 
mechanical/thermal, electrical, feed, and calibration. Special requirements of 
the feed interface include illumination with a quasi- optical quadraplexer. 


One of the primary tasks of MASK system engineering is to establish interface 
specifications between the antenna and radiometer receiver subsystems. Such 
interfaces as mechanical /thermal and electrical, including bus power and telem- 
etry and command, are routinely achieved during system definition phases of 
flight programs. Establishing those interfaces at this time in the concept phase 
is premature. However, for the present study two critical areas require pre- 
liminary interfaces to be established: the feed interface and the antenna calibra- 
tion interface. 

Heed Interface. The concept of using beam optics design formalism in the receiver 
subsystem has assumed the use of Gaussian illumination to facilitate design pro- 
cedures and to minimize diffraction losses. This concept implies that illumination 
patterns are generated in the receiver subsystem and are the responsibility of the 
receiver subsystem engineering. Thus, contrary to a conventional antenna inter- 
face such as a waveguide flange, the MASR antenna interface consists of a set of 
four complex beam parameters at the focus of the antenna. 

Corrugated horns with quarter-wave spacings carrying the HEj \ mode produce 
nearly symmetric circular Gaussian field amplitude patterns, Gaussian design 
formalism is therefore appropriate in quasi-optical feeds for MASR. Separate 
horns may optimally generate the pattern for each frequency band, and these 
bands can be combined with a quasi-optical quadraplexer. Topics 5. 1, 5. 2, and 
5. 3 address the details of the feed interface. 

Antenna Calibration. The radiometer receiver measures antenna temperature 
continuously with a AT resolution of typically- 0. 2 0 and an absolute temperature 
accuracy of ±1 °K. Antenna characteristics must be known to relate the radio- 
metrically measured antenna temperature to the desired main beam radiation 
temperature. Brief*/, the calibration procedure consists of three steps. A "cold" 
look into space will yield a residual antenna temperature caused by ohmic losses. 

A relative antenna pattern measurement will allow observations to be corrected 
for sidelobc introduced errors. Finally, comparison between radiosonde and 
MASR temperature measurements will yield an accurate absolute calibration. 

Topic 5,4 addresses the details of the antenna interface. 




****** Nor 




5-3 


5. Antenna Interface 


5.1 ILLUMINATION PATTERNS 

Corrugated horns with quarter-wave spacings carrying the HEn mode produce 
neai ly symmetric circular Gaussian field amplitude patterns. Gaussian optical 
design formalism is therefore appropriate in quasi-optical feeds for the MASR 
antenna. 


By convention, an antenna pattern of a receiving antenna is defined by reciprocity 
to be identical to the pattern generated by a transmitter having the identical 
antenna configuration. The MASR receiver antenna pattern is thus a function of 
the illumination pattern of the feed horn. 

For an antenna of infinite extent, a Gaussian illumination profile produces a far 
field pattern which has no sidelobe structure. The main beam field intensity 
decreases monotonically off axis in a characteristic Gaussian pattern. In prac- 
tice, a Gaussian illumination pattern is truncated by the antenna edge which intro- 
duces sidelobe s in the far field. Illumination patterns generated by computers 
through iterative processes to obtain minimum sidelobes are very nearly Gaussian 
profiles. 

The radiation from a horn can be made polarization symmetric through the use of 
corrugations spaced \/4 apart near the radiating aperature. d) The field pattern 
of the HEj 2 mode in the region of the horn is then given by^> 

■ Vo [HP] 

(r and a are defined in Figure 1). 

The Gaussian approximation of this function is 

g(r) = E Q exp (_r 2 /w 0 2 ) 

where w Q is the beam waist when the field is e" 1 times the field amplitude on axis 
and r is the variable distance from the z-axis. If we let these two functions be 
equal at the e -i point, we have the relations r = w Q 

and 


w 

2.405 — 2 = 1. 7538 

3 . 

w = 0.73a 
o 


Thus, 


g(r) = E q exp (-r 2 /0. 53a 2 ) 

Both the Bessel and Gaussian functions are plotted in the Figure 1. It can be seen 
that very close agreement is achieved by the proper choice of Gaussian profile. 
The tail of the Gaussian is truncated at r = a whereas the Bessel function goes 
to zero at the walls. This inconsistency is the greatest anamoly in the Gaussian 
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appr°xittiati°n of the field in the region of the horn. However, in the far field of 
the horn, for large values of z, very close approximations to Gaussian have been 
demonstrated. Side lobe levels of -40 dB and antenna beam efficiencies of 
98 percent are common using these methods. The close approximation of actual 
optimized fields m high performance antennas to Gaussian fields make the use of 
Gaussian optical design formalism not only warranted but actually preferred to 
other design approaches. 
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FIGURE 1. SCALAR HORN ILLUMINATION PATTERN COMPARED TO GAUSSIAN 
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5, Antenna Interface 


5. 2 MASR ANTENNA AND FEED GEOMETRY 

The recommended antenna configuration for MASR is a symmetrical Cassegrain 
reflector system. The main reflector truncates the Gaussian illumination at the 
99 percent beam energy points, and forward spillover is limited to -32 dB down 
from the total beam energy. 


A detailed tradeoff study for the MASR antenna system has determined that the 
symmetrical cassegrain antenna is lighter, more compact, and has a better pat- 
tern for given surface figure than the offset paraboloidal reflector. I 1 1 The basic 
antenna layout is shown in Figure 1. The primary reflector is a 4. 4 meter diam- 
eter parabolic reflector which has a prime focus at 1.54 meter from the center 
of the dish. The secondary hyperboloidal reflector is disposed between the pri- 
mary reflector and focus such that an external focus is formed 0. 16 meter behind 
the surface of the antenna. The external focus is the antenna feed point for the 
four radiometer frequency bands, 104 GHz being the lowest and 183 GHz being the 
highest frequency. 

Since the antenna feed point is located 1. 64 meter from the center of the subreflec- 
tor, the energy distribution characteristics at the external focus will be determined 
by diffraction. 

A special property of Gaussian beams is that the relative intensity measured at a 
point is equal to the fractional power for all radii beyond that point. In other 
words, the relative intensity at a point at the edge of the subreflector is equal to 
the fractional spillover. In the presently recommended design, the illumination 
of the primary is controlled by a maximum subreflector diameter of -0. 198 meter, 
while the spillover is limited by a total subreflector diameter of 0. 25 meter. 

Thus, the relative intensity at 0. 198 meter is 

~T = ex p (^i) = 0,01 

o \ w / 

N z ' 


where 


0. 198 
r = 2 

Solving for w we have 
z 

w =0. 065 m 
z 

Having determined that the e“* illumination circle of points is 6. 5 cm from the 
center of the center of the subreflector', the spillover or energy which falls beyond 
the 0. 25 meter outer diameter is 



or about -32 dB down from the total power. It has been determined that spillover 
of these levels does not seriously affect the radiometer performance (Topic 5. 3). 
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consists of four coaxiaU^ams^ ° f ca , sse S rain antenna. The feed point 
band has a beam diameter of 63 7mm to the radiometric bands. The 104 GHz 
183 GHz band has a beam diame‘ttr ™“ 39 5 I] ? tenslt J. P° ints - “<J the 

diffraction ’ ^he^in ^ 0 ' 3 ^ 11 re ^ uirem ^“t of 1 tl^subreflertoi^^^!) 1 &e r iaws e of 0 ^ er ” 

S a Sue -r^tr^nfo 1 1 i n Fig r e . 2 - ***• 

at these limits is 1 oercenf- qo * * * 1 /i th f b eam limits. The intensity 

these limits. Also note that the mumbnltfo ° f ^ beam e ? ergy is contained within 
frequencies Thus the t lllu ” 1 ™ atlon pattern requirements holds for all 

Gaussian formalism descrLd° n°TopTc ' 3 7 °^he fi^ d ^ d f. t ® rmined fromt ^ 

Of \V Z and the distance F from fV, 1 ? « , Th f d condlhons are tb e values 

tion half angle 6 subrefle ^°r to the focus which fixes the diffrac- 


w 


1/2 


= 0. 040 rad 


From Topic 3.7 we have the relation 


1/2 


= 

TTW 


or w = — 
o 7T0 


ofWm dtm^te^TtThe^ ^ Tba -lues 

the beam waist values by 3 . 035. SV P ° are obtained b V multiplying 

^ th e "antlnna^cu s. h “"aMe I'ust^talL^of b giVS beam diamet -s 

confocal parameters for each of th to t * bea ” 7 waist > beam diameters, and 
eter values arTeTsentUlTn d ° w T freqU f nCy ba »ds. The confocal param- 
described in Topic 5 3 desl g nin 8 tba quasi-optical quadraplexer which is 

Reference 

K FeallbUi^StuTv’ At ™ s P heric Sounding Radiometer Antenna 
easibiiity btudy, Final Report, Contract NAS5- 24087, December 1977 . 


TABLE 1. BEAM PARAMETERS FOR FOUR RADIOMETER BANDS 


Frequency, 
f, GHz 

Wavelength, 

A, mm 

Half Angle, 
01 /2» rad 

Beam Waist, 
w 0 , mm 

Beam Diameter 
at 1% Intensity, 
mm 

Confocal 
Parameter, 
b, mm 

183 

140 

118 

104 

1.64 

2.14 

2.54 

2.88 

0.040 

0.040 

0.040 

0.040 

13.05 

17.03 

20.2 

22.9 

39.5 

51.6 
61.2 
63.4 

326 

425 

505 

572 

1 

1 


ILLUMINATION 

AMPLITUDE 


1.0 


0.5 


0.1 


LIMITS 
0.198 m 



FIGURE 2. ILLUMINATION PROFILE REQUIREMENT FOR ALL FREQUENCIES 


183 GHz 140 GHz 118 GHz 104 GHz 



FIGURE 3. BEAM DIAMETERS AT WAIST AT EXTERNAL FOCUS 
FOR 99 PERCENT TRUNCATION 
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5. Antenna Interface 

5.3 QUASI-OPTICAL QU ADR APLEXING 

f ^ d h ° rnS f ° r fOU1 ' radiometric bands can be imaged to common 
feed point through quasi-optical techniques. Feed losses for iq^ r*w i a 

can be kept below 0. 5 dB fnd below 1 dB for th^ woz st band 83 GHZ 


sasr- 

vsn fSspx^sst as ”rv.ir,s: rss? s& 

point, to i effect images to the 183 GHz and 1 18 GHz feed horns. Lens images or 
i elay images are formed for the 104 and 140 GHz bands. ® 

The first mirror in the system was determined to be a wire grid polarizer for 
good reasons: it is the lowest loss quasi-optical element whfch i S P capable of 
separating two beams (Topic 4.6); and it takes advantage of the energy in both 

Estimated feed losses for each band are listed in Table 1. Losses for the 18^ PR* 

A horn S lTsTofO e TsYb fo balanced ^lixer which also uses a wire grfd. 

" t the throughput loss^f two” denticLTo^ns operating Iga^t'a^a^boUc" 
mroi a- conjugate foci. The total loss of two horns at 140 GHz was CL 5 dB. 


TABLE 1. ESTIMATED FEED LOSSES FOR FOUR RADIOMETER BANDS 


1 

Band 


Wire 

Grid* 

Fabry- 
Perot 2 

Fabry- 
Perot 3 

Resonant 
Grid 2 

Resonant 
Grid 3 

Horn 
Loss, dB 

Total 
Loss, dB 

* •:! 

w 

w 

M 

183 

140 

118 

104 

<0.1 

<0.1 

<0.1 

<0.1 

<0.1 

<0.1 

<0.1 

<0.1 

0.1** 

0.14 

o.r 

0.1** 

H 

0.1** 

0.1** 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

<0.45 

<0.55 

<0.95 

<0.5 

<0.85 

<0.55 

<0.95 

■ "‘i 

01 ' ! 

H ; 

R 

i 

i 

Assumed using quasi-optical balanced mixer. 





**Midband loss, ADD 1 dB at band edge. 
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5.4 ANTENNA CALIBRATION 

Antenna performance must be known to relate the r adiometrically measured 
antenna temperature to the desired main beam radiation temperature. 


Antenna Temperature. The measured antenna temperature at the radiometer 
calibration reference plane for each radiometric channel is related to the radia- 
tion brightness temperature in all directions by the equation 


= 'iff D(e, 


4>)T B (0, <t>)dS2 + (1 - k 


) T 

r o 


where 


T^ = antenna temperature 

k^ = antenna ohmic efficiency 

D = antenna directivity 

Tg - radiation brightness temperature 

= polar and azimuthal angles about the antenna 
df2 = element of solid angLe 

T q = average physical temperature of the antenna reflecting surfaces, 
weighted by the feed system illumination pattern 

The antenna directivity is defined as the ratio of the power intensity produced by 
the antenna in the given direction to the mean value of the intensitv in all directions 

The first term of equation 1 is the contribution to the antenna temperature from 
energy received by the antenna, reduced by the ohmic efficiency. The second term 
is due to radiation emitted by losses present in the antenna. 

There are two sources of antenna resistive loss. The first is due to small losses 

present in the antenna reflecting surfaces, and the second is due to the portion of 

subreflector blockage energy which is terminated in the mixer. Both of thee s 

losses contribute to k . 

r 

An alternate representation of this resistive loss due to blockage is obtained by 
considering the mixer diode as a thermal radiator at a temperature determined by 
diode characteristics and operating point. Some of the thermal energy emitted by 
the mixer will be reflected back into the feed by the subreflector, raising the 
antenna temperature (see section 5. 5). 

D(8, tj 3 ) is composed of a main beam, near sidelobes, and remote sideLobes. Near 
sidelobes are defined as those which fall on the disc of the earth (Tg = 250°K) and 




remote sidelobes are directed to cold space (Tg = 2. 5 °K). 
temperature can now be expressed as 


The antenna 


T 


A = Tm'Vlm + T ns k r r 'ns + T rs k r \ s + ' k )T 


where 


r o 


( 2 ) 


n m ” P ortion of the pattern energy in the main beam 

= ~^D(0,<|.)dn 

main beam 


and where 


m 


- weighted average main beam radiation temperature 


4irr| 


m 


D (e.4>) t (e.^jdn 


main beam 


an^ per cent* ntributio ns ^ S ’ the and rem0te sidelobe temperatures 


r 'm + 1 lns + ’l 


rs 


= 1 . 


S^b^htn^Tfl 1 Equat , ion 2 is the contribution to the antenna temperature due to 
the brightness temperature m the main beam. The second term represents enerev 
coUected by antenna near sidelobes from the eath. The third ter m P accounts for 87 
far sidelobe energy from cold space and the fourth term is due to energy emitted 
by the antenna reflecting surfaces. The energy collected from the sun £ the 
antenna remote sidelobes is negligible. 

aC , Ur f tely determine main beam brightnes s temperature T m from Ta 
other factors affecting antenna temperature must be known. This Requires A> 
knowledge of the antenna performance characteristics. 

§“^£ 7 ■ Pr ° C edur £ ■ Range measurement of the MASK antenna is impractical 
due to the large range length required to satisfy the far field criterion. ? 


R > 


2d* 

X 


For tL R MAq h R * g 1 * g ?o* d 13 the a P erture diameter, and X is the wavelength 
rr ^ h f MASK antenna at 183 GHz, the minimum range length is about 15 mills * 

Ground reflections and atmospheric effects would degrade measurement acnirarv 
However, the antenna performance can be measured after launch. aQy * 

Briefly, the calibration procedure consists of: 1) a ’’cold look*' infn 

main antenna system; 2 ) an approximate antenna pattern ^metur'men^ ^ « 

atmospheric temperature measurement by both radiosonde and MASR. 
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The cold look is used to determine the antenna ohmic efficiency k r , The pattern 
measurement yields an estimate of the ratio of r\ m to rj ns and an energy distri- 
bution in the near sidelobes, which is used to correct MASR measurements for 
sidelobe contributions in a relative manner. Finally, comparison between MASR 
and radiosonde measurements enables the total amount of pattern energy on the 
disc of the earth (rj m + q ns ) to be accurately determined. The calibration is 
complete . 

1) Cold Look . With the antenna looking at cold space, the first term in 
Equation 1 is near zero, which exposes the antenna resistive loss radiation 
(l-k r )T 0 . To minimize collected energy the earth, moon, and sun should be 
located behind the antenna. This posture could be accomplished by positioning the 
spacecraft looking away from the earth. Reflecting surface temperatures needed 
to calculate T 0 can be obtained from an antenna thermal model, allowing the deter- 
mination of kr» ' 

2) Pattern Measurement , the antenna pattern is needed over an angle 
somewhat larger than the disc of the earth to allow the calculation of energy col- 
lected by antenna sidelobes from the earth during observations. The overall 
system temperature accuracy is affected very little by the accuracy of the mea- 
sured antenna pattern, and an approximate pattern is all that is required. The 
pattern can be measured in orbit using either a ground-based transmitter or 
astronomical sources. 

A ground transmitter measurement is simple in concept and has the capability of 
wide dynamic range with accurate control of transmitted power. Atmospheric 
absorption in the oxygen and water vapor high frequency IF channels is less than 
10 dB and should not be a problem. The transmitted signal can be modulated to 
improve sensitivity against the thermal background of the earth. 

However, the demands placed on the antenna positioning system by a transmitter 
measurement are severe. To map the main beam using a transmitter would 
require an antenna positioning accuracy of about one-tenth beam width or 0. 05 mr 
at 183 GHz. The presently planned 0. 15 mr would not be sufficient. Also, due to 
the fine structure of the antenna pattern, a large number of data points would be 
required. 

Pattern measurement using the moon and the sun as millimeter wave sources 
eliminates the disadvantages of a ground-based transmitter measurement while 
yielding results accurate enough for our use. An averaged pattern, rather than a 
point-by-point mapping, can be obtained by deconvolving scans across the moon 
with antenna patterns containing various amounts and distributions of sidelobe 
energy. Constraints on a correctly deconvolved scan are that no energy come from 
beyond the disc of the moon and that there be no limb brightening at the dark edge. 
Average sidelobe levels to 45 dB below the main beam can be mapped using the 
moon. The sun can be used to extend this measurement to -60 dB. 

This procedure is described in Reference 1 where it was used to measure the pat- 
tern of the National Radio Astronomy Observatory 11 meter millimeter -wave tele- 
scope operating at 214 GHz. The pattern information obtained was used to correct 
solar maps for sidelobe introduced errors. 

The orbital motion of the satellite can be used to scan the antenna beam across the 
moon, eliminating the problems with antenna positioning inherent in a ground 
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transmitter measurement. The orbital angular rate is 0.07 mr/sec, or about 
7 s econds /beamwidth at 183 GHz, which is a convenient rate. 

Pattern measurement using large sources such as the moon and the sun does not 
yieLd information about the fine structure of the antenna pattern, but an averaged 
pattern is alL that is required to make an estimate of the temperature contribution 
of the sidelobe energy during earth observations. 

3) Radiosonde Comparison . The amount of the total antenna energy which 
falls on the earth ( T l m + r lns) rnust be known, and can be determined by comparing 
radiometric and direct atmospheric measurements. The ratio of r) m to q ns is 
known from the pattern measurements. This ratio, along with observation over 
the disc of the earth, allows the radiometric observations in the area of the radio- 
sonde measurement to be corrected for relative sidelobe contributions. At this 
point, the estimate of (rjrn + r lns) adjusted to bring the radiometric measure- 
ment^ into^agreement with the radiosonde measurements. This procedure assumes 
that T m = T ns , which is accurate. 


The amount of energy on the earth outside the main beam will be a function of beam 
position on the earth. The antenna patterns 'are used to calculate changes inri ns 
and rj rs for various beam positions. 

Temperature Accuracy Estimate . To calculate an estimate of MASR system tem- 
perature accuracy, nominal values of antenna performance and projected tempera- 
tures for the worst case at 183 GHz will be used to determine the accuracy of the 
calibration procedure previously discussed. Then the uncertainty in system per- 
formance parameters obtained during calibration will be used to calculate overall 
system temperature accuracy. 


1} Antenna Performance Estimate (183 GHz) 

Power Breakdown, % 


Resistive losses 1 

Central blockage, 2% 1 

1 

Spar blockage 2 

Spillover 1 

Random scattering, 5% 4 

1 

Sidelobe s, 4% 2 

2 . 

Main beam 8 5 


Temperature 
antenna To 
earth 

antenna To 

space _ 

earth 

earth 

space 

earth 

space 

earth 
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Performance Parameters (nominal values) 

k = 0,98 (from resistive loss and central blockage) 

T* 

r] = 0.8674 
1 m 


■n = 0, 08 16 
'ns 

■n = 0.0510 
'rs 


Z) Cold Look . To estimate the error in the --ement of (l~k r ), Equa- 

tion 2 is solved for the term 

(i-k r ) - =-[T A -(T m k r ^ m + T ns k r n ns + T rs k r n rs )] 


Since 


T 

T = T 

m 

ns rs 


1 r 

<i-V = 

f - t a' ( t . 


r = T^l A' ^signal'] 

Using nominal values and estimated errors of Tq , Ta > and ^signal* the 1 cr error 
in (1 - kr) is the rss (root sum of the squares) of the errors due to each component. 


Error 

Vi 

1 

*— < 
< 

T = 200 ± 20 °K 
0 

0. - 02 

t a +at a = t a ± °- 7 ° k 

0. 0035 

T +AT = T . . ± 0.4 °K 

s s signal 

0. 002 

( 1 - k } = 0.02 ± 0.0045 
r 



To correctly compute the expected value for Ta> the radiometric temperature 
scale, not the depressed radiometric temperature scale, should be used to add the 
contributions from T Q and T s i gna L (Topic 3.5). At 183 GHz, cold space emits 
very little radiation, and Ta is comprised largely of antenna resistive loss 
radiation. 

3) Pa ttern Measurement . It is anticipated that the pattern measurement 
using the moon and the sun will be accurate enough to determine the ratio to 
ri^g to within 20 percent. 


= 10,6 ±2 




4) Radiosonde Measurement, Using the radiosonde__data and the 
atmospheric model, the main beam radiation temperature T m can be predicted to 
within 1 °K._ This^ measurement will be conducted in an area of slowly varying tem- 
perature, T m - T ns £ 10 °K. By making observations around the test point, T ns 
can be determined to within 5 °K. 

The purpose of the joint MASR/r adiosonde observation is to determine the sum 
*lm + T )ns> fc he percentage of the pattern en ergy_f ailing on the disc of the earth. To 
solve Equation for q m + rj n ^, the values of T m and T ns are combined into an 
average earth temperature, T e , known to 1 °K due to the smalL effect of the uncer- 
tainty in T ns* 


+ \s 


T A- T rs k l n rs -( 1 -k r )T o 

T e k r 


Summing the errors as before 


Error 

A(r lm +n ns ) 

T a = 245.5 ± 0.7 °K 

0. 0036 

T = 4.5° ± 0. 5 °K 
rs 

0. 0001 

T = 200 ° ± 20 °K 
o 

0. 0020 

~T q = 259.5 ± 1 °K 

0. 0037 

k r = 0.98 ± 0.0045 

0. 0008 

q + q = 0, 949 ± 0. 0056 

m 'ns 



The calibration is complete, 

5) Overall System Accuracy /, To determine the uncertainty in T during 
observations. Equation 2 is solved for this quantity. 111 


m 


k r\ 
r 'm 


T.-T kn -T k n - T (1 - k 
A ns r 'ns r s r 1 rs o v r 


Nominal values a_nd uncertainties are known for ail terms on the right hand side, 
and the error in T m can be found as before. However, because the sum 

+ Hns + ^rs = * without error, there are only two independent errors involving 
the three q, which are the error in the ratio q^^to q ns and the error in the sum 
+ Ins* Also, the difference between T m and T ns is likely to be greater during 
observations than during calibration. 
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For example, 


160.69° 


! |245 o - 240°(0. 98){0. 0816) - 4. 2 e (0. 98)(0. 05 10) - 200°(0. 

(0.98)(0.8674)r 



Error 

AT m (°K) 

Kn^ns' 

1,65° 

Rrni* hns^ 

0.45° 

k 

0.24° 

r 


t a 

0. 82 c 

T 

0.47° 

ns 


T 

o 

0.47° 


The total system temperature accuracy is found to be ± 2 . 0 K. 
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5. Antenna Interface 
5 5 SYSTEMATIC ERRORS 

l “ r P° Sltl0n JT" mmatch between the antenna and the feed. 


ec- 


*i.* efaWiitv These constraints can oe eascu uauig optical design proce- 

»t. -a an, 

»400 °K due to nonthermal noise generation in the “ , chopper wheel reflect- 


Af = 


. . ■, r isrr-ut- arid r\ is the two way reflection path length. Foi the 

MASK antenna th/peHod of the antenna subreflector standing wave is 

88 MHz. 

Shown in Figure 1 is a sy ^ temper 

Astronomy Observatory li me ^ _ temperature variations with two dis- 

sideband spectrometer at 88. 6 GH . S J riation ^ a riod of 22. 8 MHz is due to 
tinct periods are evident. The rap larger variation with a period 

the reflection from the antenna subreflector , arnM ftew ^ absorbing surface, 
of 280 MHz is due to the reflection from the ch pp “ reflections will 

The system temperature o[ e ^h ™dione£ic channel £ ^ctTem-’ 

during calibration at the end X th ‘ len eth changes between calibrations, 

perature accuracy only if the re .P , IF “xhis effort is most noticeable in 

shifting tne phase of the standing wave xn the IF ^effort » ^ 

the narrow IF channels which only c v P f is not likely to vary, and this 

distance from the feed horns to the chopper SU |L a “ "fte reflection path length 
reflection will not degrade system ac “,.A el ^ co „ tl .oll e d.’ With a subreflector mis- 
to the mam subreflector must be NR a O telescope, the subreflector position 

^hVXeS^t^n changed TmU. in the rest of the system 

■w ill be less. 

Computer programs developed 

C._,.W available «£ » ACCOS- «£<=« »• 

rently used in the design of dn ec e , liaht The use of these powerful 

srgrK ^rr^s^^u^..., P = apP =.. p 
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1) MISSION REQUIREMENTS 

The missionof the microwave atmospheric sounding radiometer 
(MASR) is to collect data to aid in the observation* and prediction of severe 
storms and other mesoscale weather phenomena. The MASR is to be launched 
via shuttle into a geostationary altitude orbit. This orbit choice allows the con- 
tinuous measurement needed to resolve mesoscale dynamics. The instrument 
will operate in conjunction with multispectral imaging and infrared sounding 
sensors, and its function is to provide three dimensional temperature and humid 
ity fields, primarily in cloudy and overcast atmospheric regions. Further 
mention of the optical sensors will be limited to noting that the other sensors 
have stringent line of sight pointing requirements. This will limit allowable 
spacecraft attitude perturbations due to reaction torques of scanning the MASR. 
A 20 minute pitch and roll stability of 11 prad is typical of the requirements 
of the advanced atmospheric sounding and imaging radiometer (AASIR). 

As part of an IR&D program, Hughes has been assessing the observa- 
tional requirements for satellite data for all scales of weather phenomena. An 
understanding of the data requirements for mesoscale weather prediction has 
been derived from numerous discussions with NASA, NOAA, military, and 
university personnel, as well as from the literature in the field. The require- 
ments listed below represent Hughes current understanding of how the MASR 
might be used in a severe thunderstorm environment, but obviously represent 
only one of a number of scenarios for the use of the MASR. This material is 
included to indicate Hughes continuing interest in mesoscale forecasting 
problems. 

The observation requirements for mesoscale weather phenom na are 
as varied as the phenomena themselves. A particularly simple mesoscale 
phenomenon to predict is frost, in which knowledge of the synoptic conditions, 
surface temperature, and existing and projected cloud cover is adequate to 
forecast frost a few hours in advance. At the other extreme are the dota 
requirements for the severe thunderstorm environment, where knowledge of 
the synoptic environment, coupled with frequent measurements from the 
synoptic scale down to the cloud droplet scale, are needed. The severe 
thunderstorm has been chosen to size the MASR mission requirements, both 
because it places severe time constraints on frame repetition rates and 
because of the economic and social benefits to be derived from the improved 
prediction of this phenomenon. 
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What are the measurement requirement simposed^by^ which measure - 

thunderstorm environment? The time an P erstorm development. In 

ments are needed depend upon the »tate and convection have begun, 

the early stages of the storm, bef ® ti scale resolution (=400 km 

knowledge of the temperature field with synopti humidit y field and loca- 

in the horizontal and =8 levels in :he '" (~ 2 0 to 30 km in the horizontal 

tion of the clouds with mesoscale resolut ( enc o{ once per 2 hours, 

and =10 vertical levels), with a measurement f q convergence has begun 
are the measurements of should be measured with 

and the inversion layer begins to rl **' *“ hour frequency. Finally, from the 
mesoscale resolution and with a O P the end of the storm, measure- 
beginning of cumulous cell e p Ucal resolutions should be made 

ments with mesoscale he rizontal a ^ time , measurements with sub 

approximately every 15 minutes. At t ■ rates would be desirable for 

mesoscale resolution and m . ln ex tent of storrrt development, and for f°r e " 
precise location ° f r e s elution from satellite data will be limite 

■rthe^nef""; t"^ infrared imagery data. ^ 

The use of the MASK should be Umitec Ito km 2 ln 30 min- 

partly cloudy atmospheric regions Its cov a measurement over the 

utes means that it would take nearly nredict i 0 n model, which might be 
2500 km2 extent of a typical P f oba t!e storm activity in the 

The coverage area of the MASR pm^n^pro^eds’throughout 

repetition rate increased as thunderstorm P t the MASR will be 

the day. After the beginning of c "“^ 0 ^ d d k umidity under the cloud top level, 
most useful in measuring temperature and humm: Required to sound beneath 
During the height of the storm provide semiquantitative pre- 

the thunderstorm's cirrus anvil TTl^nder storm oirer an area of 500 km? with 
cipitation data in the center of ' the th accuracy requirements for the 

a 15 minute repetition rate. Estimates of . o V K £ , the temperature 

temperature and .humidity field, 
field and approximately 10 pei cent ior 

£ TvifAQR z^ata it is likely there will 
In addition to the quantitative us® 8 ° .... maD s of these data for visual 
be a desire to produce temperature and hum ld tty ™ p would fae tQ 

analysis and interpretation. * P 0 *"^ Areas 

map total vertically integrated water P &te areas Q f convergence and 

where the water vapor is increasing Another application of this 

likely locations for f "* ur ® 00 ? V tb C e temperature and humidity sounding data to 
type would be to combine both P conventional stability 

produce a stability map using one or m°r areas of convective instability. 

,b ' MASR “ 

mission objectives is given m Table 2-1. 


TABLE 2-1. MASR MISSION REQUIREMENTS 


a 

Spatial resolution 



Horizontal 

20 to 30 km 


Vertical 

8 to 1 0 levels up to 1 00 mb 

9 

Coverage rate 

750 km 2 /30 min 

9 

Temperature accuracy 

lo to 20 K 

9 

Humidity accuracy 

10 to 20% absolute 


2) SYSTEM REQUIREMENTS 

The horizontal spatial resolution requirements for temperature and 
humidity sounding of severe storms vary for different severe storms and meso- 
scale phenomena, but most can be satisfied with the 20 to 30 km resolution 
element indicated in Table 2-1, For nadir observation from synchronous alti- 
tude, this range corresponds to antenna beam widths in the range 0, 56 to 
0. 84 mr, or 0. 032° to 0. 048°. 

Due to the fundamental limitation of diffraction, a typical radiometer 
antenna of diameter D can only give a beamwidth of the order 1.3 X/D, where 
X is a wavelength. The frequencies of the microwave absorption resonances 
of O 2 and H 2 O, which are available for the sounding of atmospheric tempera- 
ture and water vapor, respectively, are shown in Figure 2-1. The 60 GHz 
O 2 and 22 GHz H 2 O resonances have been used in microwave sounders flown 
or designed to date, since these are for observing synoptic scale phenomena 
from low altitudes (-1000 km), where antenna beamwidths in the range 7° to 
10° are adequate. For geosynchronous sounding of severe storms, there is 
a need to use the higher frequency O 2 resonance at 118 GHz (2.53 mm) for 
temperature sounding, and the H 2 O resonance at 183 GHz (1.64 mm) for 
water vapor sounding. These shorter wavelengths allow the required narrow 
beamwidths to be achievable with practical antenna sizes (<4.4 m) that can 
fit into the shuttle payload bay without furling. 

The need for a narrow antenna beam thus leads to two critical require- 
ments for geosynchronous microwave sounding: 

® Large aperture reflector antenna, 3 to 4.4 m in diameter 

© Low noise radiometer operating between 118 and 183 GHz 

To map a finite area on the earth* s surface, the antenna beam has to scan in 
two orthogonal directions. Because of the very high frequencies utilized and 
the requirement for low loss reception, electrical scanning is impractical to 
implement and mechanical scanning is required. This leads to a third criti- 
cal requirement: 

© Mechanical scanning of a large antenna in two orthogonal axes. 
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This last requirement is subject to the constraints that the antenna 
scan motion must not cause excessive reaction torques on the spacecraft that 
would violate the pointing accuracy and stability requirements of other 
instruments (e.g. f AASIR). 

These three requirements determine the general characteristics of a 
microwave sounding system needed for observing severe storms from 
geosynchronous altitude. The objective of the two proposed studies is to 
verify the feasibility of developing such a system for a future flight in the 
1982-1983 time frame. 

3) MASK REQUIREMENTS 

IF Bandwidth and Channel Selection 

. ^ ° 1 ’der to better appreciate the basis for the system r equirements 

involved in the selection of operating frequencies and IF bandwidths, it is 
necessary to review briefly the principle of microwave sounding and radiative 
transfer. 

The MASR is a multichannel radiometer (or spectrometer) for perform- 
ing microwave sounding experiments from a geosynchronous satellite. Micro - 
wave sounding consists of measuring, at selected frequencies, the brightness 
temperature of passive microwave radiation emitted by the earth 1 s atmosphere 
and surface, and the derivation of meteorological parameters from these 
measurements. The basic equation that relates the brightness temperature 
o± the earth's outgoing radiation to the atmospheric parameters is the micro- 
wave radiative transfer equation. 
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where 

= brightness temper atxire of earth's outgoing microwave 
radiation as a function of frequency f 

T{z)= atmospheric temperature as a function of altitude z 



FIGURE 2-2. MICROWAVE SOUNDING SYSTEM FUNCTIONAL FLOW DIAGRAM 
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where G(f,0 , cb) is the antenna gain function, (6,<J>) are polar and azimut a 
angles about the antenna. The above equation together with the radiative 
transfer Equation 1 form the system equations shown in Figure . 

Rewriting Equation 1 in a more compact form 


“* T (f ) 


x B (f) = S T s e 


W(f, z) T (z) dz 


where 


T (f ) = I C*(f, z) dz 


is the zenith opacity and W(f, z), which depends on a{ f,z), is a weighting 
faction that weighs" the contribution to T B from the atmospheric temperature 
at various altitudes. For the case where the surface reflectance is negligi 


f 03 

_j o:{i, z)dz 


W(f, z) = a{ f, z)e 


The weighting function usually peaks at some altitude and^th is 
with frequency. In temperature sounding, ) . . , . 

crete frequencies f, about an 0 2 absorption line, each value of £; being 
selected for the desired altitude where the weighting function Wffi, z) will 

peak. 

A consideration in choosing the set of W(f L , z) is that their mutual 
overlap be minimized, since overlapping leads to 

different Tutfi) measurements and is responsible for the instability in ir.,eit 
ing Equation 3 (Reference 2). To minimize overlap, the half-width of each 
wffi, z) should be minimized, and the frequencies f L be chosen to space the 
W(f.,z) over the altitude range to be covered. 

At each RF, what is measured is the average brightness temperature 
T b over a spectral band Af = B IF . The appropriate weighting function for 

T„ is a band averaged weighting function 
B 


Wff. z) 


= 4. f w(f. 



1 


In general, as Af increases, W(f, z) becomes broader in half-width. There- 
fore, it is desirable to keep Af or Bjf small enough so that W(f, z) will 
approach the monochromatic W(f, z), which has the minimum width. Hence, 
the maximum value of Bif for the temperature sounding channels is limited 
by the need to maintain vertical resolution for the sounder. 

pit T ° f aCh , IF / n T \ ble 1 of the radiometer RFP, there corresponds two 
RF symmetrical about the 118.75 GHz 0 2 line center. To each RF frequency, 
there corresponds a weighting function W( f, z) that can be calculated from the 
known oxygen absorption characteristics about 118.75 GHz and the known oxv- 
gen concentration profile in the atmosphere. However, if *(f,z) [ s assumed 

which "I ^Td me t T m K Ut 118 * 75 GHZ> th6n T(f) and W < f ’ z) are symmetrical, 

which leads to T B (f) being symmetrical as well. Each MASR temperature 
sounding channel will observe the average of the brightness temperatures at 
3 18 . "/IF and 118.75 + f IF , and there is effectively only one weighting 
function for each of the T-channels. The values H ave shown in Table 1 of 
the radiometer RFP indicate the average height of each channel weighting 
function. Channels T-l to T-4 are designed to probe the troposphere, while 
the other channels are to probe the lower stratosphere. Figure 2-3 shows 
some calculated weighting functions on the low frequency side of the 
118. 75 GHz line, which are taken from Reference 3. A few of these weight - 
mg functions are at frequencies approximately the same as some of the lower 
sideband frequencies in Table 1 of the RFP, and the average heights of the 
weig. ting functions in Figure 2-3 generally agree with the corresponding 
average heights given in Table 2-1. 

Factors Affecting System Performance 

The Figure 2-4 diagram summarizes how the various system require- 
ments interact and lead to the radiometer performance needed. The require- 
ments of area coverage,, temporal and (horizontal) spatial resolution of severe 
storm phenomena, combine together to determine the dwell time available 
for sounding each FOV cell. Contiguous coverage of an area L, km by L km 
within a frame time tF, with a nadir ground resolution size of X km, leads 
to an .FOV dwell time, or integration time 


t . 

1 x . 2 F 


is formula is approximate in that it neglects the time to turn around, step, 
and calibrate at the end of each scan, and it assumes that the L x L area is 
located in the subsatellite region. Contiguous coverage of the 183 GHz chan- 
nels is required to generate water vapor maps. This requirement leads to an 
integration time of (ti)ig 3 - 1 second as specified in Table 2 of the RFP. 

This value is obtained from the above equation with X = (35, 800 km) 

(1. 3 X./D), if we substitute in the maximum antenna diameter of 4.4 m, and 
Ij - 750 km, t F = 30 min as typical mission requirements (Table 2-1). The 


F‘-V 
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integration time and beamwidth for the 183 GHz channel determine a scan 
rate 0 SC , which rate the antenna scan mechanism must be designed to 
move the 183 GHz beam. But the same mechanism will move the other 
channel beams at the same rate, so the FOV dwell times for the other chan- 
nels are simply determined by (t^ = (1.3 \/D)/e sc , or (tj*. is proportional 
to Hence, ( fc ; ) 1 1 8 183/118 = 1.55 times (ts)io-a, as indicated in Tables 1 

and 2 of the RFP. 1 

As indicated in Figure 2-4, the requirements of temperature and spec- 
tral resolution determine the required AT and Bjp, respectively. These 
requirements differ for different types of sounding experiments. For temper- 
ature profile sounding, the accuracy required in measuring Tg(f) has to be 
typically many times the accuracy being sought for the temperature profile. 
This is because the mathematical "inversion" techniques used to invert the 
integral Equation 1 to recover T{z) are highly sensitive to data noise, as the 
computations deal with inversion of ill-conditioned matrices. For example, 
in Tiros-N, where temperature profiles accurate to 1° to 1. 5°K are sought, 
the microwave sounding unit operating at the 60 GHz O 2 band has a 
AT = 0. 3°K, Similarly, NEMS at 60 GHz had a AT = 0. 2° to 0. 3°K, and 
was designed to retrieve atmospheric temperatures to accuracies of 1° to 
3 K. The random errors in Tg(f) that can be tolerated is also dependent on 
the particular inversion algorithm to be used. The ATg accuracy needed can 
be found through "computational experiments" by simulating the inversion 
process with artificial noise injected into the data. From the ATg accuracy, 
one can specify the noise equivalent AT^ desired. 

The selection of Bjp for each temperature sounding channel is a com- 
v f j‘ :rnise between the desire of increasing Bjjr to reduce AT versus that of 
*! *c reusing Bjp to provide sharper weighting functions. In temperature 
-ot; wring, spectral resolution is related to vertical spatial resolution; so the 
Bfp tradeoff is between temperature and vertical resolution. 

In water vapor sounding, the required AT is controlled by the -10 to 
20 percent accuracy desired in measuring the water vapor profile n wv (z). 

The question of what random errors in the set of 183 GHz Tg(f^)'s will cor- 
respond to 10 to 20 percent error in n wv (z^}'s is dependent on the numerical 
algorithm used to go from the set of TBffjJ's to the set of n wv (z^) 1 s. A way 
to answer this is again to perform computational experiments on the actual 
algorithm to be used, by injecting artificial noise into the water T„ data 
vector. ^ 
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APPENDIX B. REFERENCE AVERAGING IN PRESENCE OF GAIN FLUCTUATIONS 


1, Introduction 

The temperature resolution of a microwave radiometer is determined by 
errors from two sources: (1) random fluctuations inherent in the square- 

law detection of noise power in a frequency band; and (2) temporal variations 
in gain. The use of "reference averaging" to improve a Dicke radiometer, 
as first proposed by Bremer and Prince, ^ has the effect of reducing the 
first type of errors. The technique reduces the variance in the reference 
voltage by increasing its integration time far postdetection smoothing, 
which leads to a decrease in the overall variance of the differential output 
voltage. In a switching radiometer with 50% duty cycle, instead of integrating 
over one-half the dwell time (T^) for a resolution element or "pixel", we can 
integrate over Nt^/ 2 by averaging the reference samples in N consecutive 
pixels. Such an extension of integration time is possible for the reference 
samples because the reference temperature is essentially constant. On the 
other hand, the integration time for the antenna samples must be limited to 
one-haif the pixel dwell time in order to resolve the variation of scene 
temperature from pixel to pixel. 

As the reference averaging time Nt^/ 2 increases, the temperature error 
due to detector fluctuations decreases, which corresponds to a decrease in 
the K-factor. If the chopping fraction remains at 50%, it can be shown that 
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K = J 2 + | . Mathematically , as N — K -^2, which corresponds to 

about 30% reduction in K from its original value of 2. However, physically, 
it is not clear how large a value of N can or should be used, which would 
determine how much improvement is physically realizable. Before this 
question can be answered, it is necessary first to understand how reference 
averaging affects the second type of errors; in particular, to know also 
how the errors due to gain variations depend on the parameter H. This paper 
is aimed at clarifying these questions. Our results will show that the 
gain variation error increases with increasing S, so there is, in fact, an 
optimal value of N that minimizes the combination of detector fluctuation 

and gain variation errors . 

Above, we have described reference averaging with symmetric chopping. 

A further reduction in detector noise is possible by incorporating asym- 
metric chopping with reference averaging. If the fraction of time for 
antenna sampling is increased from .5 to some value q (.5 < q < 1.0), the 
antenna integration time is lengthened from .5-r p to qf p , while the referen 
integration time in each pixel is shortened from .5T p to (l-q)T p . However, 
the latter effect is counteracted by averaging the references over H pixels 
to give a total reference integration time of N(l-q)T p . With a large 
enough H, the overall effect is to reduce both the variances of the reference 

and antenna voltages. The K-factor now becomes a function of « and q, and 

( 1 ) _ . . 


is expressed by 


f J 1 , For a given N, K is minimum when 

K “ / q + N(l-q) 


. If q is set equal to q 


for each N value, then as 



N is increased, the K- factor decreases as K = 1 + 


. Thus, as N **♦ 00 , 


u-j 
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K -♦ 1, which is the value for a total power radiometer. As in the symmetric 
chopping case, one might again inquire how large a value of N can be used, 
in an attempt to ascertain the maximum improvement achievable. However, in 
any practical situation, there will always be gain variabilities, in which 
case the question of maximum improvement must be reformulated differently. 

We already know the detector fluctuation error as a function of N and q, 
but we also need to know the gain variation error as a function of N and q, 
so we can determine the values of N and q that will minimize the combined 
error due to both effects. This will then give the maximum resolution 
improvement achievable by reference averaging with asymmetric chopping. 

In this paper, we will start with analyzing only the case of reference 
averaging with symmetric chopping. Our approach is to start with a conven- 
tional Dicke system with symmetric square-wave chopping and square-wave 
detection, and examine what happens to the gain variation error as reference 
averaging is incorporated with increasing values of N. The chopping fraction 
is kept constant at 50%, while N is varied. This facilitates comparison with 
the conventional Dicke case, and allow;? easy interpretation of the results. 
The effect of adding asymmetric chopping to reference averaging will be con- 
sidered in later studies. However, some of the characteristics we will 
derive for the symmetric chopping case will be qualitatively similar to that 
for asymmetric chopping, though the details will differ; for example, the 
combined error due to detector fluctuation and gain variation when plotted 
against N will show a minimum, for any constant value of q, similar to wtat 
we derive for q = .5. 


The effect of gain variabilities on a reference averaging radiometer 
was^ discussed in Reference 1, but was not treated' with adequate depth to 







B-3 




reveal the tradeoff between gain variation and detector noise errors. 
Consequently, there was some lack of appreciation of the limitation on the 

improvement realizable from reference averaging. 

( 2 ) 

Ross Graves has studied the problem of reference averaging in the 
presence of gain variation as a statistical estimation problem; i.e., as a 
problem in obtaining a best estimate to the true value of the time-varying 
reference voltage at a specified time (the observed pixel time), given data 
samples of the reference voltage at other times which are contaminated by 
the detector noise. He treated the generalized problem of using N pixels of 
reference samples to estimate an m-parameter reference-voltage function of 
time (m < N) , and studied also the generalized technique of reference 
averaging with non-uniform weighting of the samples. His study utilized 
very elegant techniques to arrive at a number of interesting results, and 
the reader is referred to Reference 2 for a discussion of his numerous results. 

In this paper, we will examine the reference averaging system in its 
originally proposed form, where the one-parameter estimator is a constant, 
equal to the unweighted, simple average of the reference samples over N pixels. 
''Symmetric" averaging is assumed, where the N consecutive pixels of reference 
samples to be averaged are time-symmetrically distributed about the center of 
the observed pixel dwell time. Noting that gain variations are independent of 
detector fluctuations, we may set the latter to zero while we analyze the 
errors due solely to gain variations. This is the reverse of setting the gain 
variations to zero while analyzing the errors due solely to detector fluctua- 
tions, as was done in deriving the K- factor improvement for ideal reference 
averaging. The errors derived for these two opposite cases may be combined 
by RSS to yield the radiometer's temperature resolution. Since the detector 
f luct*’' .on noise is "turned off", we may leave aside the statistical 
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estimation aspects of the problem for the time being, and concentrate on 
studying the signal processing operations implied by the reference averaging 
algorithm. We want to find out how the radiometer output response to the 
antenna temperature signal is modified by gain variation with time which may 
be deterministic and/or random; and how these modifications depend on the 
reference averaging system parameters, such as N, t ? and the chopping period. 
The abjective is to relate antenna temperature error to gain variabilities, 
and to examine the dependence of that relationship on the parameters of 
reference averaging. Another objective is to derive equations that can be 
used for future computer simulation of reference averaging systems in 
general, to study system parameter tradeoffs and error sensitivities. 













2 * Response of Dicke vs. Reference Averaging Radiometer 

To compare the reference averaging and Dicke systems, it is help- 
ful to cast both signal processing schemes into a common format, so 
that they can be readily compared. Toward this end, we view both 
systems as having postdetection sections equivalent to the block diagram 
of Figure 1. The modulated output from the square-law detector, denoted 
by V(t), is multiplied by a reference waveform C(t) that is synchronized 
with the modulation frequency and phase. The product C(t)V(t) is 
filtered and smoothed by a pure integrator, to yield the radiometer 
output V Q . The multiplier reference waveform C(t) is sometimes referred 
to as the "correlation” in Dicke radiometry. (3) It will be shown that 
by properly defining the correlation C(t) a digital reference averaging 
system can be made equivalent to the system of Figure 1. The function 
C(t) is different for different values of N, and we will study the 
sequence of C N (t) 's. 

The modulated output from the square-law detector is given by 

V(t) = Bkg(t)['X(t) + Tjjjj] + n(t) (X) 

where B is the IF or predetection bandwidth, k is Boltzmann's constant, 
g(t) xs ’’e power gain, is the receiver noise temperature, T(t) is the 
noise temperature at the switching terminal, and „(t) is the zero-mean random 
process representing the fluctuations in the square-law detector. Later on, 
n(t) will be set equal to zero, when we analyze the errors due solely to gain 
variations. The notation is simplified by defining the product 

G(t) = Bkg (t) 
and rewriting V(t) as 

VCt) = G(t)[T(t) + Trk ] + n(t) (3) 
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we will refer to G<M as the "gain", with the understanding that it is really 
the product of the electronic gain and kB; so variations in G(t) can arise 
from variations in the electronic gain and/or bandwidth. Therefore, the 
error formulas to be derived for variations in G(t.l can be applied to 
either of these two ef.-ects, or a combination of both. 


Consider a Dicke radiometer with 50% chopping fraction and a chopping 

period t o . T(t) will be a square-wave function of time that alternately 

equals the antenna temperature T A and reference temperature T R every t c /2 

seconds. In general, T & is a function of time as the radiometer beam 

scans the scene, while T R is essentially constant, though it may have 

some residual variation depending on how well the reference source 

temperature is stabilized. To simplify the analysis, the radiometer beam 

will be assumed to move from pixel to pixel in an idealized step-dwell 

mode, with the stepping time being negligible compared to the dwell time. 

Dnder this assumption, a typical plot of T(t.) might look like that shown 

on Figure 2(a), The labeled intervals along the top denote pixel dwell 

time slots, each having a width equal to the pixel dwell time T For 

P 

simplicity of illustration, the waveform is drawn for the case t = t /2 

c p # • 

T a is constant within each pixel dwell time, but its value will generally 

change from pixel to pixel following the scene brightness variation. The 

values of and T R within the j-th pixel time are denoted by T A> and T R 

respectively. The corresponding values of V(t) when T(t) equals T^ and*T R 

will be denoted by V. and V_ , respectively. 

J j 

Referring to Figure (1) , the radiometer output V q has the general 

form 



/■ 


= constant x JC(t) V(t)dt 
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In what follows, the correlation C(t) for the Dicke and reference 
averaging systems will be determined and compared. In the Dicke case, 
the radiometer output that indicates the brightness temperature of the 
i-th pixel (see Figure 2) is simply given by 

V ‘ C5) 

o k t R. 

where \ and V R are, respectively, the average of the antenna and 

reference segments of V(t) , both taken within the i-th pixel dwell 

time. Combining the integrals in V A and V R _ , we can also express V Q as 

3 * . 


V « — f c. (t)v(t)dt 

° T pX 


( 6 ) 


where (^(t) is a square-wave function that alternates between +1 and 

-1 every t c /2 seconds, as shcn 
outside the pixel dwell time. 


-1 every rj2 seconds, as shown in Figure 2(b), and is equal to aero 


We will use the symbol C N (t) to denote the correlation for reference 
averaging over H pixels. Since the Dicke is a special case of the 
reference averaging system with N = 1, the Dicke correlation is denoted 
by C x (t). We will find that the C^t) functions are easier to describe 
by sketching them, rather than writing out their algebraic equations; 
therefore, we will sketch out each function in defining it. 

In discussing reference averaging, we will always assume "symmetric" 
averaging, where the N consecutive pixels of reference samples to be 
averaged are time-symmetrically distributed about the center of the observed 
pixel time slot. It is expected that a digital reference averaging system 



will most likely use symmetric averaging, since this can be implemented 

easily with a small amount of memory and can yeild great advantage in 

reducing the system sensitivity to gain variation. For symmetric 

N-l 

averaging, if N is odd, we will take an integral number (— ) of 
reference pixels before and after the observed pixel time slot. But 
if N is even, we will take an integral number plus half of a reference 
pixel on either side of the observed pixel time. 

Consider first the specific case of reference averaging with N = 3. 
The radiometer output that indicates the brightness temperature of the 
i-th pixel is computed by 




+ V 0 + V. 

K.. 

1 


R, 


J 


(7) 


where V A is the same as previously defined and V R ^ is the average of 

A i .3 

the reference voltages in the j-th pixel. CDmbining the integrals 

from all four terms in (7) , it can also be expressed as 


V o = ~ J C 3 (t)V(t)dt 
3t p 

where Cg(t) is the square-wave function shown in Figure 2(c). It can be 
described as a train of rectangular pulses of width t c /2 and heights 
equal to +1, -j , or zero, distributed in a symmetric pattern over a 
time interval of 3 t^. Outside this interval, Cg(t) — 0, so its overall 
width is 3 

Now, consider the case of reference averaging with arbitrary value of 
N. Assuming that N is odd, the radiometer output for the i-th pixel is 


given by 


V 

o 


m 


V - i V V n 

A i N — R iHj 
j= -m 
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where m - l£ we again combine the integrals from all the terms, 

Cm 

we cnn write 


7 

o 



C N (t)V(t)dt 



( 10 ) 


where C N (t) for N>3 is illustrated in Figure 2(d). Comparing C 3 (t) 

and C N (t) for N>3, we note the following: The positive pulses in both 

functions are located inside the central pixel time slot, and their 

pulse heights remain equal to unity. The height of the negative pulses 

in C^Ct) are reduced to a value - while tha number of negative pulses 

is increased, being spread out over a larger Interval. The overall width 

of C (t) is increased to a value NT . C (t) may be viewed as a train of 
N P ** 

T 1 

rectangular pulses of width and heights equal to + 1 , j or zero, 
symmetrically distributed over a time interval of NT p 5 outside this 
interval, C^Ct) “ 0, 

C N (t) depends not only on N, but also on the chopping period and 
pixel dwell time. We may indicate this by writing C N (t) = C N (t; Tc> T p ), 
Substituting Equation (3) into (10), we obtain 


V 0 - f I C N (t;T c ,T p )G(t)[T(t) + T^dt 
NT 

P 

+ r L ~[ V t;T c’ V n(t)dt (11) 

P ‘'Nt 

P 

This generalized system equation can be used to set up a future computer 
simulation of reference averaging radiometers, to study parameter trade- 
offs and error sensitivities . It describes the linear response of the 




radiometer output to the modulated temperature signal, shows how this 
response is disturbed by gain variations and detector fluctuations, 
and indicates how the response is modified by changes in the design 
parameters N. t q and t . 
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3, A n tenna Te mperature E rrors due to Gain Variability 

In what follows, we will set r.(t) to zero, and analyze the system 
behavior in the presence of only gain variations. With n(t) = 0, 
Equation (IX) reads 


V = ^ 



Vt)G(t)[T(t) + T^Jdt 


cu 


Consider the above equation for the Dicke case, for which N = 1. Within 

the pixel dwell time, T. and T are corje;f*;*nf- • cn 

A u are constant, so T(t) inside the interval 

of integration may be represented by 

TCt )=i(T A - Wt)+f(i A + y (13 

where C^t) is the Dicke correlation shown in Figure 2(b). Substituting 
(13) into (12) *?ith N = 1, we have 



Defining G as tne average gain within the pixel i 


pixel interval. 


G * ~ / G (t)dt 


CIS) 


and y 1 as the integral 


Y l m r^J c iCt)G(t)dt 

T 


(16) 
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and making use of the fact that C^Ct) = X within the pixel interval, 
we may rewrite (14) as 

V Q = G(T a - V + Yl + Tjjjj) 

The second term in (17) is usually small in a conventional Dicke, but 
will not necessarily remain small with reference averaging. Therefore, 
we will retain this term throughout our analysis and will evaluate and 
compare it for the Dicke and reference averaging case. 

If the second term were neglected, as is usually done, we would 
arrive at the following classic equation for the gain variation error 
in a Dicke radiometer: 


( 17 ) 





( 18 ) 


Evidently, the gain referred to in this familiar formula should be the 
average gain G within the observed pixel dwell time, and 6G is the 
deviation of the actual G at the time of observation from its value at 
the last calibration. 


With the second term present, any variations in G(t) will, in 
general, cause changes in G and Yl via Equations (15) and (16). The 
complete equation for the antenna temperature! error caused by gain 
variations should, therefore, be as follows: 



fiG(X A - Tr ) + 6 Yl 



Hr 

2 



where 6G has the same meaning. as above, and 6y^ is the variation in 


( 19 ) 
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Now, consider Equation (12) for the case of reference 


averaging over 

N pixels. It will be convenient to divide the integration in Equation (12) 


T 

into two parts; integration inside the observed pixel where |t~t J <— £, 

T o 2 

plus integration outside the observed pixel where It-t I >—£. t is the 

1 o l 2 o 


time at the center of the observed pixel interval, 

T 

In the region |t-t | <vp-, T(t) can be represented by Equation (13) 


which we rewrite as 


T(t) - (T a - T r ) [% G x (t) + %] + T r 

where is the antenna temperature in the observed pixel, and C 1 (t) 

- T 
xs the Dicke correlation. In the region |t-t |< -£• 

1 o' 2 


c N ct>:% c x (t) + %] « % o^t) + % 


Because of this property, the central part of the integral in (12) is 
given by 


2 _ f 

T *l 


; ■# 


C N (t)G(t) 


_ T ( fc ) + yj 


dt 




lh C x (t) + %]G(t)dt 


+(t r + V 




G N (t)G(t)dt 




+ d R 



C (t)G(t)dt ( 



where G and Yl are previously defined in (15) and 


(16) 


In the region |t-t l> 

o' 2 9 


X(t) = 


l A. 

,1 


“R 


during antenna ha If -cycles 
during reference half-cycles 


""" ,A ] ** j-tt pixel other thee tt. 

observed pixel. , _ , r „ ot £ „ ^ ^ 


But since, for it-t l> -2. 

1 o f 2 


C N Ct) N l 


N 


during antenna half-cycles 
during reference half-cycles 


the following relation holds: 


Therefore, the remaining portion of the i 


integral in Equation (12) is 


■/ 


C N (t)G(t) 


t-t > -2 
o 1 2 


[ T(t) + V,] 


dt 


( 23 ) 


(24) 


(25) 


t r + T m)r 


C N (t)G(t)dt 


t-t > -£• 
o 2 


(26) 
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Addins (22) and (26), we find that 


V 

o 


/ C N (t)G(t) 
Nt 

P 


[ T(t) + T Rxn] 


dt 



where we define 


r N 



C N (t)G(t)dt 


If we rewrite Equation (17) for the Dicke case as follows: 

V o ' ( G + t) ( T A ' + Y X ^T r + T ffl j (17) 

we see that (27) also describes the Dicke case when y N - v Ihe slope 
of the linear relation between V Q and is the same for the Dicke as 
for the reference averaging case, but the intercept is different. 

Taking the differential of Equation (27), we obtain the following 
expression for the antenna temperature error due to gain variations in 
a reference averaging radiometer: 


6 V - 




+ T„„V1 


Vi 



where SG, 6y^ and 6y^ are the variations in the gain dependent quantities 

G, y x and y N from their values at the last calibration. The above formula 

for N - 1 applies to the Dicke radiometer, and is the same as Equation (19). 

One should distinguish between the two terms appearing above. The 

first term is multiplied by (T^ - T R ) , and can therefore be suppressed by 

keeping T R close to T^. We may call it a "differential radiometer type" 

gain variation error. The second term, on the other hand, is multiplied 

by T r + T^, the system noise temperature during the reference look. It 

behaves like a "total power radiometer type" gain variation error, and is 

a more sensitive term. Moreover, it contains the dependence of 6T on 

A 

the parameter N. In particular, the difference between the 6T for a 

A 

Dicke and 5T^ for reference averaging over N pixels lies entirely in the 
second term. Therefore, our attention will be focused on analyzing the 
behavior of the second terms in Equations (27; and (29), 
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4 ‘ g ain Variation Error Dependence o n N and Chopping Fr^ncngy 

Up to this point, our formalism is applicable to arbitrary G(t) 
functions. For any given model of G(t), we first evaluate G, y and 
y N from Equations (15), (16) and (28), then substitute these inputs 
into Equation (27) to find the V q response; and then determine the 
antenna temperature error from Equation (29) by considering the 
variation of these inputs. This analytical procedure can be readily 

coded in a computer program, to handle G(t) functions of arbitrary 
complexity . 

In the following analysis, we apply our general formalism to a 

simple model for G(t), which can describe low frequency gain changes 

approximately. Our objective is to investigate the dependence of the 

gain variation error on the parameters N, t and t , and to obtain 

c p 

formulas that can be used to find order-of -magnitude estimates of the 
increase in gain variation error traded for reduction of detector noise 
error when reference averaging is used. 

We will assume that the variations in G(t) within an interval Nt 
are slow enough that they can be approximated by a cubic equation of the 
form 

3 

V — ^ 

G(t) a n (t - V n 

n=0 


where a n ere constants, and t Q is the time at the center of the observed 
pixel interval, which is also the center of the reference averaging 
interval. The range of N values for which this approximation is valid 
will depend on T p and on the frequencies of the significant components 
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in the gain variation spectrum. Should this cubic model be inadequate 
to cover the interested range of N values for the gain variation fre- 
quencies encountered, one could repeat this analysis with higher order 
models for G(t), 

In the analysis to follow, integrals involving C N (t) will be 
evaluated, and the following assumptions will be made: 

T 

m That -*■ = n is an integer, n being the number of chop 

v y T C C 

c 

cycles per pixel time. Consequently, 


Cjj(t) x constant dt = 0 


This equality follows from the fact that C N (t) has n Q positive pulses of 
height + 1 and Nn, negative pulses of height - all the pulses being 
of equal widths, 

(2) That Cjy(t) is a symmetric function of t-t Q , where t Q is the 
center of the observed pixel dwell time and also the center of the 
integration range. This requires, first of all, symmetric reference 
averaging, which we have assumed all along. In addition, it requires 
letting the pixel-to-pixel transition time fall, not at an antenna- 
reference chopper transition time, but a quarter chop cycle later or 
before, as shown in Figure 2. The same phasing also applies to the 
end points of the total integration range. The reference averaging 
algorithm should take in the appropriate fraction of a chop on each 
end of the integration range to make C N (t) symmetric about t Q . As a 
consequence of this symmetry, • 



C N (t) • (t - t o ) dt = 0 


J 

NT 

p 

for n equal to an odd integer. 

T 

(3) That -£■ << 1 , which allows us to use certain approximations 
P r 

to evaluate integrals involving C w (t). In practice, -—may range from 

W P 

to y-j , depending on chopper design. 

Substituting the G(t) given in (30) into Equation (28), 



V 0 <t-t o ) 2 dt 


since the a , a , and a terms do not contribute because of (31) and 
o 1- J 

(32). Changing the integration variable to t-t Q , and noting the 
symmetry of the integrand, the above expression may be rewritten as 


Y N = 



NT 72 
P 


c N (t)t dt 


For N = 1, 


Vi 



C, (t)t Z dt 






t! 


where the integral has been evaluated approximately by a technique similar to 
that of Attachment 1. The negative value for the integral assumes that C^t) 
starts out at t = 0 with H, which would be the case if there are an 
Odd number of chop cycles per pixel. If there are an even number of 
Chop cycles per pixel, 0 l( t) would start out at t = 0 with - 1 , and 
the sign of the integral would be positive, but the magnitude remains 

the same. 

, , r. a s the pulse-train function of period 

In Attachment 1, we defined pC ) ^ 

T , which alternately equals +1 and 0 every -j- seconds. p(t) is ^ 

il Figure A-l of Attachment 1. Identifying the pulse width. At = T , ^Ct) 

can be expressed in terms of p(t) as follows: 


P(t) 


N 


(*-> 


% 

3 < t < -f 


(36) 


c H (t) 




t > 
u ^ 2 


Therefore, the 


inte graX that appears in (34) is gxven by 


NT 

—E. 

2 


,2 


nt 
— 2 . 
2 


2 1 

C N (t)t 2 dt = p(t)t dt 


f .(.-*) 


t 2 dt 


0 


0 


■r 


T 

je. 

2 


NT„ T 
•+ 


I— E.+-— 


P (t)t dt 


-if 


p (t) t A dt 





In Attachment 1, we proved the following approximation, 
t. 


/ 


‘b .3 3.2 2 

9 - t, t, - t /T 

p(t)t 2 dt « — — 3 b a 


ft) 


( 38 ) 


By this formula. 


T 


(39) 


'•T? 


/ 


NT T 

—£ J. 

2 2 


p(t)t dt = 


N' 


+ £ • r -A( 


T \ 

_E_ 


3 ' 2 t pA 16 


(40) 


Substituting these into (37) gives 


/ 


Nt 
— E 
2 


C N (t)t‘dt = 


N 2 -l N+l . T c 


T 3 

. JE_ 
16 


(41) 


Equation (34) becomes 


Y = - ¥-£- i + N+l . V 
Y N 4 l 3 2 T 

\ P, 


( 42 ) 


Returning to Equation (29) , let us distinguish that part of 6T, 


which depends on N by the symbol 6T. 


( 43 ) 


8<? a - \ 6Y n ' I t r + t ru 
G+- 


' ( t r + t rh) 


8y h / . ' 

* G \ 5 i T V 


where has been neglected in comparison with G because 
6 


Yl 


a t 2 T 
2 p n c 

2 


8 T 



p 


< < G 


(44) 


which follows from -=■ < < 1, and from the expected upper bound on gain 

T P 

variability of, say, 

v2 



a 

2 \ 2 


< few percent 


(45) 


Combining (42) and (43), we find that 6T^ for reference averaging is 
given by 

6 $ a _ 6 Ve! + +T ) 

6T A 4G \ 3 + 2 


(46) 


which shows explicitly the dependence on the parameters N and — , the 

.. - P 

latter being inversely proportional to chopping frequency. Compare this 

with 6T. for a Diche radiometer, which is given by 
A 


A 8a„T " /T 


(47) 
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and is reduced to a 


A Q 

In the Dicke case, 6T. is proportional to — - , 

A p 

small value by keeping ^ < < 1. In the reference averaging case, 61^ 

P 

T 

c 

contains a term independent of — which rapidly increases with increasing N. 

P 

The difference between the two cases lies in the factor appearing 
in the first parenthesis. We will call this dimensionless factor 


F = F ( N ’ ^) = 1 r + Tr 


(48) 


A 


The ratio of 6T^ for reference averaging to 6T^ for a conventional 
Dicke is given by 




*[l» ~ 


H 2 -l , N+l 

— + ~T 


(49) 


For example, if — - .1, going from a conventional Dicke to a reference 
P . • • • • : ; ; 

averaging system with N=6 will Increase 61^ by a factor of 120. Hence, 

while the term 5T^ may be ignored in a Dicke system, it must be dealt with 

in a reference averaging system. 

Figure 3 gives a plot of F /n, . 


Note that ~ affects F only for low values of N. For N > 3, the 


dependence on — • may be ignored, and 
P 


F « 


N 2 -l 


(50) 
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To give some feeling for the magnitude of 6T^, Equation (46) will 
be evaluated below for some typical system parameters and gain variabilities 
But first, it is useful to express 6a 2 in terms of more easily visualized 
parameters to describe gain variability. 

Sometimes gain stability may be specified by such a statement as 
"gain variations should be less than 1% at .01 Hz". This may be interpreted 
by a simple model of the form 


G(t) - G q + cos (tut + 0) 


where G q is constant, ~ < 1% , = .01 Hz and 0 is random. It can 

o 


be shown that if 


Nr < - 

p U> 


Nt Nt 

G (t) in (51) can be approximated in the interval ^ < t-t Q < 

by the first four terms of its power series expansion about t Q , 
with errors less than 4% of G^: 


G(t) = G(t o ) + G'(t o )(t-t o ) + — jp- (t-t/ + — (t-t Q ) 3 C 


Identifying the coefficients with those in Equation (30), we find that 

G"(t o ) 
a 2 2 

Now, as 0 varies randomly; G"(t o ) will vary over the range 
- U) 2 G < G"(t ) < tu 2 G 

ID ' O ID 


Therefore, the worst-case peak-to-peak variation in a is 



Substituting this into (46), and considering only the absolute value 
of 6T a , 


6 * Ui \ T p 2 ( N 2 -! + M+l M L +T 

6 T a- 4 g l 3 2 OV* > 


As an example, consider a reference averaging system with the following 


parameters : 


t r + t en = 3500 k 


T = 1 sec 

P 


and evaluate Equation (56) for gain variabilities described by 


S. = ..01 Hz 


G 

-rr = ,1%, .5% and 1% 

G 

The results are shown in Figure 4 as plots of versus N# According to 

condition (52), these curves are valid for values of N up to 


N < — « 32. 

TT 


For the case T. = T n , as in a null balanced radiometer, the -first 
term in Equation (29) vanishes, and the totaL temperature error due to 
gain variation is equal to ST^. This residual error is quite negligible 
when there is no reference averaging, since its value is ~.001°K, as 




indicated by the values of 6T^ for N = 1. But when reference averaging 
is used, this residual error grows rapidly with increasing N, as shown 
in Figure 4 . 


For cases where ^ T R , the two terms in Equation (29) must be 2 

combined to find the total gain variation error. The manner they are 
to be combined depends on how 6G and are correlated and on the 

sign and magnitude of - T R ; therefore, for these cases, it is p. 

advisable to treat each case separately, according to the specific 

characteristics oi the gain and scene temperature variations, rather 

than to generalize. 7 - 



TTTA 




5 . 


Detector Noise vs. Gain Variation Error Tradeoff and Estimates of 
Maximum Possible Improvement ~~~ — - 

Since gain variation error increases with N while detector noise 

error decreases with N, there must be some optimal value of N where the 

total antenna temperature error is minimum. To illustrate this tradeoff 

and optimisation, we will work out the case for a null-balanced reference 

averaging radiometer where T. - T n . 

A R 

Let 6T An and 5T Ag denote respectively the detector noise and gain 
variation errors, which are assumed to be uncorrelated. From Reference 1 
^An S^ ven by 


: 

I 


where 





SN 


= T + T 
K RN 


= T a + T 


RN 


is the system noise temperature. On the other hand, the gain variation 

error 6T Ag is gLven by Equation (56); so the total antenna temperature 

error 6T A consists of 
A 


6T a = / + < 6 V 2 


= _§S 

7 % 


414 

2+ § + ^% 
n 16G 2 


-1% 


N 2 -l , N+l T c \ 

— + --J % 
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Above, we have generalized the definition of the K- factor to include 

gain variation error. This generalized K- factor now depends not only 

T c 

on N, but also the dimensionless parameters — , Bt , and the gain 

variability dimensionless parameters m and -f 1 . 

V G 

As an example, we again choose the system parameters given in 
(57) and assume a gain fluctuation of ,VL at .01 Hz. For the MkSR 
water vapor 183 GHz channels, the typical IF bandwidths are 2.5 and 
1.25 GHz. With t « 1 sec, we have BT p = 2.5 x 10 and 1.25 x 10 , 
respecitvely. Assuming the aforementioned system parameters and gain 
variability, we evaluated the K-factor in Equation (60) as a function 
of N, for the too values of Bt . The results are shown in Figure 5. 


Br n » 2.5 x 10® 


Biwl.26 x IQ® 


FLUCTUATION 
AT 0.01 Hz 


ZERO GAIN FLUCTUATION 


1 2 3456789 10 20 SO 

NUMBER OF PIXELS, N 

FIGURE 5. GENERALIZED K FACTOR {COMBINED ERROR FROM 
DETECTOR NOISE AND GAIN FLUCTUATIONS) VERSUS N 









ATTACHMENT 1 TO APPENDIX B 


-> 


11 


The following integral will be evaluated approximately: 


/ 


p(t)t dt 


(A-l) 


where 


P(t) 


. ( 


\ 0 


alternately every At seconds 


as illustrated in Figure 2 of Appendix B, and 


(A-2) 


/ I 

?■- ‘J 


where 


t, - t = 2nAt 
d a 


At > 0 

n c positive integer > > 1 


CA-3) 


(A-4) 


The last condition means At < < 


fc b “ fc a 


* 


K 

I! ■ 


il 

l .. 


f $ 


Let t^, t^, t^, ... be the time at the centers of the intervals 
when p(t) = 1, as indicated in Figure 2. Consider the approximating 
sum that enters in the definition of the following integral: 

/ t 2 dt « t x 2 At + (t x +At) 2 At + t 2 2 At + (t 2 +At) 2 At + ... 

tt 


a 


[*.*< 

hi 


At + t 2 At + 


At + 2t 2 At + ... j At 


-] 

• • • 


+ (At)"(t b - t a ) 


(A-5) 


US 

j IV 

l 
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The sum in the second bracket is approximately equal to the following 
integral : 






At 

2 






(t b - t a )At 


Whereas, the sum in the first bracket is approximately equal to the 

integral (A-l), which we are seeking. Therefore, Equation (A-5) is 
approximately 


<A-6) 


\ 


i 



p(t)t dt + At 


't 

a 



f 


‘•-t 

i 


+ 


(At)-(t b -t a ). 

2 


(A-7) 


:: \ 


Substituting (A-6) into (A-7), we obtain 


; f. 
V-L 



"n 

. fj 

v a 


: r | 


b a 


2 ~~ - (At)-' 


tAt- 2 

ft a 


uh 
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Laser Beams and Resonators 

H. ICOGELNIIC and T. LI 


Abstract — ; This paper is a review of the theory of laser beams and 
resonators. It is meant to be tutorial in nature and useful in scope. No 
attempt is made to be exhaustive in the treatment. Rather, emphasis is 
placed on formulations and derivations which lead to basic understand- 
ing and on results which bear practical significance. 

1. Introduction 

T HE COHERENT radiation generated oy lasers or 
masers operating in the optical or infrared wave- 
length regions usually appears as a beam whose 
transverse extent is large compared to the wavelength. 
The resonant properties of such a beam in the resonator 
structure, its propagation characteristics in free space, and 
its interaction behavior wijh various optical elements and 
devices have been studied extensively in recent years. 
This paper is a review of the theory of laser beams and 
resonators. Emphasis is placed on formulations and 
derivations which lead to basic understanding and on 
results which are of practical value. 

Historically, the subject of laser resonators had its 
origin when Dicke [l ], Prokhorov [2], and Schawlow and 
Townes [3] independently proposed to use the Fabry- 
Perot interferometer as a laser resonator. The modes in 
such a structure, as determined by diffraction effects, 
were first calculated by Fox and Li [4]. Boyd and Gordon 
[5], and Boyd and Kogelnik [6] developed a theory for 
resonators with spherical mirrors and approximated the 
modes by wave beams. The concept of electromagnetic 
wave beams was also introduced by Goubau and Schwe- 
ring [7], who investigated the properties of sequences of 
lenses for the guided transmission of electromagnetic 
waves. Another treatment of wave beams was given by 
Pierce [8]. The behavior of Gaussian laser beams as they 
interact with various optical structures has been analyzed 
by Goubau [9], Kogelnik [10], [ll], and others. 

The present paper summarizes the various theories and 
is divided into three parts. The first part treats the passage 
of paraxial rays through optical structures and is based 
on geometrical optics. The second part is an analysis of 
laser beams and resonators, taking into account the wave 
nature of the beams but ignoring diffraction effects due 
to the finite size of the apertures. The third part treats the 
resonator modes, taking into account aperture diffrac- 
tion effects. Whenever applicable, useful results are pre- 
sented in the forms of formulas, tables, charts, and 
graphs. 
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2. Paraxial Ray Analysis 

A study of the passage of paraxial rays through optical 
resonators, transmission lines, and similar structures can 
reveal many important properties of these systems. One 
such “geometrical” property is the stability of the struc- 
ture [6], another is the loss of unstable resonators [12]. 
The propagation of paraxial rays through various optical 
structures can be described by ray transfer matrices. 
Knowledge of these matrices is particularly useful as they 
also describe the propagation of Gaussian beams through 
these structures; this will be discussed in Section 3. The 
present section describes briefly some ray concepts which 
are useful in understanding laser beams and resonators, 
and lists the ray matrices of several optical systems of 
interest. A more detailed treatment of ray propagation 
can be found in textbooks [13] and in the literature on 
laser resonators [14]. 
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Fig. 1 . Reference planes of an optical system. 
A typical ray path is indicated. 


2.1 Ray Transfer Matrix 

A paraxial ray in a given cross section (z= const) of an 
optical system is characterized by its distance x from the 
optic ( z ) axis and by its angle or slope x! with respect to 
that axis. A typical ray path through an optical structure 
is shown in Fig. 1. The slope x' of paraxial rays is assumed 
to be small. The ray path through a given structure de- 
pends on the optical properties of the structure and on the 
input conditions, i.e., the position xi and the slope x t ' of 
the ray in the input plane of the system. For paraxial rays 
the corresponding output quantities x« and x 2 ' are linearly 
dependent on the input quantities. This is conveniently 
written in the matrix form 


Xi ; 


A 

B 

Xi 

x/ 


C 

D 

•Tl' 


C-3 


TABLE I 

Ray Transfer Matrices of Six Elementary Optical Structures 



where the slopes are measured positive as indicated in the 
figure. The ABCD matrix is called the ray transfer matrix. 
Its determinant is generally unity 

AD - BC = 1. (2) 

The matrix elements are related to the focal length / of 
the system and to the location of the principal planes by 



hi — 


D- 1 
C 


(3) 



where h\ and h > are the distances of the principal planes 
from the input and output planes as shown in Fig. 1. 

In Table I there are listed the ray transfer matrices of 
six elementary optical structures. The matrix of No. 1 
describes the ray transfer over a distance d. No. 2 de- 
scribes the transfer of rays through a thin lens of focal 
length f Here the input and output planes are immediately 
to the left and right of the lens. No. 3 is a combination 
of the first two. It governs rays passing first over a dis- 
tance d and then through a thin lens. If the sequence is 
reversed the diagonal elements are interchanged. The 
matrix of No. 4 describes the rays passing through two 
structures of the No. 3 type. It is obtained by matrix 
multiplication. The ray transfer matrix for a lenslike 
medium of length d is given in No. 5. In this medium the 
refractive index varies quadratically with the distance r 
from the optic axis. 

n = 71 o — hn«r"-. (4) 


An index variation of this kind can occur in laser crystals 
and in gas lenses. The matrix of a dielectric material of 
index n and length d is given in No. 6. The matrix is 
referred to the surrounding medium of index 1 and is 
computed by means of Snell’s law. Comparison with No. 
1 shows that for paraxial rays the effective distance is 
shortened by the optically denser material, while, as is 
well known, the “optical distance" is lengthened. 

2.2 Periodic Sequences 


Light rays that bounce back and forth between the 
spherical mirrors of a laser resonator experience a periodic 
focusing action. The effect on the rays is the same as in a 
periodic sequence of lenses [15] which can be used as an 
optical transmission line. A periodic sequence of identical 
optical systems is schematically indicated in Fig. 2. A 
single element of the sequence is characterized by its 
ABCD matrix. The ray transfer through n consecutive 
elements of the sequence is described by the nth power 
of this matrix. This can be evaluated by means of Sylves- 
ter’s theorem 


A B » 
C D 


1 

sin 0 



A sin 7i0 — sin (n — 1)0 B sin ?iC-) 

C sin 7i0 D sin 7i0 — sin(n — 1)0 
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where 


cos C-) = $(,t + D). (0) 

Periodic sequences can be classified as either stable or 
unstable. Sequences are stable when the trace D) 
obeys the inequality 

-1 < + D) < 1. (7) 

Inspection of (5) shows that rays passing through a stable 
sequence are periodically refocused. For unstable sys- 
tems, the trigonometric functions in that equation be- 
come hyperbolic functions, which indicates that the rays 
become more and more dispersed (he further they pass 
thicugh the sequence. 
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Fig. 3. Spherical-mirror resonator and the 
equivalent sequence of lenses. 



Fig. 2. Periodic sequence of identical systems, 
each characterized by its ADCD matrix. 


2.3 Stability of Laser Resonators 

A lase resonator with spherical mirrors of unequal 
curvatur- is a typical example of a periodic sequence that 
can be either stable or unstable [6], In Fig. 3 such a 
resonator is shown together with its dual, which is a 
sequence of lenses. The ray paths through the two struc- 
tures are the same, except that the ray pattern is folded in 
the resonator and unfolded in the lens sequence. The focal 
lengths J\ and /- of the lenses are the same as the focal 
lengths of the mirrors, i.e., they are determined by the 
radii of curvature and R» of the mirrors (f=R h ’2 t 
f<=R< 2). The lens spacings are the same as the mirror 
spacing cl. One can choose, as an element of the peri- 
odic sequence, a spacing followed by one lens plus another 
spacing followed by the second lens. The A BCD matrix 
of such an element is given in No. 4 of Table I. From this 
one can obtain the trace, and write the stability condition 
(7) in the form 


0 < ('“f)( 1 "S <i - (8) 


To show graphically which type of resonator is stable 
and which is unstable, it is useful to plot a stability dia- 
gram on which each resonator type is represented by a 
point. This is shown in Fig. 4 where the parameters d/R x 
and d R» are drawn as the coordinate axes; unstable 
systems are represented by points in the shaded areas. 
Various resonator types, as characterized by the relative 
positions of the centers of curvature of the minors, are 
indicated in the appropriate regions of the diagram. Also 
entered as alternate coordinate axes are the parameters g, 
and gi which play an important role in the diffraction 
theory of resonators (see Section 4). 
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Fig. 4. Stability diagram. Unstable resonator 
systems lie in shaded regions. 


3. Wave Analysis of Beams and Resonators 

In this section the wave nature of laser beams is taken 
into account, but diffraction effects due to the finite size 
of apertures are neglected. The latter will be discussed in 
Section 4. The results derived here are applicable to 
optical systems with "large apertures," i.e., with apertures 
that intercept only a negligible portion of the beam power. 
A theory of light beams or "beam waves" of this kind was 
first given by Boyd and Gordon [5] and by Goubau and 
Schwering [7]. The present discussion follows an analysis 
given in [II]. 

3.1 Approximate Solution of the Wave Equation 

Laser beams are similar in many respects to plane 
waves; however, their intensity distributions are not uni- 
form, but are concentrated near the axis of propagation 
and their phase fronts are slightly curved. A field com- 
ponent or potential u of the coherent light satisfies the 
scalar wave equation 

+ k'U' — 0 ( 0 ) 

where k= 2ir/X is the propagation constant in the medium. 


C-5 



For light traveling in the z direction one writes 

u = P(x, y, z ) exp(-jA'z) (10) 

where ^ is a slowly varying complex function which 
represents the differences between a laser beam and a 
plane wave, namely: a nonuniform intensity distribu- 
tion, expansion of the beam with distance of propagation, 
curvature of the phase front, and other differences dis- 
cussed below. By inserting (10) into (9) one obtains 


dhb d*tp dtp 

“+-“-* 2 jk~ 
dx'- d>f dz 


= o 


(U) 


where it has been assumed that tp varies so slowly with z 
that its second derivative Shp/dz- can be neglected. 

The differential equation (11) for tp has a form similar 
to the time dependent Schrddinger equation. It is easy to 
see that 

*-"p (i2> 

is a solution of (11), where 

i* = ** + »*. (13) 


The parameter P(z) represents a complex phase shift which 
is associated with the propagation of the light beam, and 
q(z) is a complex beam parameter which describes the 
Gaussian variation in beam intensity with the distance r 
from the optic axis, as well as the curvature of the phase 
front which is spherical near the axis. After insertion of 
(12) into (11) and comparing terms of equal powers in r 
one obtains the relations 


II 

(14) 

-O 

II 

1 

1 

(15) 


<1 


where the prime indicates differentiation with respect to z. 
The integration of (14) yields 

q s ~qi + z (16) 


which relates the beam parameter q 2 in one plane (output 
plane) to the parameter <71 in a second plane (input plane) 
separated from the first by a distance z. 



Fig. 5. Amplitude distribution of the fundamental beam. 


When (17) is inserted in (12) the physicat meaning of these 
two parameters becomes clear. One sees that R(z) is the 
radius of curvature of the wavefront that intersects the 
axis at z, and u(z) is a measure of the decrease of the 
field amplitude E with the distance from the axis. This 
decrease is Gaussian in form, as indicated in Fig. 5, and 
w is the distance at which the amplitude is 1/e times that 
on the axis. Note that the intensity distribution is Gaus- 
sian in every beam cross section, and that the width of 
that Gaussian intensity profile changes along the axis. 
The parameter w is often called the beam radius or “spot 
size,” and 2 iv, the beam diameter. 

The Gaussian beam contracts to a minimum diameter 
2u’ 0 at the beam waist where the phase front is plane. If 
one measures z from this waist, the expansion laws for 
the beam assume a simple form. The complex beam 
parameter at the waist is purely imaginary 


7TUV 
<70 = j — — 


(18) 


and a distance z away from the waist the parameter is 


irttfo 2 

q - Qq + 2 = j — h 2. 

A 


(19) 


After combining (19) and (17) one equates the real and 
imaginary parts to obtain 


and 


w-(z) = Wa- 


ft (z) = z 



(20) 


( 21 ) 


3.2 Propagation Laws for the Fundamental Mode 

A coherent light beam with a Gaussic. i intensity pro- 
file as obtained above is not the only solution of (11), 
but is perhaps the most important one. This beam is often 
called the “fundamental mode” as compared to the higher 
order modes to be discussed later. Because of its impor- 
tance it is discussed here in greater detail. 

For convenience one introduces two real beam param- 
eters R and w related to the complex parameter q by 


Figure 6 shows the expansion of the beam according to 
(20). The beam contour w{z) is a hyperbola with asymp- 
totes inclined to the axis at an angle 


\ 

*100 


( 22 ) 


This is the far-field diffraction angle of the fundamental 
mode. 

Dividing (21) by (20), one obtains the useful relation 



\z ~u>- 
1TW 0 - XR 


(23) 
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Fig. 6. Contour of a Gaussian beam. 


which can be used to express wo and s in terms of w and R : 

w " s=w! /[ i+ S) ! ] (24) 

z= V[ i+ OT i2s) 

To calculate the complex phase shift a distance z away 
from the waist, one inserts (19) into (15) to get 

P f = - — = - (26) 

Q s + iOrWo'VA) 

Integration of (26) yields the result 

;P(z) = ln[l -./(As/irWo 5 )] 

= InVl 4- (Az/ 7mto : ) 2 — j arc tan (\z/ irw 0 2 ) . (27) 


ilar properties, and they are discussed in this section. 
These solutions form a complete and orthogonal set of 
functions and are called the “modes of propagation.” 
Every arbitrary distribution of monochromatic light can 
be expanded in terms of these modes. Because of space 
limitations the derivation of these modes can only be 
sketched here. 

a) Modes in Cartesian Coordinates: For a system with* 
a rectangular (.v, y, z) geometry one can try a solution for 
(11) cf trie form 




p + ~ (a - 8 + IP) 

2q 


Ji 


(HO) 


where g is a function of .v and z, and It is a function of y 
and z. For real g and h this postulates mode beams whose 
intensity patterns scale according to the width 2w(c) of a 
Gaussian beam. After inserting this trial solution into 
(11) one arrives at differential equations for .if and //of the 
form 


d-Ii m 

dx- 


dii M 

2.r f- 2 m// w = 0. 

dx 


( 31 ) 


This is the differential equation for the Hermite poly- 
nomial H m { x) of order in. Equation (1 1) is satisfied if 

B • A = II,. ( s/2 II. (\/2 (32) 


The real part of P represents a phase shift difference <£ be- 
tween the Gaussian beam and an ideal plane wave, while 
the imaginary part produces an amplitude factor wa/w 
which gives the expected intensity decrease on the axis due 
to the expansion of the beam. With these results for the 
fundamental Gaussian beam, (10) can be written in the 
form 



where 

<I> = arc tan(X2/7rt/’o 2 ). (29) 

It will be seen in Seotion 3.5 that Gaussian beams of this 
kind are produced by many lasers that oscillate in the 
fundamental mode. 

3.3 Higher Order Modes 

In the. preceding section only one solution of (II) was 
discussed, i.e., a light beam with the property that its 
intensity profile in every beam cross sectiqn is given by 
the same function, namely, a Gaussian. The width of this 
Gaussian distribution changes as the beam propagates 
along its axis. There are other solutions of (11) with sim- 
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where m and n are the (transverse) mode numbers. Note 
that the same pattern scaling parameter «;(z) applies to 
modes of all orders. 

Some Hermite polynomials of low order are 
ffoCr) = 1 
//,(*) = -r 
H*(x) — 4x- — 2 

Hz(x) = 8r 3 - V2x. (63) 

.Expression (28) can be used as a mathematical descrip- 
tion of higher order light beams, if one inserts the product 
gdi as a factor on the right-hand side. The intensify pat- 
tern in a cross section of a higher order beam is, thus, de- 
scribed by the product of Hermite and Gaussian functions. 
Photographs of such mode patterns are shown in Fig. 7. 
They were produced as modes of oscillation in a gas laser 
oscillator [16]. Note that the number of zeros in a mode 
pattern is equal to the corresponding mode number, and 
that the area occupied by a mode increases with the mode 
number. 

The parameter R( ?) in (28) is the same for all modes, 
implying that the phase-front curvature is the same and 
changes in the same way fi r modes of all orders. The 
phase shift <I>, however, is a function of the mode numbers. 
One obtains 

‘Km, 2 ) — (m + n + l) arc tan(A2/7rta 0 5 ). (34) 
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Fig. 7. Mode patterns of a gas laser oscil- 
lator (rectangular symmetry). 


This means that the phase velocity increases with increas- 
ing mode number. In resonators this leads to differences 
in the resonant frequencies of the various modes of oscil- 
lation. 

h) Modes in Cylindrical Coordinates: For a system with 
a cylindrical ( r , 0 , z) geometry one uses a trial solution 
for (l 1) of the form 

W p+ 5’ , + 4 )}- (35) 

After some calculation one finds 

„ = (v2 0-v(2-3 (30) 

where L.J is a generalized Laguerre polynomial, and p 
and / are the radial and angular mode numbers. L p '(x ) 
obeys the differential equation 

d*L p ' (1L„ l 

x 4 (/ 4- 1 - x) -f- + pL p ' = 0. (37) 

dx 2 ax 

Some polynomials of low order are 

WW = 1 

Ldix) = / 4- 1 — ** 

L 2 ‘(x) = J(/+ 1 HI 4- 2) - (/ 4- 2)x 4- k 2 - (38*, 


As in the case of beams with a rectangular geometry, the 
beam parameters uiz) and R(z) are the same for all cylin- 
drical modes. The phase shift is, again, dependent on the 
mode numbers and is given by 

4>(p, l;z) = (2p 4- / 4- 1) arc tiui(Xz/irt 0 o*). (39) 


3.4 Beam Transformation by a Lens 

A lens can be used to focus a laser beam to a small spot, 
or to produce a beam of suitable diameter and phase- 
front curvature for injection into a given optical structure. 
An ideal lens leaves the transverse field distribution of a 
beam mode unchanged, i.e., an incoming fundamental 
Gaussian beam will emerge from the lens as a funda- 
mental beam, and a higher order mode remains a mode 
of the same order after passing through the lens. However, 
a lens does change the beam parameters R(z ) and uiz). 
As these two parameters are the same for modes of all 
orders, the following discussion is valid for all orders; 
the relationship between the parameters of an incoming 
beam (labeled here with the index 1) and the parameters 
of the corresponding outgoing beam (index 2) is studied in 
detail. 

An ideal thin lens of focal length / transforms an incom- 
ing spherical wave with a radius R 1 immediately to the 
left of the lens into a spherical wave with the radius R 2 
immediately to the right of it, where 


11 1 

Ri R\ f 


(40) 


Figure 8 illustrates this situation. The radius of curvature 
is taken to be positive if the wavefront is convex as 
viewed from z= sc. The lens transforms the phase fronts 
of laser beams in eactly the same way as those of spherical 
waves. As the diameter of a beam is the same immediately 
to the left and to th-' right of a thin lens, the (/-parameters 
of the incoming and outgoing beams are related by 


1 _ 1 1 

92 q\ f 


(41) 


where the q’ s are measured at the lens. If q 1 and r/ 2 are 
measured at distances d y and d 2 from the lens as indicated 
in Fig. 9, the relation between them becomes 

(1 — dt/f)qi 4- (rfi 4- dt — didt/f) 

(lo — ——————————— ————— • 

-(<?.//) + (! ~dx/S) 


This formula is derived using (16) and (41). 

More complicated optical structures, such as gas lenses, 
combinations of lenses, or thick lenses, can be thought of 
as composed of a series of thin lenses at various spacings. 
Repeated application of (16) and (41) is, therefore, suffi- 
cient to calculate the effect of complicated structures on 
the propagation of laser beams. If the A BCD matrix for 
the transfer of paraxial rays through the structure is 
known, the q parameter of the output beam can be cal- 
culated from 
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Fig. 8. Transformation of wavefronts by a thin lens. 
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Fig. 9. Distances and parameters for a 
beam transformed by a thin lens. 


trip. If the complex beam parameter is given by q u im- 
mediately to the right of a particular lens, the beam 
parameter q 2 , immediately to the right of the next lens, 
can be calculated by means of (16) and (41) as 


92 ?x + d f 


Self-consistency requires that qi=q 2 =q, which leads to 
a quadratic equation for the beam parameter q at the lenses 
(or at the mirrors of the resonator): 


1 1 1 

1 1 = 0 . 

q- fq fd 


The roots of this equation are 


1 1 - / 11 

F" -y (+)j Vw~7r 


(45) 


(46) 


Aqi + B 
C?T+ D 


(43) 


This is a generalized form of (42) and has been called the 
A BCD law [10]. The matrices of several optical structures 
are given in Section II. The A BCD law follows from the 
analogy between the laws for laser beams and the laws 
obeyed by the spherical waves in geometrical optics. The 
radius of the spherical waves R obeys laws of the same 
form as (16) and (41) for the complex beam parameter q. 
A more detailed discussion of this analogy is given in [l 1 ]. 


3.5 Laser Resonators ( Infinite Aperture) 

The most commonly used laser resonators are com- 
posed of two spherical (or flat) mirrors facing each other. 
The stability of such “open” resonators has been discussed 
in Section 2 in terms of paraxial rays. To study the modes 
of laser resonators one has to take account of their wave 
nature, and this is done here by studying wave beams of 
the kind discussed above as they propagate back and forth 
between the mirrors. As aperture diffraction effects are 
neglected throughout this section, the present discussion 
applies only to stable resonators with mirror apertures 
that are large compared to the spot size of the beams. 

A mode of a resonator is defined as a self-consistent 
field configuragion. If a mode can be represented by a 
wave beam propagating back and forth between the 
mirrors, the beam parameters must be the same after one 
complete return trip of the beam. This condition is used 
to calculate the mode parameters. As the beam that repre- 
sents a mode travels in both directions between the mirrors 
it forms the axial standing-wave pattern that is expected 
for a resonator mode. 

A laser resonator with mirrors of equal curvature is 
shown in Fig. 10 together with the equivalent unfolded 
system, a sequence of lenses. For this symm .trical struc- 
ture it is sufficient to postulate self-consistency for one 
transit of the resonator (which is equivalent to one full 
period of the lens sequence), instead of a complete return 


where only the root that yields a real beamwidth is used. 
(Note that one gets a real beamwidth for stable resonators 
only.) 

From (46) one obtains immediately the real beam 
parameters defined in (17). One sees that R is equal to the 
radius of curvature of the mirrors, which means that the 
mirror surfaces are coincident with the phase fronts of 
the resonator modes. The width 2 w of the fundamental 
mode is given by 

To calculate the beam radius w 0 in the center of the reso- 
nator where the phase front is plane, one uses (23) with 
z=d/ 2 and gets 

wo- = — Vd(2R - d). (48) 

27T 

The beam parameters R and w describe the modes of 
all orders. But the phase velocities are different for the 
different orders, so that the resonant conditions depend on 
the mode numbers. Resonance occurs when the phase 
shift from one mirror to the other is a multiple of ir. 
Using (28) and (34) this condition can be written as 

led — 2 (m A- n + 1) arc ta,n(\d/2irw 0 2 ) = *■(? 4- 1) (49) 

where q is the number of nodes of the axial standing-wave 
pattern (the number of half wavelengths is q+ 1), 1 and m 
and n are the rectangular mode numbers defined in Sec- 
tion 3.3. For the modes of circular geometry one obtains 
a similar condition where (2p+/+l) replaces (m+n+1). 

The fundamental beat frequency v 0 , i.e., the frequency 
spacing between successive longitudinal resonances, is 
given by 

vo = c/2d (50) 

1 This q is not to be confused with the complex beam parameter. 
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Fig. 10. Symmetrical laser resonator and the equivalent sequence 
of lenses. The beam parameters, <71 and are indicated. 
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Fig. II. Mode parameters or interest for a resonntor with 
mirrors of unequal curvature. 


where c is the velocity of light. After some algebraic 
manipulations one obtains from (49) the following for- 
mula for the resonant frequency v of a mode 

»/vb =(<H- 1)4 — (m+n+1) arc cos(l— d/R) (51) 

7T 

For the special case of the confocal resonator \d - R ~ b), 
the above relations become 

w- = Xft/V, Wo 2 = \b/2ir] 

»!h “ (q + 1) + -f n -}- I). (52) 


The parameter b is known as the confocal parameter. 

Resonators with mirrors of unequal curvature can be 
treated in a similar manner. The geometry of such a 
resonator where the radii of curvature of the mirrors are 
R ! and R 2 is shown in Fig. 1 1. The diameters of the beam 
at the mirrors of a stable resonator, 2ir, and 2w», are 
given by 


H'l 


4 


(\R \! 7r) 2 


Ri - d 

~r^7i 


d 

ill + Rn-d 


1C ; 1 


- (X/ij/ir)* 


Ri -d d 

Ri — d Ri + Ri — d 


(53) 


The diameter of the beam waist 2w 0 , which is formed 
either inside or outside the resonator, is given by 


,^\ 2 d(R y - d)(R, - d)(Ri + R*~ d) 


Wo = 


(Ry + R,~ 3rf)* 


(54) 


The distances ty and /« between the waist and the mirrors, 
measured positive as show '. in the figure, are 


d(Ro - d ) 

Ri + Ri - 2d 
d(Rt - d) 

Ri + R > - 2d ’ 


( 55 ) 


The resonant condition is 
i'/vo = (q + 1) 4 (w + n + 1) 


x 

arc eos\/(l — d/Ri)(l — d/R ■>) 


(50) 


where the square root should be given the sign of ( 1 — df Ry), 
which is equal to the sign of (1 —d/R») for a stable resona- 
tor. 

There are more complicated resonator structures than 
the ones discussed above. In particular, one can insert a 
lens or several lenses between the mirrors. But in every 
case, the unfolded resonator is equivalent to a periodic 
sequence of identical optical systems as shown in Fig. 2. 
The elements of the A BCD matrix of this system can be 
used to calculate the mode parameters of the resonator. 
One uses the ABCD law (43) and postulates self-con- 
sistency by putting (ji=q>=q. The roots of the resulting 
quadratic equation are 

7 - ( +> is m 

which yields, for the corresponding beam radius w, 

w s = (2\b/tv) / ( 58 ) 


3.6 Mode Matching 

It was shown in the preceding section that the modes of 
laser resonators can be characterized by light beams with 
certain properties and parameters which are defined by 
the resonator geometry. These beams are often injected 
into other optical structures with different sets of beam 
parameters. These optical structures can assume various 
physical forms, such as resonators used in scanning 
Fabry-Perot interferometers or regenerative amplifiers, 
sequences of dielectric or gas lenses used as optical trans- 
mission lines, or crystals of nonlinear dielectric material 
employed in parametric optics experiments. To match 
the modes of one structure to those of another one must 
transform a given Gaussian beam (or higher order mode) 
into another beam with prescribed properties. This trans- 
formation is usually accomplished with a thin lens, but 
otiier more complex optical systems can be used. Although 
the present discussion is devoted to the simple case of the 
thin lens, it is also applicable to more complex systems, 
provided one measures the distances from the principal 
planes and uses the combined focal length /of the more 
complex system. 

The location of the waists of the two beams to be 
transformed into each other and the beam diameters at 
the waists are usually known or can be computed. To 
match the beams one has to choose a lens of a focal length 
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/ that is larger than a characteristic length f 0 defined by 
the two beams, and one has to adjust the distances be- 
tween the lens and the two beam waists according to rules 
derived below. 

In Fig. 9 the two beam waists are assumed to be 
located at distances d\ and d« from the lens. The complex 
beam parameters at the waists are purely imaginary; they 
are ; 


0i = jvw r/X, 02 = jwWa*/\ (59) 


where 2w x and 2 w» are the diameters of the two beams at 
their waists. If one inserts these expressions for q x and qi 
into (42) and equates the imaginary parts, one obtains 


di—f _ wr 
di *-* f w 2 2 


(60) 


Equating the real parts results in 

Wi = P~f 0 2 


where 


(61) 


fo = TTWiWi/X. (62) 

Note that the characteristic length f Q is defined by the waist 
diameters of the beams to be matched. Except for the 
term fa 2 , which goes to zero for infinitely small wave- 
lengths, (61) resembles Newton’s imaging formula of 
geometrical optics. 

Any lens with a focal length f>f Q can be used to per- 
form the nut<chfo>; transformation. Once/is chosen, the 
distances d x and d« have to be adjusted to satisfy the 
matching formulas [10] 


di~ f ± — VP - fo 2 , 


d 2 =/±— VP-fo 2 . (63) 

Wi 

These relations are derived by combining (60) and (61). 
In (63) one can choose either both plus signs or both 
minus signs for matching. 

It is often useful to introduce the confocal parameters 
by and b 2 into the matching formulas. They are defined 
by the waist diameters of the two systems to be matched 

bi = 2mui 2 /\ t bi = 2tw 2 2 /\. (64) 

Using these parameters one gets for the characteristic 
length/o 

fo 2 = 5&1&2, (65) 

and for the matching distances 

di=f± *&i VWh 2 ) - 1, 

d 2 =f± hh V(P7W)-=1. ( 66 ) 

Note that in this form of the matching formulas, the 
wavelength does not appear explicitly. 

Table II lists, for quick reference, formulas for the two 
important parameters of beams that emerge from various 


TABLE II 

Formulas for the Confocal Parameter and the Location of 
Beam Waist for Various Optical Structures 



optical structures commonly encountered. Th, are the 
confocal parameter b and the distance t whicu gives the 
waist location of the emerging beam. System No. 1 is a 
resonator formed by a flat mirror and a spherical mirror 
of radius R. System No. 2 is a resonator formed by two 
equal spherical mirror a . System No. 3 is a resonator 
formed by mirrors of unequal curvature. System No. 4 
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Fig. 12. The confocal parameter b : as a rune- 
lion of the lens-waist spacing </,. 



Fig. 1 3. The waist spacing as a function of 
l lie lens-waist spacing d,. 


is, again, a resonator formed by two equal spherical mir- 
rors, but with the reflecting surfaces deposited on plano- 
concave optical plates ol index n. These plates act as 
negative lenses and change the characteristics of the 
emerging beam. This lens effect is assumed not present in 
Systems Nos. 2 and 3. System No. 5 is a sequence of 
thin lenses of equal focal lengths /. System No. 6 is a 
system of two irises with equal apertures spaced at a 
distance d. Shown are the parameters of a beam that 
will pass through both irises with the least possible beam 
diameter. This is a beam which is “confocal** over the 
distance d. This beam will also pass through a tube of 
length d with the optimum clearance. (The tube is also 
indicated in the figure.) A similar situation is shown in 
System No. 7, which corresponds to a beam that is 
confocal over the length d of optical material of index n. 
System No. 8 is a spherical mirror resonator filled with 
material of index n. or an optical material with curved 
end surfaces where the beam passing through it is as- 
sumed to have phase fronts that coincide with these sur- 
faces. 

'■■hen one designs a matching system, it is useful to 
know the accuracy required of the distance adjustments. 
The discussion below indicates how the parameters b~> and 
r/'j change when b, and f are fixed and the lens spacing d, 
to the waist of the input beam is varied. Equations (60) 
and (61) can be solved for with the result [9] 

. ,■ . f 

<>2 J ~ - . ((> 7 ) 

U -ih VY 


This means that the parameter b 2 of the beam emerging 
from the lens changes with d x according to a Lorentzian 
functional form as shown in Fig. 12. The Lorentzian is 
centered at d x =f and has a width of b u The maximum 
value of 6o is 4 f*/b t . 

If one inserts (67) into (60) one gets 



1 - rf,// _ 

- <U/fY +(b h nty 


(08) 


which shows the change of d 2 with d x . The change is 
reminiscent of a dispersion curve associated with a 
Lorentzian as shown in Fig. 13. The extrema of this curve 
occur at the halfpower points of the Lorentzian. The slope 
of the curve at t/i=/is (2 f/b )-. The dashed curves in the 
figure correspond to the geometrical optics imaging re- 
lation between d u d-i t and/ [20]. 


' 7 Circle Diagrams 

The propagation of Gaussian laser beams can be repre- 
bL 'ted graphically on a circle diagram. On such a diagram 
one can follow a beam as it propagates in free space or 
passes through lenses, thereby affording a graphic solu- 
tion of the mode matching problem. The circle diagrams 
for beams are similar to the impedance charts, such as the 
Smith chart. In fact there is a close analogy between 
transmission-line and laser-beam problems, and there are 
analog electric networks for every optical system [17], 
The first circle diagram for beams was proposed by 
Collins [18]. A dual chart was discussed in [19]. The 
basis for the derivation of these charts are the beam prop- 
agation laws discussed in Section 3.2. One combines 
(17) and (19) and eliminates q to obtain 



Th. relation contains the four quantities u\ R, w 0 , and r 
which were used to describe the piopay tion of Gaussian 
beams in Section 3.2. Each pair of the- lantities can be 
expressed in complex variables \V anc 


TTlt'o' 

z = — - j* = W 2 - jz. (70) 


where b is the confocal parameter of the beam. For these 
variables (69) defines a conformal transformation 

»’ = 1 /*. ( 71 ) 

The two dual circle diagrams are plotted in the complex 
planes of W and Z, respectively. The IF-plane diagram 
[IS) is shown in Fig. 14 where the variables X W- and 
1 R are plotted as axes. In this plane the lines of constant 
b 2 = jthv' X and the lines of constant z of the Z plane 
appear as circles through the origin. A beam is represented 
by a circle of constant h, and the beam parameters w and 
R at a distance z from the beam waist can be easily read 
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Fig. 14, Geometry lor the W'-plane circle diagram. 
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Fig. 15. The Gnussian beam chart. Both IP-plane and Z-pIane 
circle diagram are combined into one. 


from the diagram. When the beam passes through a lens 
the phase front is changed according to (40) and a new 
beam is formed, which implies that the incoming and 
outgoing beams are connected in the diagram by a 
vertical line of length l /f The angle shown in the figure 
is equal to the phase shift experienced by the beam as 
given by (29): this is easily shown using (23). 

The dual diagram [19] is plotted in the Z plane. The 
sets of circles in both diagrams have the same form, and 
only the labeling of the axes and circles is different. In 
Fig. 15 both diagrams are unified in one chart. The 
labels in parentheses correspond to the Z-plane diagram, 
and b is a normalizing parameter which can be arbitrarily 
chosen for convenience. 

One can plot various other circle diagrams which are 
related to the above by conformal transformations. One 
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such transformation makes it possible to use the Smith 
chart for determining complex mismatch coefficients for 
Gaussian beams [20]. Other circle diagrams include those 
for optical resonators [21] which allow the graphic deter- 
mination of certain parameters of the resonator modes. 

4. Laser Resonators (Finite Aperture) 

4. 1 General Mathematical Formulation 

In this section aperture diffraction effects due to the 
finite size of the mirrors are taken into account; these 
effects were neglected in the preceding sections. There, 
it was mentioned that resonators used in laser oscillators 
usually take the form of an open structure consisting of a 
pair of mirrors facing each other. Such a structure with 
finite mirror apertures is intrinsically lossy and, unless 
energy is supplied to it continuously, the electromagnetic 
field in it will decay. In this case a mode of the resonator 
is a slowly decaying field configuration whose relative 
distribution does not change with time [4]. In a laser 
oscillator the active medium supplies enough energy to 
overcome the losses so that a steady-state field can exist. 
However, because of nonlinear gain saturation the me- 
dium will exhibit less gain in those regions where the 
field is high than in those where the field is low, and so the 
oscillating modes of an active ^senator are expected to 
be somewhat different from the decaying modes of the 
passive resonator. The problem of an active resonator 
filled with a saturable-gain medium has been solved re- 
cently [22], [23], and the computed results show that if 
the gain is not too large the resonator modes are essen- 
tially unperturbed by saturation effects. This is fortunate 
as the results which have been obtained for the passive 
resonator can also be used to describe the active modes of 
laser oscillators. 

The problem of the open resonator is a difficult one 
and a rigorous solution is yet to be found. However, if 
certain simplifying assumptions are made, the problem 
becomes tractable and physically meaningful results can 
be obtained. The simplifying assumptions involve essen- 
tially the quasi-optic nature of the problem ; specifically, 
they are 1) that the dimensions of the resonator are large 
compared to the wavelength and 2) that the field in the 
resonator is substantially transverse electromagnetic 
(TEM). So long as those assumptions are valid, the 
Fresnel-Kirchhoff formulation of Huygens’ principle can 
be invoked to obtain a pair of integral equations which 
relate the fields of the two opposing mirrors. Further- 
more, if the mirror separation is large compared to mirror 
dimensions and if the mirrors are only slightly curved, 
the two orthogonal Cartesian components of the vector 
field are essentially uncoupled, so that separate scalar 
equations can be written for each component. The solu- 
tions of these scalar equations yield resonator modes 
which are uniformly polarized in one direction. Other 
polarization configurations cu„ be constructed from the 
uniformly polarized modes by linear superposition. 
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Fig. 16. Geometry of a spherical-mirror resonator with finite 
mirror apertures and the equivalent sequence of lenses set in 
opaque absorbing screens. 


In deriving the integral equations, it is assumed that a 
traveling TEM wave is reflected back and forth between 
the mirrors. The resonator is thus analogous to a trans- 
mission medium consisting of apertures or lenses set in 
opaque absorbing screens (see Fig. 16). The fields at the 
two mirrors are related by the equations [24} 

y “>£<»(«,) - f s i )E^(s 2 )dS 2 

= f K<»(s a> Si)A T <»(s,)d5i (72) 

J Si 

where the integrations are taken over the mirror surfaces 
S 2 and Si, respectively. In the above equations the sub- 
scripts and superscripts one and two denote mirrors one 
and two; Si and s 2 are symbolic notations for transverse 
coordinates on the mirror surface, e.g., .\t=(.x',, y{) and 
s 2 =(x 2 , Ti) or *i=0'i> anci s«=(r 2> <£ 2 ); £ u> and £ m 
are the relative field distribution functions over the mir- 
rors; 7 (1> and y m give the attenuation and phase shift 
suffered by the wave in transit from one mirror to the 
other; the kernels K w and K (2> are functions of the dis- 
tance between .v, and s-< and, therefore, depend on the 
mirror geometry; they are equal [/C n) (.v 2t .s'i)- K (2) (s ■, s 2 ) ] 
but, in general, are not symmetric [/£ ( "(.v 2 , .Vi);^/f (1) (.r,, .v.), 

K*(s u y»)|. 

The integral equations given by (72) express the field 
at each mirror in terms of the reflected field at the other; 
that is, they are single-transit equations. By substituting 
one into the other, one obtains the double-tiansit or 
round-trip equations, which state that the field at each 
mirror must reproduce itself after a round trip. Since the 
kernel for each of the double-transit equations is sym- 
metric [24], it follows [25] that the field distribution 
functions corresponding to the different mode orders are 
orthogonal over their respective mirror surfaces; that is 


/. 

/. 


B m <»(8i)BS»(8ddS l = 0, 


EJ'-KsJEnt'-KsJdSi = 0, m ^ » (78) 


where m and n denote different mode orders. It is to be 
noted that the orthogonality relation is non-Hermitian 
and is the one that is generally applicable to lossy sys- 
tems. 

4.2 Existence of Solutions 

The question of the existence of solutions to the 
resonator integral equations has been the subject of 
investigation by several authors [26]— [28 ]. They have 
given rigorous proofs of the existence of eigenvalues and 
eigenfunctions for kernels which belong to resonator 
geometries commonly encountered, such as those with 
parallel-plane and spherically curved mirrors. 

4.3 Integral Equations for 
Resonators with Spherical Mirrors 

When the mirrors are spherical and have rectangular or 
circular apertures, the two-dimensional integral equations 
can be separated and reduced to one-dimensional equa- 
tions which are amenable to solution by either analytical 
or numerical methods. Thus, in the case of rectangular 
minors [4] [6], [24], [29], [30], the one-dimensional 
equations in Cartesian coordinates are the same as those 
for infinite-strip mirrors; for the .v coordinate, they are 

7 , <n M< ,) (xi) = J K(x 1 , xt)u™(x s )<lx* 

= f '/C(.ri, .r 2 )i< <,> (.ri)r/.r, (74) 

J-., 

where the kerne! K is given by 


E U'i, .t* 2 ) 


i/i 


nxp 


j jk 

\ ~ „ , + f/** 

l Id 


f-2.r J-. 


(75) 


Similar equations can be written for the y coordinate, so 
that £(.v, y) = u(x)v(y) and y = y,y,,. In the above equa- 
tion a i and a > are the half-widths of the mirrors in the x 
direction, d is the mirror spacing, k is 2ir/\, and X is the 
wavelength. The radii of curvature of the mirrors R \ and 
R., are contained in the factors 


(!i = 1 — 


<h 


«» 

d 

It: 


(70) 


For the case of circular mirrors [4], [31 }, [32] the equa- 
tions are reduced to the one-dimensional form by using 
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cylindrical coordinates and by assuming a sinusoidal 
azimuthal variation of the field ; that is, E(r, <f>)- Rt(r)e~^. 
The radial distribution Functions and i?j (2> satisfy the 

one-dimensional integral equations: 

\ 

7t (l) iZi <1} (n) : VrT= f ’ IUrur-W'KrlV^dr* 
i 0 

7 i (2> i2 1 { 2) (>*2) = f Kt(r u r^R t M (ri)y/r^dri (77) 
J 0 

where the kernel Ki is given by 

j l+l { n r a \ 

Ki(r 1 , r 2 ) = — —J Vn r z 

•exp + ?='•=•■)} (78) 

and Ji is a Bessel function of the first kind and /th order. 
In (77), fli and a 2 are the radii of the mirror apertures and 
d is the mirror spacing; the factors gi and g 2 are given by 
(76). 

Except for the special case of the confocal resonator 
[5] (g^g.aO), no exact analytical solution has been 
found for either (74) or (77), but approximate methods 
and numerical techniques have been employed with suc- 
cess for their solutions. Before presenting results, it is 
appropriate to discuss two important properties which 
apply in general to resonators with spherical mirrors; 
these are the properties of “equivalence” and “stability.” 

4.4 Equivalent Resonator Systems 

The equivalence properties [24], [33] of spherical- 
mirror resonators are obtained by simple algebraic manip- 
ulations of the integral equations. First, it is obvious that 
the mirrors can be interchanged without affecting the 
results; that is, the subscripts and superscripts one and 
two can be interchanged. Second, the diffraction loss and 
the intensity pattern of the mode remain invariant if both 
gi and g 2 are reversed in sign*, the eigenfunctions E and 
the eigenvalues 7 merely take on complex conjugate 
values. An example of such equivalent systems is that of 
parallel-plane (gi=g2*=l) and concentric (gi=g2= — 1) 
resonator systems. 

The third equivalence property involves the Fresnel 
number N and the stability factors Gi and G 2 , where 


Gi = 0i — 

a 2 

<? 2 =0 2 -- (79) 

ai 

If these three parameters are the same for any two resona- 
tors, then they would have the same diffraction loss, the 


same resonant frequency, and mode patterns that are 
scaled versions of each other. Thus, the equivalence rela- 
tions reduce greatly the number of calculations which are 
necessary for obtaining the solutions for the various 
resonator geometries. 

4.5 Stability Condition and Diagram 

Stability of optical resonators has been discussed in 
Section 2 in terms of geometrical optics. The stability 
condition is given by (8). In terms of the stability factors 
Gi and G 2 , it is 

0 < GiG 2 < 1 
or 

0 < Q\Qz < 1. (80) 

Resonators a r e stable if this condition is satisfied and 
unstable otherwise. 

A stability diagram [6], [24] for the various resonator 
geometries is shown in Fig. 4 where g x and g 8 are the co- 
ordinate axes and each point on the diagram represents a 
particular resonator geometry. The boundaries between 
stable and unstable (shaded) regions are determined by 
(80), which is based on geometrical optics. The fields of 
the modes in stable resonators are more .concentrated 
near the resonator axes than those in unstable resonators 
and, therefore, the diffraction losses of unstable resona- 
tors are much higher than those of stable resonators. The 
transition, which occurs near the boundaries, is gradual 
for resonators with small Fresnel numbers and more 
abrupt for those with large Fresnel numbers. The origin 
of the diagram represents the confocal system with mirrors 
of equal curvature (R x =R»=d) and is a point of lowest 
diffraction loss for a given Fresnel number. The fact that 
a system with minor deviations from the ideal confocal 
system may become unstable should be borne in mind 
when designing laser resonators. 

4.6 Modes of the Resonator 

The transverse field distributions of the rcsonutor 
modes are given by the eigenfunctions of the integral 
equations. As yet, no exact analytical solution has been 
found for the general case of arbitrary G 1 and G 2 , but 
approximate analytical expressions have been obtained to 
describe the fields in stable spherical-mirror resonators 
[5], [6], These approximate eigenfunctionsjarc the same 
as those of the optical, beam modes w' .ich are discussed in 
Section 2; that is, the field distributions are given approxi- 
mately by Hermite-Gaussian functions for rectangular 
mirrors [5], [6], [34], and by Lagucrre-Gaussian func- 
tions for circular mirrors [6], [7]. The designation of the 
resonator modes is given in Section 3.5. (The modes are 
designated as TEM mn , for rectangular mirrors and 
TEM,,r, for circular mirrors.) Figure 7 shows photo- 
graphs of some of the rectangular mode patterns of a 
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Fig. 17. Linearly polarized resonator mode configurations 
for square and circular mirrors. 



Fig. 18. Synthesis of different polarization configurations 
from the linearly polarized TEM 0 , mode. 

laser. Linearly polarized mode configurations for square 
mirrors and for circular mirrors are shown in Fig. 17. 
By combining two orthogonally polarized modes of the 
same order, it is possible to synthesize other polarization 
configurations; this is shown in Fig. 18 for the TEMoi 
mode. 

Field distributions of the resonator modes for any 
value of G could be obtained numerically by solving the 
integral equations either by the method of successive ap- 
proximations [4], [24], [3t] or by the method of kernel 
expansion [30], [32]. The former method of solution is 
equivalent to calculating the transient behavior of the 
resonator when it is excited initially with a wave of arbi- 
trary distribution. This wave is assumed to travel back 
and forth between the mirrors of the resonator, under- 
going changes from transit to transit and losing energy by 
diffraction. After many transits a quasi steady-state con- 
dition is attained where the fields for successive transits 


differ only by a constant multiplicative factor. This steady- 
state relative field distribution is then an eigenfunction of 
the integral equations and is, in fact, the field distribu- 
tion of the mode that has the lowest diffraction loss for 
the symmetry assumed (e.g., for even or odd symmetry in 
the case of infinite-strip mirrors, or for a given azimuthal 
mode index number / in the case of circular mirrors); the 
constant multiplicative factor is the eigenvalue associated 
with the eigenfunction and gives the diffraction loss and 
the phase shift of the mode. Although this simple form 
of the iterative method gives only the lower order solu- 
tions, it can, nevertheless, be modified to yield higher 
order ones [24], [35]. The method of kernel expansion, 
however, is capable of yielding both low-order and high- 
order solutions. 

Figures 19 and 20 show the relative field distributions 
of the TEMoo and TEMot modes for a resonator with a 
pair of identical, circular mirrors (N=l, ax~a» t gi=g 2 
=g) as obtained by the numerical iterative method. 
Several curves are shown for different values ofg, ranging 
from zero (confocal) through one (parallel-plane) to 1.2 
(convex, unstable). By virtue of the equivalence property 
discussed in Section 4.4, the curves are also applicable to 
resonators with their g values reversed in sign, provided 
the sign of the ordinate for the phase distribution is also 
reversed. It is seen that the field is most concentrated 
near the resonator axis for g=0 and tends to spread out 
as |g| increases. Therefore, the diffraction lo&,» is ex- 
pected to be the least for confocal resonators. 

Figure 21 shows the relative field distributions of some 
of the low order modes of a Fabry-Perot resonator with 
(parallel-plane) circular mirrors ( N= 10, ai = « 2 , gi = g 2 = i) 
as obtained by a modified numerical iterative method 
[35]. It is interesting to note that these curves arc not 
very smooth but have small wiggles on them, the number 
of which are related to the Fresnel number. These wiggles 
are entirely absent for the confocal resonator and appear 
when the resonator geometry is unstable or nearly un- 
stable. Approximate expressions for the field distribu- 
tions of the Fabry-Perot resonator modes have also been 
obtained by various analytical techniques [36], [37]. They 
are represented to first order, by sine and cosine func- 
tions for infinite-strip mirrors and by Bessel functions for 
circular mirrors. 

For the special case of the confocal resonator (g\=g* 
= 0), the eigenfunctions are self-reciprocal under the 
finite Fourier (infinite-strip mirrors) or Hankel (circular 
mirrors) transformation and exact analytical solutions 
exist [5], [38] -[40]. The eigenfunctions for infinite-strip 
mirrors are given by the prolate spheroidal wave func- 
tions and, for circular mirrors, by the generalized prolate 
spheroidal or hypcrspheroidal wave functions. For large 
Fresnel number? these functions can be closely approxi- 
mated by Hermite-Gaussian and Laperre-Gaussian 
functions which are the eigenfunction > for the beam 
modes. 
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Fig. 19. Relative field distributions of the T£M 00 mode for 
a resonator with circular mirrors (N~ 1). 



Fig. 20. Relutive field distributions of the TEM 0I mode for 
n resonator with circular mirrors (iY= I), 
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4.7 Diffraction Losses and Phase Shifts 

The diffraction loss a and the phase shift /3 for a par- 
ticular mode are important quantities in that they deter- 
mine the Q and the resonant frequency of the resonator 
for that mode. The diffraction loss is given by 

« = 1 “ I -V | 2 (81) 

which is the fractional energy lost per transit due to dif- 
fraction effects at the minors. The phase shift is given- by 

0 = angle of y (82) 

which is the phase shift suffered (or enjoyed) by the wave 
in transit from one mirror to the other, in addition to the 
geometrical phase shift which is given by 2ird/\. The 
eigenvalue y in (SI) and (82) is the appropriate 7 for the 
mode under consideration. If the total resonator loss is 
small, the Q of the resonator can be approximated by 


2 t rrf 
\cxt 


(83) 


where ct h the total resonator loss, includes losses due to 
diffraction, output coupling, absorption, scattering, and 
other effects. The resonant frequency v is given by 

I'/i’o =* (<7 + 1 ) + 0/7T (84) 


where </, the longitudinal mode order, and v 0t the funda- 
mental beat frequency, are defined in Section 3 . 5 . 
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The diffraction losses for the two lowest order (TEMoo 
and TEMoO modes of a stable resonator with a pair of 
identical, circular mirrors (a l =ao, gt=g 2 =g) are given 
in Figs. 22 and 23 as functions of the Fresnel number N 
and for various values of g. The curves are obtained by 
solving (77) numerically using the method of successive 
approximations [31]. Corresponding curves for the phase 
shifts are shown in Figs. 24 and 25. The horizontal por- 
tions of the phase shift curves can be calculated from the 
formula 

/3 = (2p + 1 + 1) arc cos 
= (2 p + 1+1) arc cos g, for (j\ - g± (85) 

which is equal to the phase shift for the beam modes 
derived in Section 3.5. It is to be noted that the loss curves 
are applicable to both positive and negative values of g 



Fig. 24. Phase shift per transit for the TliMm mode of a 
stable resonator with circular mirrors. 



Fig. 25. Phase shift per transit for the mode of a 

stable resonator with circular mirrors. 

while the phase-shift curves are for positive # only; the 
phase shift for negative g is equal to 180 degrees minus 
that for positive £. 

Analytical expressions for the diffraction loss and the 
phase shift have been obtained lor the special cases of 
parallel-plane (#=1.0) and confocal (# = 0) geometries 
when the Fresnel number is either very large (small dif- 
fraction loss) or very small (large diffraction loss) |36|, 
[38], [39], [41], [42]. In the case of the parallel-plane 
resonator with circular mirrors, the approximate expres- 
sions valid for large N, as derived by Vainshtein [36], are 


5(3/ 5) 

£V " 8h "' [(3/ + 5)‘ + 5-]- 

(8«) 


(87) 


C-18 


October 1966 / Vol. 5, No. 10 / APPLIED OPTICS 1565 


where 5 = 0.824, M**V*vN, and k p , is the (/;+l)th zero 
of the Bessel function of oiuer /. For the confocal resona- 
tor with circular mirrors, the corresponding expressions 
are [39] 


27r(87T.V)S> l+W C- 

p!(p + f+ 1) 


r[ 1+0 (i)] (88) 


0 = (2 p + l + 0 


m 


Similar expressions exist for resonators with infinite-strip 
or rectangular mirrors [36], [39]. The agreement be- 
tween the values obtained from the above formulas and 
those from numerical methods is excellent. 

The loss of the lowest order (TEMoo) mode of an 
unstable resonator is, to first order, independent of the 
mirror size or shape. The formula for the loss, which is 
based on geometrical optics, is [12] 


1 ~ V ' 1 ~ 

1 ± 1 + VI - (0>ff2)~ , 


(90) 


where the plus sign in front of the fraction applies for g 
values lying in the first and third quadrants of the stability 
diagram, and the minus sign applies in the other two quad- 
rants. Loss curves (plotted vs. N) obtained by solving the 
integral equations numerically have a ripply behavior 
which is attributable to diffraction effects [24], [43]. How- 
ever, the average values agree well with those obtained 
from (90). 


5. Concluding Remarks 

Space limitations made it necessary to concentrate the 
discussion of this article on the basic aspects of laser 
beams and resonators. It was not possible to include such 
interesting topics as perturbations of resonators, resona- 
tors with tilted mirrors, or to consider in detail the effect 
of nonlinear, saturating host media. Also omitted was a 
discussion of various resonator structures other than 
those formed of spherical mirrors, e.g., resonators with 
corner cube reflectors, resonators with output holes, or 
fiber resonators. Another important, but omitted, field is 
that of mode selection where much research work is cur- 
rently in progress. A brief survey of some of these topics 
is given in [44]. 
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APPENDIX D. PROPERTIES OF MATERIALS "ROM 100 TO 200 GHz 


1. INTRODUCTION 

Surface resistivities of metals and dissipation factors (loss tangents) 
of dielectrics cause ohmic losses in the radiometer signal channel. Metallic 
losses are predictable in both waveguides and quasi-optical assemblies 
because surface resistivities of metals and behavior are well known. How- 
ever, dielectric data in the frequency region from 100 to 200 GHz are incom- 
plete and spotty. This study shows the need for a dielectric materials 
measurement program. 

Section 2 lists the reflective losses encountered from various types 
of metals. The theoretical values are more optimistic than measured values 
by a typical factor of 3. Section 3 lists dielectric data obtained from an 
exhaustive survey with all sources referenced. In cases where 100 to 
200 GHz data are not available, data from other frequencies are listed. 
Section 4 discusses two effective methods of reducing reflective losses from 
dielectric surfaces. 


2. REFLECTIVE LOSSES FROM METALS 

Reflective losses from metals can be shown to be given by 


R 


_ 2u6 


where 

R r = fraction of energy absorbed by metal surface 

6 = skin depth of metal 

w = frequency 
c = velocity of light" 

penetrated tato b * ^ dlSt “ ce atl electromagnetic 


wave field 


6 = 


Hi 

2 


'1 + 


fe) - 


1/2 


'J 


-1 


where 


M- = magnetic permeability 
€ = dielectric permittivity 

o' = surface resistivity of metal 
If the material 


is a good conductor, cr » we and 


6 = 


cop,cr 



"" Fink > Mectronic Engineer's HandhnnL- 
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This approximation holds for frequencies beyond 1000 GHz for metals of 
typical resistivities. 

at 2 inn aD i»rq 1 ?nn t nw Ur / aC ! resistivities > s ^in depths, and reflective losses 
at 2 , 100, and 200 GHz for ten common metals used in microwave and milli- 
meter wave circuits. It should be observed that theoretical reflective losses 
trom surfaces of copper and silver are the order of 0. 1 percent. In practice 
these are approximately 0. 3 percent at 200 GHz. P 

TABLE 1. REFLECTIVE LOSSES OF COMMON METALS 


Aluminum 


Beryllium 


{66 cu 34 zn] 


Chromium 
Copper 
G 

Nickel 


Platinum 


Silver 


Assumes p r = 
Verified by i 


Skin Depth, Reflective Loss 

Resistivity (£2-m) x 10' 8 2 GHz 100 GHz 200 GHz ~~2 GHz 100 GHz 200 GHz 


0.000152 0.0012 0.0015 

0.000202 0.0014 0.0020 

0.000185 0.0013 0.0018 


0.059 


0.516 


0.203 


0.537 


0.148 


0.176 


0.042 


0.365 


0.143 


0.380 


0.000151 0.0011 0.00151 

0.000124 0.0008 0.00124 

0.000148 0.0010 0.00147 


0.000307 I 0.002 I 0.0031 


0.000120 0.0008 0.0012 


0.000319 I 0.0022 I 0.0032 


comparison with published values; see Fink Electronics Engineers' Handbook, 













































3. DIELECTRIC DISSIPATION FACTOR (LOSS TANGENT) 
Equations and Definitions 


The complex relative permittivity £ and complex index of refraction 


n are related by 


£ = (ft) 2 


where 




- J*’ 


and 


A 

n 


n - j 


c a 
4 7TV 


n(l - jK) 


g 

c = speed of light = 3x10 m/ sec 
K = extinction coefficient = co'/47rvn 
a = absorption coefficient per meter 
n = real index of refraction 

Now 


4 = 


n 


-( 


co 

4irv 


\ jnco 
/ " 2ttv 


equating components gives 
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and 


e' = 


eng 

Zttv 


tan 8 = — = 


S-rrvnco' 


(4uvn)^ - (ccy)^ 


8 = loss angle in fir ad 


When 


a « 


4-n-vn 

c 


then 


e » 



and 


l- very useful approximations 


tan 8 = 


eg 

2irvr 


In general, these approximations are made whenever a « v/c and 
fi 1 = fi Q (magnetization can be neglected for organic polymers) Then if 
absorption is assumed to be the only source of attenuation, we may write 
an expression for transmission t 


t = exp [-o-f] 


where £ is the thickness of dielectric material in question. 
Table Z lists the dielectric properties of materials. 
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TABLE 2. DIELECTRIC PROPERTIES OF MATERIALS, GHz 


Material 

Polyethylene 

Polypropylene 

Polystyrene 


Rcxolite 1 14221 
(cross linked polystyrene) 


T.P.X. (basically a 
poly 4-methyl 1-pentene) 

Teflon (treated PTFE, 
polytetrafluoroethyleno) 


PTFE (unsintered) 

RTFE I3) (sintered) 
Paraffin 

Lucite 

Plexiglas 


Mylar (polyethylene 
terephthalate] 

Fused Silica (85% density 
slip cast) 

Fused Silica (spectrosil) 


Glass (Corning 7070) 
Quartz (crystal-anisotropic) 

Ethyl Cellulose 


Dielectric 
Constant 
c ' 

Dielectric 

Loss 

c 

Loss Tangent 
c 7e tan 8 

Absorption 
Coefficient^ ) 
lnepcr-cm‘1) 

2.24 

25 

0.0005 

25 

0.00021 

25 

0.0017 

25 

2.31 

143 







2.31 

343 

0.00048 

300 

0.00021 

300 

0.02 

300 

2.33 

850 

0.0009 

850 

0.0004 

850 

0.1087 

850 

2.25 

35 







2.25* 

120 

0. 1012 

120 

0.00053 

120 

0.02 

120 

2.54 

25 

0.0u3 

25 

0.0012 

25 

0.010 

25 

2.52 

140 

0.005 

140 

0.002 

140 

0.0924 

140 

2.53 

210 

0.00316 

120 

0.0013 

120 

0.05 

120 

2.57 

343 

0.0063 

168 

0.0025 

16B 

0.14 

168 

2.52 

850 

0.00227 

850 

0.0009 

850 

0.25 

850 

2.54 

10 

0.00127 

10 

0.0005 

10 

0.0017 

10 

2.47 

140 

0.00494 

140 

0.002 

140 

0.0922 

140 

2.50 

210 







2.54 

343 







2.525 

850 

0.0076 

850 

0.003 

850 

0.85 

850 

2.126 

35 

0.001 

35 

0.00048 

35 

0.0051 

35 

2.15* 

120 

0.0016 

120 

0.00076 

120 

0.028 

120 

2.08 

25 

0.00125 

25 

0.0006 

25 

0.0045 

25 

2.058 

35 

0.00062 

50 

0.0003 

50 

0.0045 

50 

2.039 

60-*90 







2 05 

140 

0.00615 

140 

0.003 

140 

0.126 

140 

2.38 

210 







2.07 

343 







2.042 

850 

0.00143 

850 

0.0007 

850 

0.95 

850 

1.950 

35 







1.952 

35 

0.0001 

35 

0.00005 

35 

0.00049 

35 

2.052 

300 

0.0002 

300 

0.0001 

300 

0.009(2) 

300 

1.933 

300 

0.00084 

300 

0.000435 

300 

0.38 

300 

2.2 

25 

0.00066 

25 

<0.0003 

25 

<0.0023 

25 

2.19 

120 

0 00595 

120 

0.0027 

120 

0.10 

120 

2.2 

168 

0.00296 

168 

0.00135 

168 

0.07 

168 

2.3 

2-300 



0.001 




2.57 

25 

0,0082 

25 

0.0032 

25 

0.027 

25 

2.56 

140 

0.038 

140 

0.015 

140 

0.696 

140 

2.50 

210 







2.59 

25 

0.017 

25 

0.0067 

25 

0.056 

25 

2.59 

120 

0.019 

120 

0.0075 

120 

0.30 

120 

2.60 

143 

0,023 

168 

O.OOB9 

16B 

0,49 

16B 

2.61 

343 







3.35 

140 

0.0335 

140 

0.01 

140 

0.537 

140 

3,145 

55 

0.0133 

55 

0.0044 

55 

0.C9 

55 

3.39 

94 

0,0085 

94 

0.0026 

94 

0.09.13 

94 

3,78 

10 

0.00064 

10 

0.00017 

10 

0.0007 

10 

3.85 

2-300 

0.00039 


0.0001 


0.003. 

75 

3.82 

60-90 







3.9 

25 

0,0121 

25 

0.0031 

25 

0.032 

25 

4.0 

2-300 

0.0095 


0.0024 


_ 


4.43W 

35 

0.00014 

35 

1 0.000031 35 

0.00048 

35 

4.45(4 

2-300 



0.0002 


— 


4.64(5 

35 







2.65 

25 

0,0795 

25 

0.03 

25 

0.256 

25 

3.71 

140 

0.371 

140 

0.1 

140 

5.65 

140 


Reference 


b 

b and j 
h 


i 

a 

a and d 
b and d 
h 

i 


e and k 


d 

b and d 
b 


Note: Temperature is usually at or near 25°C or 3Q0°K. 

'Assumed valt e necessary to calculate absorption. ..... . • ,. . 

1) Absorption coefficient per unit length is without regard to effects at the medium boundaries. The neper here, as defined in submillimetre optics, is only 
half the magnitude of the neper used in E.E.; therefore, 1 Neper ■ 4.343 dB is the conversion. 

2) This is an upper bound corresponding to the limit of the particular experimental technique used in reference j. 

3) "Sintered" PTFE is obtained by stretching a sample at 250°C and allowing it to cool under tension. 

4) Ordinary ray Co- 
51 Extraordinary e e . 

6) Arbitrarily chosen frequency between 2 and 300 GHz. 
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Table 2 (continued) 
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4. DIELECTRIC REFLECTION LOSSES 


Reflective losses from a dielectric surface are determined by 


L 

r 


n - 1 
n + 1 


2 


where n is the refractive index of the material or n/c. Typical losses for an 
air -quartz surface are the order of 1 dB. A quartz lens having two surfaces 
will have loss contributions from both surfaces or 2 dB. 

Reflective losses can be minimized by two methods (see Figure D-l). 
A common method in optics is a 1/4 coating on the surface of proper refrac- 
tive index. For example, if a fused silica lens n s = 1. 962 and teflon 
n t ~ 1.43, the reflective losses of the quartz can be reduced by coating the 
quartz with a layer of teflon. The reflective loss with coating is then 


L 


(R 10 + R 21 


- 4R 1Q R 21 


1 + R 10 R 21 


where 


R 


10 


n t - 1 
n t + 1 


and 




n + n, 
s t 


which reduces the loss to 0. 019, or 0. 08 dB per surface. 

Another method of eliminating surface reflections is to machine 
grooves in the optical surface to simulate the quarter wave spacing. This 
technique applies especially to materials like Rexolite which mahine easily, 
but can also be achieved with Teflon. The technique is to adjust the ratio of 
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INCIDENT RAY ANGLES (0) AT INTERFACE OF ANY 
SPHERICAL PLANO-CONVEX LENS 
(ASSUME PARAXIAL COLLIMATED BEAM) 


MATCHED GROOVE 

FIGURE D-1. REDUCTION OF REFLECTION LOSSES BY DIELECTRIC FILM 
AND GROOVE MATCHINGS 


i'i 

groove width to groove spacing to achieve the correct "effective" matching 
index of refraction, and to adjust the depth of the grooves to obtain the 
quarter wave spacing. * 


nn> 

*S. B. Cohn and T. Morita, "Measured Performance of Dielectric Lenses, 
IRE Transactions on Antennas and Propagation, January 1956, pp. 3 1-33. 
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